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Objectives 


>  Experimental  visualization  of  the  stress  and  deformation  fields  near  dynamieally 
propagating  crack  tips  subjected  to  thermal,  mechanical  and  thermo-mechanical  loading 
in  graded  materials. 

>  Development  of  rate-dependent  constitutive  models  for  graded  multifunctional  materials 
under  thermo-mechanical  loading. 

>  Identification  of  damage  mechanics  in  graded  multifunctional  materials  subjected  to  high 
strain  rate  loading. 

>  Development  of  asymptotic  stress  fields  for  stationary  and  propagating  cracks  subjected 
to  combined  thermo-mechanical  loading  in  graded  materials. 

>  A  critical  evaluation  of  the  influence  of  material  gradation,  thermal  loading  and  loading 
rates  on  the  damage  mechanisms  of  graded  materials. 

ACCOMPLISHMENTS 

A  project  entitled  ''Dynamic  Failure  of  Integrated  Durable  Hot  Structure  for  Space  Access 
Vehicles'"  was  funded  by  the  AFOSR  for  a  period  of  3  years  (Jan  2006  to  Nov  2008).  During  this  project  a 
fundamental  study  was  conducted  to  evaluate  dynamic  fracture  and  failure  in  functionally  graded 
materials  at  room  and  at  high  temperatures.  Room  temperature  studies  included  experimental  and 
analytical  evaluation  of  dynamic  fracture  in  model  materials.  Analytical  studies  were  also  conducted  to 
develop  thermo-mechanical  stress,  strain  and  displacement  fields  around  stationary  cracks  and  around 
cracks  propagating  under  steady  state  conditions.  We  also  initiated  preliminary  studies  on  dynamic 
constitutive  behavior  as  well  as  fracture-initiation  toughness  of  a  realistic  FGM  (Ti/TiB)  at  various 
temperatures  and  high  strain  rates  of  loading. 

This  research  work  produced  (8)  journal  and  (6)  proceedings  papers.  One  of  the  papers  received 
the  “Best  Paper  Award”  from  the  Journal  of  Strain  and  also  received  the  Flyde  Prize  from  the  British 
Society  of  Strain  Measurement.  Some  of  the  details  of  the  work  conducted  on  this  project  are  given 
below. 

1.  Asymptotic  Stress  Fields  for  Thermo-mechanically  Loaded  Cracks  in  FGMs 

The  problem  of  a  stationary  crack  in  functionally  graded  materials  (FGM)  subjected  to  a 
combination  of  thermal  and  mechanical  loading  was  considered.  An  asymptotic  analysis  coupled  with 
Westergaard’s  stress  function  approach  was  used  to  characterize  the  stress  field  around  the  crack-tip. 
Thermal  and  mechanical  properties  (e.g.  elastic  modulus,  coefficient  of  thermal  expansion  and  thermal 
conductivity)  were  assumed  to  vary  exponentially.  The  crack  was  assumed  to  be  inclined  to  the  direction 
of  the  property  gradation.  The  thermal  loading  was  taken  to  be  a  uniform  heat  flow  in  a  direction  inclined 
to  the  crack.  The  principal  of  superposition  from  linear  elasticity  was  used  to  solve  the  problem,  whereby 
the  problem  was  divided  into  a  number  of  sub-problems.  The  first  four  terms  in  the  expansion  of  the 
stress  field  were  derived  to  explicitly  bring  out  the  influence  of  nonhomogeneity  on  the  structure  of  the 
stress  field.  It  was  observed  that  the  presence  of  heat  flow  produced  no  additional  singularity  and  hence 
the  classical  inverse  square  root  singularity  still  prevails  around  the  crack  tip.  Using  these  stress  fields, 
contours  of  constant  maximum  shear  stress  were  generated  and  the  effect  of  thermal  loading  on  the  crack- 
tip  stress  field  was  discussed. 

2.  Asymptotic  Analysis  and  Reflection  Photoelasticity  to  Study  Transient  Crack  Propagation  in 
FGMs 


The  behavior  of  a  rapidly  moving  transient  crack  in  functionally  graded  materials  (FGMs)  was 
investigated  theoretically  and  experimentally.  First,  a  systematic  theoretical  analysis  was  presented  for  the 
development  of  the  transient  elastodynamic  local  stress,  strain  and  displacement  field  expansions  near  a 
growing  mixed  mode  crack  tip  in  FGMs.  The  crack  propagation  direction  was  assumed  to  be  inclined  to 
the  direction  of  the  property  variation.  The  displacement  potential  approach  in  conjunction  with 
asymptotic  analysis  was  utilized  to  derive  explicit  expressions  for  stress,  strain  and  in-plane  displacement 
fields.  To  further  understand  the  transient  crack  growth  behavior  a  series  of  dynamic  fracture  experiments 
were  performed  with  in-house  fabricated  functionally  graded  material.  The  full-field  stress  data  around 
the  crack  was  recorded  using  dynamic  photoelasticity  and  high-speed  digital  photography.  Due  to 
opaqueness  of  FGMs,  birefringent  coatings  were  employed  to  obtain  the  full-field  isochromatics  around 
the  crack-tip.  The  stress  field  expansions  developed  in  the  first  part  of  the  study  were  used  to  interpret  the 
experimental  observations.  The  results  of  the  experiments  showed  that  the  higher  order  transient 
expansion  provides  an  accurate  representation  of  crack-tip  fields  under  severe  transient  conditions. 

3.  Dynamic  Crack  Tip  Stress  and  Displacement  Fields  under  Thermo-Mechanical  Loading 

Thermomechanical  stress  and  displacement  fields  for  a  propagating  crack  in  functionally 
graded  materials  (FGMs)  were  developed  using  displacement  potentials  and  asymptotic  analysis. 
The  shear  modulus,  mass  density,  and  coefficient  of  thermal  expansion  of  the  FGMs  were 
assumed  to  vary  exponentially  along  the  gradation  direction.  Temperature  and  heat  flux 
distribution  fields  were  also  derived  for  an  exponential  variation  of  thermal  conductivity.  The 
mode  mixity  due  to  mixed-mode  loading  conditions  around  the  crack  tip  was  accommodated  in 
the  analysis  through  the  superposition  of  opening  and  shear  modes.  Using  the  asymptotic  stress 
fields,  contours  of  constant  maximum  shear  stress  were  developed  and  the  results  were  discussed 
for  various  crack-tip  thermo-mechanical  loading  conditions. 

4.  Dynamic  Constitutive  Behavior  of  Ti/TiB  layer  FGM  under  Thermo-mechanical  Loading 

An  experimental  investigation  was  conducted  to  evaluate  the  thermo-mechanical 
constitutive  behavior  of  a  functionally  graded  material  (FGM)  under  dynamic  loading. 
Cylindrical  specimens  were  machined  from  titanium  /  titanium  mono-boride  (Ti/TiB)  layered 
FGM  plate  using  electrical  discharge  machining  (EDM).  A  Split  Hopkinson  Pressure  Bar 
(SHPB)  apparatus  with  infrared  spot  heaters  was  used  to  investigate  the  effect  of  temperature  on 
mechanical  response  of  the  FGM  material.  A  series  of  experiments  were  conducted  at  different 
temperatures  and  the  stress  strain  relation  for  different  temperatures  was  obtained.  The  material 
showed  high  thermal  softening  at  elevated  temperature  resulting  in  a  reduction  in  compressive 
strength  and  an  increase  in  failure  strain. 

5.  Effect  of  Dynamic  Thermo-Mechanical  Stress  Field  on  Crack  Propagation  Direction  in 
Functionally  Graded  Materials 

Asymptotic  analysis  in  conjunction  with  displacement  potentials  was  used  to  develop 
thermo-mechanical  stress  fields  for  a  mixed  mode  propagating  crack  in  a  functionally  graded 
material  (FGM).  First,  asymptotic  temperature  fields  were  derived  for  an  exponential  variation  of 
thermal  conductivity  and  later  these  temperature  fields  were  used  in  deriving  stress  fields.  Using 
asymptotic  thermo-mechanical  stress  fields  the  variation  of  maximum  shear  stress, 
circumferential  stress  and  strain  energy  density  as  a  function  of  temperature  around  the  crack-tip 
were  generated.  Finally,  utilizing  the  minimum  strain  energy  density  criterion  and  the  maximum 


circumferential  stress  eriterion,  the  eraek  growth  direction  for  various  crack-tip  speeds,  non¬ 
homogeneity  eoeffieients  and  temperature  fields  were  determined. 

6.  Dynamic  Fracture  Initiation  Toughness  of  Ti/TiB  Layered  FGM  under  Thermo-mechanical 
Loading 

Initial  experiments  were  conducted  to  study  dynamic  fracture  initiation  toughness  of  Ti/TiB 
layered  FGM  was  investigated  using  the  modified  SHPB  apparatus.  A  three-point  bending  test  specimen 
was  machined  from  Ti/TiB  plate  using  EDM  technique.  The  crack  was  made  using  EDM  wire  of  0.002” 
which  results  in  a  very  small  notch  radius.  For  the  elevated  temperature  experiment,  the  specimen  was 
first  heated  up  using  an  induction  heating  system  and  the  incident  bar  was  made  in  contact  manually  just 
before  propelling  the  projectile.  A  series  of  experiment  were  conducted  at  a  fixed  high  loading  rate  and 
different  temperatures  to  study  the  effect  of  temperature  and  loading  rate  on  the  fracture  initiation 
toughness. 

7.  Experimental  Facilities  Developed 

During  the  grant,  facilities  for  high  temperature  testing  of  FGMs  were  established.  These 
included  set-ups  with  infrared  heating  system  and  induction  coil  heating  system  that  were  coupled  to  our 
existing  SHPB  apparatus.  Also,  two  and  three  dimensional  full  field  digital  image  correlation  (DIG) 
techniques  including  three  high  speed  cameras  (Photron  SAl)  have  been  procured. 
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ABSTRACT 


A  fundamental  study  is  conducted  to  evaluate  dynamic  fracture  and  failure  in 
functionally  graded  materials  (FGM)  at  room  and  high  temperatures.  The  study 
includes  both  analytical  and  experimental  evaluation  of  dynamic  fracture  in  model 
and  commercially  available  materials. 

Asymptotic  analysis  in  conjunction  with  displacement  potentials  is  used  to 
develop  thermo-mechanical  stress  helds  for  a  mixed  mode  propagating  crack  in 
homogenous  and  FGMs.  First,  asymptotic  temperature  helds  are  derived  for  an 
exponential  variation  of  thermal  conductivity  and  later  these  temperature  helds 
are  used  in  deriving  stress  helds.  Using  asymptotic  thermo-mechanical  stress  helds 
the  variation  of  maximum  shear  stress,  circumferential  stress  and  strain  energy 
density  as  a  function  of  temperature  around  the  crack-tip  are  generated.  Finally, 
utilizing  the  minimum  strain  energy  density  criterion  and  the  maximum  circum¬ 
ferential  stress  criterion,  the  crack  growth  direction  for  various  crack-tip  speeds, 
non-homogeneity  coefficients  and  temperature  helds  are  determined. 

An  experimental  investigation  is  conducted  to  evaluate  the  thermo-mechanical 
constitutive  behavior  of  a  FGM  under  dynamic  loading.  Gylindrical  specimens  are 
machined  from  titanium  /  titanium  mono-boride  (Ti/TiB)  layered  FGM  plate  us¬ 
ing  electrical  discharge  machining  (EDM).  A  split  Hopkinson  pressure  bar  (SHPB) 
apparatus  with  infrared  spot  heaters  is  used  to  investigate  the  ehect  of  tempera¬ 
ture  on  mechanical  response  of  the  FGM.  A  series  of  experiments  are  conducted 
at  diherent  temperatures  and  the  stress  strain  relation  for  diherent  temperatures 
is  obtained.  The  material  showed  high  thermal  softening  at  elevated  temperature 
resulting  in  a  reduction  in  compressive  strength  and  an  increase  in  failure  strain. 

Quasi-static  and  dynamic  fracture  initiation  toughness  of  Ti/TiB  layered  FGM 
is  investigated  using  a  three  point  bend  specimen.  The  modihed  SHPB  apparatus 


in  conjunction  with  induction  coil  heating  system  is  used  during  elevated  temper¬ 
ature  dynamic  loading  experiments.  A  simple  and  accurate  technique  has  been 
developed  to  identify  the  time  corresponding  to  the  load  at  which  the  fracture 
initiates.  A  series  of  experiments  are  conducted  at  different  temperatures  ranging 
from  room  temperature  to  800  °C,  and  the  effect  of  temperature  and  loading  rate 
on  the  fracture  initiation  toughness  is  investigated.  The  material  fracture  tough¬ 
ness  is  found  to  be  sensitive  to  temperature  and  the  fracture  initiation  toughness 
increases  as  the  temperature  increases.  Furthermore,  the  fracture  initiation  tough¬ 
ness  is  strain  rate  sensitive  and  is  higher  for  dynamic  loading  as  compared  to 
quasi-static  loading. 

A  model  transparent  graded  materials  is  used  to  investigate  the  steady  state 
and  transient  crack  propagation  in  a  functionally  graded  material.  High-speed 
digital  photography  combined  with  photoelasticity  technique  is  used  to  record  the 
full-held  stress  data  around  the  propagating  crack.  By  analyzing  the  photoelastic 
fringe  patterns,  the  propagation  crack  tip  velocity  and  displacement  are  obtained. 
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PREFACE 


The  present  study  addresses  the  dynamic  fracture  and  failure  of  functionally 
graded  materials  (FGMs)  under  thermo  -  mechanical  loading.  Analytical  and 
experimental  investigation  on  the  model  and  commercially  available  FGMs  are 
presented. 

The  hrst  chapter  provides  a  detailed  summary  of  a  relevant  review  of  litera¬ 
tures  followed  by  a  brief  introduction  and  scope  of  the  present  study. 

Ghapter  2  presents  analytical  solution  for  homogenous  material  subjected  to 
thermo-mechanical  loading.  Thermo-mechanical  stress  held  equations  are  devel¬ 
oped  for  a  mixed-mode  crack  propagating  at  constant  velocity  in  homogeneous  and 
isotropic  materials  using  asymptotic  approach  along  with  displacement  potentials. 
Using  these  thermo-mechanical  stress  helds,  various  components  of  the  stresses  are 
developed  and  the  effects  of  temperature  on  these  stress  components  are  discussed. 

Asymptotic  analysis  in  conjunction  with  displacement  potentials  has  been 
used  to  develop  thermo-mechanical  stress  helds  for  a  mixed  mode  propagating 
crack  in  a  FGMs  and  discussed  in  detail  in  chapter  3.  Using  asymptotic  thermo¬ 
mechanical  stress  helds  the  variation  of  maximum  shear  stress,  circumferential 
stress  and  strain  energy  density  as  a  function  of  temperature  around  the  crack- 
tip  are  generated. Utilizing  the  minimum  strain  energy  density  criterion  and  the 
maximum  circumferential  stress  criterion,  the  crack  growth  direction  for  various 
crack-tip  speeds,  non-homogeneity  coefficients  and  temperature  helds  are  deter¬ 
mined. 

Ghapter  4  discusses  an  experimental  investigation  conducted  to  evaluate  the 
thermo-mechanical  constitutive  behavior  of  a  FGMs  under  dynamic  loading.  A 
split  Hopkinson  pressure  bar  (SHPB)  apparatus  with  infrared  spot  heaters  is  used 
to  investigate  the  ehect  of  temperature  on  mechanical  response  of  the  FGM.  The 


V 


material  showed  high  thermal  softening  at  elevated  temperature  resulting  in  a 
reduction  in  compressive  strength  and  an  increase  in  failure  strain. 

Chapter  5  presents  the  quasi-static  and  dynamic  fracture  initiation  tough¬ 
ness  of  Ti/TiB  layered  FGM.  The  modihed  SHPB  apparatus  in  conjunction  with 
induction  coil  heating  system  is  used  during  elevated  temperature  dynamic  load¬ 
ing  experiments.  A  simple  and  accurate  technique  has  been  developed  to  identify 
the  time  corresponding  to  the  load  at  which  the  fracture  initiates.  The  material 
fracture  toughness  is  found  to  be  sensitive  to  temperature  and  the  fracture  initia¬ 
tion  toughness  increases  as  the  temperature  increases.  Furthermore,  the  fracture 
initiation  toughness  is  strain  rate  sensitive  and  is  higher  for  dynamic  loading  as 
compared  to  quasi-static  loading. 

Chapter  6  discusses  an  experimental  investigation  of  the  steady  state  and  tran¬ 
sient  crack  propagation  in  a  transparent  functionally  graded  material.  High-speed 
digital  photography  combined  with  photoelasticity  technique  is  used  to  record  the 
full-held  stress  data  around  the  propagating  crack.  By  analyzing  the  photoelastic 
fringe  patterns  the  dynamic  stress  intensity  factor  and  propagation  velocity  of  the 
crack  tip  are  obtained. 

Finally,  the  summary  of  the  present  work  and  recommendation  for  the  future 
work  are  presented  in  chapter  7. 
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CHAPTER  1 


Introduction 

1.1  Introduction 

The  way  to  outer  space  requires  a  high  level  of  technology  to  design  safe,  reli¬ 
able,  economical  and  novel  transportation  systems.  When  developing  a  new  design, 
the  engineers  consider  the  conceivable  failures  of  the  new  system.  Today  the  main 
challenging  technical  problem  in  designing  space  access  vehicles  is  the  development 
of  the  reentry  heat  shield  that  can  be  reused  with  minimal  refurbishment  and  test¬ 
ing.  On  reentry,  a  shock  wave  forms  ahead  of  the  spacecraft,  which  heats  the 
compressed  air  in  the  region  near  the  spacecraft  to  as  much  as  20,000  °C,  creating 
thermal  loads  in  excess  of  one  megawatt  per  square  meter.  This  enormous  ther¬ 
mal  load  heats  the  aircraft  structure  up  to  approximately  1,700  °C.  This  requires 
an  external  thermal  protection  systems  (TPS)  which  usually  in  the  form  of  rigid 
surfaces  in  areas  of  high  or  moderate  working  temperature.  In  addition  to  high 
temperature  and  high  convective  fluxes,  TPS  must  be  able  to  withstand  mechani¬ 
cal  stresses  associated  with  vibrations  at  launch,  acoustic-frequency  vibrations  and 
structural  movement  of  the  vehicle,  as  well  as  landing  impact.  Usually  TPS  are 
composite  layers  of  two  different  materials  with  ceramic  on  the  top  due  to  its  high 
thermal  properties  and  metals  in  the  other  end  for  structural  integrity.  However, 
the  property  mismatch  and  high  thermal  gradient  across  these  two  materials  re¬ 
sults  in  delamination  and  interface  failure.  A  very  recent  unfortunate  example  of 
this  type  of  failure  was  the  Space  Shuttle  Columbia.  All  of  the  materials  involved 
in  the  Columbia  accident  were  classihed  as  belonging  to  the  TPS. 

In  recent  years  a  new  system  called  functionally  graded  material  (PGM)  has 
been  developed  to  circumvent  for  the  interface  problems  that  occurs  in  bi-material 
systems.  FGMs  are  materials,  which  have  continuously  or  discreetly  varying  me- 
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chanical  and  thermal  properties  in  a  specific  direction.  This  type  of  material  was 
first  proposed  in  Japan  in  1987  for  space  access  vehicles  as  an  alternative  to  thermal 
barrier  coatings  [1].  FGMs  are  potential  candidates  to  be  used  as  integrated  hot 
structure  in  space  access  vehicles  as  the  material  combination  and  the  composition 
of  gradation  in  these  materials  can  be  tailored  to  optimize  their  performance  to 
meet  multiple  functions  like  thermal  resistance  and  structural  integrity.  Typically, 
FGMs  are  made  of  a  metal  and  a  ceramic  as  opposite  faces  with  the  intermediate 
zones  consisting  of  varying  volume  fractions  of  constituents.  Titanium  /  titanium 
mono-boride  (Ti/TiB),  FGM  which  is  currently  in  use,  consists  of  the  desirable 
properties  of  ceramic,  such  as  hardness,  corrosion  resistance,  and  high  melting  tem¬ 
perature,  without  losing  the  required  properties  of  metallic  titanium,  such  as  good 
fracture  toughness,  machinability  and  weldability.  For  efficient  design  of  structures 
using  FGM  a  fundamental  understanding  of  the  properties  of  the  material  at  dif¬ 
ferent  loading  conditions  is  required.  For  designing  a  failure  criteria  the  fracture 
parameters  and  constitutive  behavior  of  this  material  is  important. 

1.2  Related  Literatures 

The  elasticity  of  non-homogeneous  solids  had  been  an  interesting  subject  well 
before  FGMS  are  developed.  Solutions  to  different  types  of  mixed  boundary  value 
problems  in  non-homogeneous  solids  have  been  obtained  by  Plevako,  [2],  Kas- 
sir  [3,  4],  Popov  [5]  and  Singh  [6,  7].  The  earliest  study  into  the  crack  problems 
in  non-homogenous  solids  can  be  traced  down  to  early  1970s.  Kassir  [3]  investi¬ 
gated  the  problem  of  an  axisymmetric  penny  shaped  crack  in  a  cylinder  subjected 
to  pure  torsion  by  varying  the  shear  modulus  in  axial  direction.  Later,  Dhaliwal 
and  Singh  [8]  obtained  an  expression  for  the  stress  intensity  factor  for  a  crack  in  a 
semi-infinite  non-homogeneous  solid  subjected  to  crack  face  anti  plane  shear  load¬ 
ing.  Glements  et  al.  [9]  also  solved  antiplane  crack  problems  for  inhomogeneous 
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materials  and  obtained  the  energy  release  rates  for  various  crack  orientations  and 
property  variations.  These  initial  investigations  were  focussed  on  the  tearing  mode 
of  fracture,  because  of  its  simplicity  among  the  three  modes  of  fracture.  Gerasoulis 
and  Srivastav  [10]  are  the  first  to  analysis  the  opening  mode  fracture  in  nonho- 
mogeneous  solids,  which  is  of  the  more  practical  importance,  by  assuming  plane 
strain  conditions.  Later  on  Delale  and  Erdogan  [11],  Eischen  [12]  and  Erdogan  et 
al.  [13]  solved  crack  problems  for  non-homogeneous  materials  under  quasi-static 
loading. 

With  the  introduction  of  FGMs,  research  on  fracture  mechanics  of  nonhomo- 
geneous  solids  gained  additional  impetus.  Jin  and  Noda  [14],  Konda  and  Erdo¬ 
gan  [15]  and  Erdogan  [16]  obtained  the  stress  intensity  factors  for  cracks  in  FGMs 
for  different  crack  orientation,  types  of  loading  and  levels  of  nonhomogeneity.  Jin 
and  Batra  [17,  18],  studied  the  fracture  toughness  and  resistance  curve  (R-curve) 
of  metal-ceramic  FGMs  using  the  crack  bridging  concept  and  rule  of  mixtures  and 
showed  that  the  fracture  toughness  is  signihcantly  increased  when  a  crack  grows 
from  the  ceramic-rich  region  to  the  metal-rich  region  in  the  FGM.  All  these  in¬ 
vestigations  concluded  that  the  inverse-square  root  singularity  at  the  crack  is  not 
affected  by  nonhomogeneity. 

All  the  above  investigations  are  limited  to  quasi-static  loading  conditions.  For 
propagating  cracks  in  FGMs,  Nakagaki  et  al.  [19]  developed  a  hnite  element  simu¬ 
lation  of  fracture  occurring  in  a  ceramic-metal  FGM  under  a  dynamic  load.  They 
observed  that  the  crack  opening  displacement  (GOD)  was  higher  in  the  FGM  when 
crack  propagated  from  the  ceramic  rich  side  to  the  metal  rich  side  as  compared 
to  the  reverse  direction  of  propagation.  Later  on  Parameswaran  and  Shukla  [20] 
and  Ghalivendra  et  al.  [21]  developed  the  structure  of  the  first  stress  invariant 
and  the  out  of  plane  displacement.  In  their  study  they  brought  out  the  effects 


3 


of  nonhomogeniety  through  an  asymptotic  analysis.  Rousseau  and  Tippur  [22] 
have  investigated  the  dynamic  fracture  of  continuously  graded  particulate  (glass  - 
epoxy)  FGMs  for  cracks  aligned  along  the  gradient  using  interferometry.  Lee  [23] 
developed  nonhomogeneity  specihc  terms  for  individual  stress  and  displacement 
components  using  displacement  potentials.  Recently,  Shukla  and  Jain  [24]  and 
Chalivendra  and  Shukla  [25]  developed  transient  held  equations  for  cracks  propa¬ 
gating  at  arbitrary  velocities.  Chalivendra  [26]  developed  an  asymptotic  analysis 
of  the  transient  out  of  plane  displacement  helds  for  a  curved  crack  propagating  at 
arbitrary  velocity  in  FGMs.  Kim  and  Paulino  [27]  evaluated  the  mixed-mode  stress 
intensity  in  the  form  of  an  equivalent  domain  integral,  in  combination  with  the  h- 
nite  element  method.  In  another  numerical  study,  Jin  and  Dodds  [28]  simulated 
the  crack  growth  resistance  in  a  ceramic/metal  FGM  using  a  cohesive  zone  ahead 
of  the  crack  front.  A  review  paper  by  Shukla  et  ah  [29]  presents  a  comprehensive 
summary  of  dynamic  fracture  studies  in  FGM. 

There  are  very  few  studies  on  the  crack  problems  of  FGM  under  thermal  or 
thermo- mechanical  loading.  Kawasaki  and  Watanabe  [30],  evaluated  the  thermo¬ 
mechanical  properties  of  metal/ceramic  which  simulated  the  real  environment  of 
the  heated  inner  wall  of  a  rocket  combustor.  They  used  Disk-shaped  graded  sam¬ 
ples  of  a  material  combination  of  partially  stabilized  zirconia  and  stainless  steel,  in 
which  the  ceramic  surface  of  the  sample  was  heated  with  a  burner  flame  and  the 
back  surface  was  cooled  with  flowing  water.  They  determined  the  critical  tempera¬ 
ture  of  the  hrst  crack  formation,  which  was  always  observed  on  the  ceramic  surface 
during  cooling,  and  they  found  that  it  is  almost  constant  regardless  of  specimen 
size  and  compositional  prohle.  Based  on  the  experiment,  they  concluded  that  the 
critical  temperature  is  material-dependent.  They  also  discussed  the  mechanism 
of  crack  formation  and  deflection,  including  spallation,  on  the  basis  of  the  stress 
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distributions  in  the  specimen  during  testing  and  the  fracture  mechanics  approach. 
According  to  their  study,  FGM  structure  is  advantageous  for  TBC  application 
with  the  particular  emphasis  on  the  capability  of  vertical-crack  arrest  and  the  pos¬ 
sibility  of  its  propagation  control  by  tailoring  the  graded  structure  on  the  basis  of 
fracture  mechanics.  Noda  [31]  studied  the  thermal  stresses  and  thermal  stress  in¬ 
tensity  factor  (SIF)  in  the  FGM  subjected  to  steady  temperature  fields  or  thermal 
shock.  He  also  presented  analytical  and  numerical  studies  on  the  optimal  compo¬ 
sition  prohle  problems  of  the  FGM  in  decreasing  thermal  stresses.  Based  on  his 
study  he  forwarded  the  following  remarks:  1.  When  the  continuously  changing 
composition  between  ceramics  and  metals  selected  properly,  thermal  stresses  and 
the  thermal  SIFin  FGM  can  be  drastically  decreased.  2.  The  crack  propagation 
path  can  be  predicted  by  fracture  mechanics  analysis. Kawasaki  and  Watanabe  [32], 
conducted  a  well  controlled  burner-heating-test  in  order  to  study  the  thermal  frac¬ 
ture  behavior,  the  cyclic  thermal  fracture  behavior  and  spallation  life  of  different 
metal/ceramic  functionally  graded  thermal  barrier  coatings  (TBG).  In  PSZ/INIOO 
FGM,  the  fracture  toughness  improved  with  an  increase  in  the  metal  phase  content. 
In  PSZ/Inco718  FGM  the  fracture  toughness  was  lower  than  that  of  PSZ/InlOO 
FGM,  owing  to  roughly  dispersed  metal  particles.  They  observed  that  function¬ 
ally  graded  TBG  possess  the  desirable  effect  for  improvement  of  spallation  life 
under  cyclic  thermal  loads.  In  addition,  they  observed  that  the  spallation  life  in 
functionally  graded  coatings  depends  on  the  composition  prohle. 

All  the  analytical  studies  presented  in  the  previous  section  provide  a  closed 
form  solution  for  determining  SIF,  however  for  extracting  fracture  parameters  from 
experimental  studies  under  thermo-mechanical  loading  an  asymptotic  stress  held 
around  the  crack  tip  is  required.  In  this  aspect  recently,  Jain  et  al.  [33]  developed 
a  held  equation  for  a  stationary  crack  in  FGMs,  subjected  to  thermo-mechanical 
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loading.  They  used  an  asymptotic  analysis  coupled  with  Westergaard’s  stress 
function  approach  to  characterize  the  stress  field  around  the  crack  tip.  Using 
the  stress  helds  they  generated  the  contour  of  constant  maximum  shear  stress  and 
discussed  the  effect  of  thermal  loading  on  the-crack  tip  stress  fields.  They  observed 
that  the  presence  of  a  thermal  load  produces  no  additional  singularity  and  hence 
the  classical  inverse  square  root  singularity  still  prevails  around  the  crack  tip.  The 
development  of  an  asymptotic  expansion  held  equation  for  FGMs  under  thermo¬ 
mechanical  dynamic  loading  is  a  problem  yet  to  be  investigated. 

In  other  directions,  there  are  several  studies  related  to  the  stress  and  displace¬ 
ment  response  of  FGMs  under  different  loading  conditions.  For  instance  Wang  et 
al.  [34]  present  a  solution  for  the  displacement  and  stress  in  an  FGM  subjected  to 
a  vertical  point  load  in  a  continuously  nonhomogeneous  transversely  isotropic  half¬ 
space  with  Young’s  and  shear  moduli  varying  exponentially  with  depth.  Horgan 
and  Ghan  [35]  investigate  the  effects  of  material  nonhomogeneity  on  the  response  of 
linearly  elastic  isotropic  hollow  circular  cylinders  under  uniform  internal  or  external 
pressure.  Li  et  al.  [36]  examined  the  mechanical  behavior  of  layered  plates  made  of 
metal-ceramic  composites  with  the  volume  fraction  of  ceramic  reinforcement  vary¬ 
ing  through  thickness  direction  under  impulse  loading.  Ghi  and  Ghung  [37]  studies 
the  mechanical  behavior  of  a  functionally  graded  material  plate  under  transverse 
loading. 

There  are  also  a  few  theoretical  studies  on  displacement  and  stress  response 
of  FGMs  under  thermo-mechanical  loading.  Praveen  and  Reddy  [38]  investigated 
the  response  of  functionally  graded  ceramic-metal  plates  under  mechanical  and 
thermal  loading  using  a  plate  hnite  element  method.  They  investigated  the  static 
and  dynamic  response  of  FGM  plates  by  varying  the  volume  fraction  of  the  ceramic 
and  metallic  constituents  using  a  simple  power  law  distribution.  They  found  that 
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the  deflection  and  stresses  response  of  the  plates  with  material  properties  between 
those  of  the  ceramic  and  metal  is  not  intermediate  to  the  responses  of  the  ceramic 
and  metal  plates  in  the  case  of  both  thermal  and  mechanical  loading.  However, 
in  the  absence  of  thermal  loading,  they  found  that  the  dynamic  response  of  the 
graded  plates  is  intermediate  to  that  of  the  metal  and  ceramic  plates. 

Dai  et  ah,  [39]  presented  a  meshfree  model  for  the  active  shape  control  as  well 
as  the  dynamic  response  repression  of  FGM  plate  containing  distributed  piezoelec¬ 
tric  sensors  and  actuators  subjected  to  thermo-mechanical  loading.  They  assumed 
that  the  FGM  plate  was  made  of  ceramics  and  metals  whose  volume  fractions  vary 
continuously  in  the  thickness  direction  according  to  the  power  law.  They  used  the 
element-free  Galerkin  method  to  derive  the  shape  functions  using  the  moving  least 
squares  (MLS)  method.  They  studied  the  mechanical  loading  as  well  as  thermal 
gradient.  They  found  that  the  relations  between  the  deflection  and  the  volume 
fraction  exponent  are  quite  different  under  the  two  loadings.  The  theoretical  stud¬ 
ies  presented  above  attempts  to  And  the  displacement  and  stress  response  of  FGMs 
at  different  loading  conditions,  however,  the  studies  are  not  supported  by  exper¬ 
imental  evidence  and  experimental  investigation  on  the  constitutive  behavior  of 
metal/ceramic  FGM  are  not  reported. 

To  date  there  is  not  a  detailed  study  on  the  dynamic  thermo-mechanical  be¬ 
havior  of  FGM,  particularly  on  the  metal-ceramic  based  materials.  Fundamen¬ 
tal  investigation  into  the  thermomechanical  response  and  dynamic  failure  of  the 
graded  materials  is  necessary  before  they  can  be  incorporated  into  the  design  of 
future  space  access  vehicles  that  can  operate  reliably  in  combined  and  extreme 
environments.  The  proposed  research  is  intended  to  All  this  gap  by  a  detailed 
experimental  and  analytical  investigation  of  dynamic  thermo-mechanical  failure  of 
graded  materials. 
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1.3  Present  Study 


The  present  research  is  aimed  at  the  critical  investigation  of  the  dynamic 
fracture  and  failure  properties  of  functionally  graded  material  under  thermo¬ 
mechanical  loading.  The  study  involves  both  analytical  and  experimental  tech¬ 
niques  and  can  be  categorized  in  the  following  list. 

•  Analytical  study  to  develop  a  thermo-mechanical  dynamic  stress  held  around 
the  crack  tip  in  functionally  graded  material  using  high  order  stress  held 
equations. 

•  To  experimentally  investigate  the  dynamic  constitutive  behavior  of  Ti/TiB 
functionally  graded  material  under  thermo-mechanical  loading. 

•  To  experimentally  investigate  the  dynamic  fracture  initiation  toughness  in 
FGM  under  thermo-mechanical  loading. 
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CHAPTER  2 


Dynamic  Thermo-Mechanical  Stress  Fields  for  a  Mixed-Mode 

Propagating  Cracks 

2.1  Abstract 

Thermo-mechanical  stress  held  equations  are  developed  for  a  mixed-mode 
crack  propagating  at  constant  velocity  in  homogeneous  and  isotropic  materials 
using  asymptotic  approach  along  with  displacement  potentials.  Asymptotic  tem¬ 
perature  held  equations  are  hrst  developed  for  steady  state  temperature  conditions 
using  insulating  crack-face  boundary  conditions.  These  temperature  held  equa¬ 
tions  are  later  used  to  derive  the  hrst  three  terms  of  thermo-mechanical  stress  held 
equations  for  a  steady  state  propagating  mixed-mode  crack.  Using  these  thermo¬ 
mechanical  stress  helds,  various  components  of  the  stresses  are  developed  and  the 
ehects  of  temperature  on  these  stress  components  are  discussed. 

2.2  Introduction 

In  many  practical  engineering  problems,  the  structures  are  subjected  to 
thermo- mechanical  loads.  Under  these  thermo- mechanical  loads  the  cracks  in 
these  structures  can  initiate  and  cause  catastrophic  failure.  The  crack-tip  initi¬ 
ation,  rapid  crack  growth,  crack  branching  and  arrest  are  of  signihcant  impor¬ 
tance  to  understand  material’s  failure  under  combined  thermo-mechanical  loads. 
In  the  classical  studies  of  thermoelastic  crack  problems,  Sih  [1]  and  Kassir  and 
Bergman  [2]  investigated  quasi-static  stress  helds  for  a  crack  in  inhnite  medium 
when  it  is  subjected  to  special  thermal  loadings.  Later  Wilson  and  Yu  [3]  employed 
hnite  element  analysis  and  J-integral  approach  to  determine  crack-tip  stress  inten¬ 
sity  factors  for  hnite  specimen  geometries  under  thermal  loads. 

In  continuation  of  the  above  studies,  Lee  and  Sim  [4]  determined  mode-I  ther- 
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mal  shock  stress  intensity  factor  using  Bueckner’  weight  function  method  for  a 
surface  cracked  inhnite  strip  under  sudden  conductive  cooling.  Using  a  general 
hnite  element  model,  Chen  and  Weng  [5],  investigated  a  coupled  transient  ther¬ 
moelastic  problem  for  an  edge-  cracked  plate  without  an  inertia  term.  Katsareas 
et  al.  [6]  determined  shock  stress  intensity  factors  using  boundary-only  element 
method  for  surface  cracked  inhnite  strip  and  hnite  edge  cracked  plate.  Recently, 
Hosseini- Tehran!  et  al.  [7]  investigated  dynamic  mode-I  stress  intensity  factors  for 
an  edge  crack  in  a  hnite  specimen  geometry  using  boundary  element  method  in 
conjunction  with  Laplace  transforms. 

The  above  studies  are  focused  on  the  determination  of  stress  intensity  factor 
values  from  numerical  models  using  integral  transforms.  The  closed-form  solutions 
obtained  using  numerical  and  integral  transform  methods  cannot  be  used  in  ex¬ 
tracting  fracture  parameters  from  experimental  stress  or  deformation  helds.  To 
meet  this  requirement,  asymptotic  expansion  of  stress  held  equations  are  essen¬ 
tial  [8].  Hence  in  this  paper,  thermo-mechanical  stress  helds  are  developed  for  a 
mixed-mode  propagating  crack  of  constant  velocity  in  homogeneous  isotropic  ma¬ 
terials  using  asymptotic  approach  along  with  displacement  potentials.  Asymptotic 
temperature  held  equations  are  hrst  developed  for  a  steady  state  temperature  con¬ 
ditions  using  insulated  crack-tip  boundary  conditions.  These  temperature  helds 
are  later  used  to  derive  hrst  three  terms  of  thermo-mechanical  stress  helds  for 
steady  state  propagating  crack.  Using  these  thermo- mechanical  stress  helds,  vari¬ 
ous  stress  components  are  developed  and  the  ehect  of  temperature  on  these  stresses 
is  discussed. 


13 


2.3  Theoretical  Formulation 

Hooke’s  law  for  a  plane  strain  thermo-mechanical  problem  can  be  written  as 
o'xx  =  ((A  -I-  2^)exx  +  Aeyy  —  (3A  -|-  2fi)aT)  (2.1) 

<^YY  =  ((A  -|-  2/i)eyy  -|-  Xexx  —  (3A  -|-  2jj,)aT)  (2-2) 


txy  —  fJ>^XY 


(2.3) 


where  X  and  Y  are  reference  coordinates,  and  where  i  =  X,Y  and  j  =  X,Y 
are  in-plane  stress  and  strain  components,  A  and  /i  denote  Lame’s  constant  and 
shear  modulus  respectively,  a  is  coefficient  of  thermal  expansion  and  T  represents 
the  change  in  temperature  in  the  inhnite  medium.  The  equations  of  motion  for  a 
plane  problem  in  homogeneous  isotropic  material  are  given  in  Eq.  2.4 


da 


XX 


+ 


dr 


XY 


dX  dY 


d‘^u  dayY  drxY  d‘^v 

dY  ^  dX  ^ ''IFt 


(2.4) 


where  u  and  v  are  functions  of  X,  Y  and  t  and  represent  the  displacements  in  the 
X  and  Y  directions  respectively  with  t  representing  the  time.  For  plane  strain 
deformation,  the  displacements  u  and  v  are  derived  from  dilatational  and  shear 
wave  potentials  0  and  ip.  These  potentials  can  be  expressed  as 


(90  dip  d(p  dip 

ax^W'  '’^W~dx 


(2.5) 


Substituting  Eqs.  2.1  -  2.3  in  to  Eq.2.4  and  utilizing  Eq.2.5  and  after  simplihcation 
the  equations  of  motion  in  terms  of  displacement  potentials  can  be  expressed  as 

(2/i  +  A)  v'  0  -  a(3A  +  2/i)T  =  p— 0  (2.6) 


fi\7^ip 


d^ 


(2.7) 
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Figure  2.1.  Propagating  crack  tip  orientation  with  respect  to  reference  coordinate 
system 

where  =  ^  + 

For  a  propagating  crack  as  shown  in  Fig.  2.1,  the  reference  coordinates  (XandY) 
can  be  transformed  to  crack  tip  coordinates  using  the  relations  x  =  Xc  —  t,  y  =  Y , 
where  c  is  constant  crack  tip  speed.  By  transforming  to  crack-tip  coordinates  and 
rearranging  terms,  the  above  equations  of  motion  can  be  written  as 


2(9V  ^  3(5  +  2^_^ 

Oil  ~  oi  ^  T  —  0 


dx"^  dy"^  5  +  2 


(2.8) 


20^  ,  d'^ii 

=  0 


dx"^  dy"^ 


(2.9) 


“i  =  \/i  -  iim’  “*  =  ^ 

It  is  assumed  that  in  the  above  transformation,  the  helds  0  and  0  do  not  depend 
explicitly  on  time  in  the  moving  coordinate  system  and  their  time  dependence  is 
only  through  the  transformation  x  =  X  —  ct. 


2.3.1  Temperature  Fields  around  the  Crack  Tip 

In  this  analysis  the  thermo-elastic  cooling  and  the  transient  effects  are  ne¬ 
glected.  The  steady  state  heat  conduction  equation  can  be  written  as 


dX  \  dX)  dY 


(2.10) 


Where  k  is  the  coefficient  of  thermal  condnctivity 

Assnming  that  k  is  constant  in  the  small  region  aronnd  the  crack  tip,  and  again 
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transforming  the  above  equation  to  the  crack-tip  moving  coordinate  system  {x  = 
Xc  —  t,y  =  Y),  it  can  be  written  as 


=  0 


(2,11) 


wneie  v  —  dx^  ^  Oy^ 

As  this  stages  asymptotic  analysis  is  performed  to  solve  the  above  Eq.  2.11  .  In 
the  asymptotic  analysis,  first  a  new  set  of  coordinates  is  introduce  as 

X  y 


Vi  = 


V2  = 


e  e 

where  e  is  an  arbitrary  parameter  and  is  assumed  to  be  0  <  e  <  1 
Equation  2.11  can  be  written  in  a  new  scaled  coordinates  {r]i,  772) as 

v"r  =  o 


(2.12) 


(2.13) 


where  v"  =  ||  +  ^ 

For  the  asymptotic  analysis  T  is  represented  as  a  power  series  expansion  in  e  as 


92 


T(x,y)  =  r(i(i6,i(2e)  = 


(2.14) 


m=0 


Substituting  the  inhnite  series  expansion  into  Eq.  2.13  results  in 


E 

m=0 


V  <9^1  <9^2 


(2.15) 


For  Eq.  2.15  to  be  valid,  the  partial  differential  equations  corresponding  to  each 
power  of  e  e,  e^/^...)  should  vanish  independently.  This  leads  to  a  set  of  partial 
differential  equations.  Assuming  the  crack  surfaces  are  insulated  (i.e  dT/d6=0  at 
6  =  7i),  the  solutions  for  the  hrst  three  temperature  terms  can  be  written  as 

T  =  qop^^^sin  06*^  +qipcos{9)  +  q2p^^‘^  sin  06*^  (2.16) 

-1  /Q 

where  p={pl  +  pi)  go,  qi  and  g2  are  real  constants 

By  transforming  T  into  crack-tip  coordinates  x  and  y,  the  temperature  helds  near 
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the  crack-tip  can  be  given  as 


T  =  qor^^'^  sin  ^os  (6*)  -|-  q2r^^'^  sin 

1  /  ‘y 

where  r={x‘^  -|-  y"^)  '  and  tan(0)=  - 


(2.17) 


2.3.2  Asymptotic  Expansion  of  Crack  Tip  Fields 

The  Eqs.  2.8  and  2.9  are  now  written  in  these  scaled  coordinates  as  below 


2(920  ^  q2^  (3(5  +  2) 

0^1  o  9  T  9  —  Cl “77 - TvTC  9  —  0 

'  dyf  dyl  (5  +  2) 


(2.18) 


2^  <920  _ 


dyf  '  dyl  ^ 

At  this  stage  it  is  assumed  that  0,  0  and  T  are  represented  as  a  power  series 


expansion  in 


OO 


0(a;,  y)  =  0(er/i,  er/2)  =  X]  e  ^  0^(7^i,  r/2) 

m=0 

(2.20) 

OO 

^  (m+3) 

0(a;,  y)  =  0(er/i,  er/2)  =  2^  e  2  0m(r/i,  V2) 

m=0 

(2.21) 

00 

^ ^  (m  +  l) 

T{x,  y)  =  T{eyi,  €^2)  =  "  70(^i,  V2) 

m=0 

(2.22) 

As  given  by  Jin  and  Noda  [9]  the  solution  for  the  temperature  held  is  obtained  by 
considering  the  singularity  of  heat  hux  near  the  crack  tip.  This  requires  in  that  the 
derivative  of  temperature  held  is  singular  resulting  in  the  potential  for  temperature 
to  be  one  power  less  than  the  displacement  potentials. 

Now  substituting  Eqs.  2.20  -  2.22  into  Eqs.  2.18  and  2.19  gives  the  following 
equations 


OO 

g(m+3)/2 


m=0 


Oi] 


drjl  dri2 


a- 


OO 


(3(5  +  2) 

(<5  +  2) 

^  '  m=0 


,{m+5)/2rp  _  g 


(2.23) 
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,(m+3)/2 


m=0 


For  Eqs.  2.23  and  2.24  to  be  valid,  the  partial  differential  equations  corresponding 
to  each  power  of  e  e,  e^/^...)  should  vanish  independently.  This  leads  to  a 

set  of  partial  differential  equations.  The  solutions  for  m  =  0  and  m  =  1  are  well 
known  and  available  in  literature  [10]  and  can  be  written  as 


cos 


[m  +  3) 


+  CmPl 


(m+3)/2 


Sin 


[m  +  3) 


Oi)  (2.25) 


'lpm{Ps,0s,t)  =  sin  cos 

(2.26) 

where  pi=ipj  +  tan(6'«)=  Ps={r]f  +  tan(6'^)=  ^  and  A^, 

Bm,  Cm,  and  Dm  are  real  constants 

The  solution  for  equation  corresponding  to  higher  power  of  (m  =  2)  consists  of  two 
parts:  classical  solution  and  particular  solution.  The  classical  solution  is  similar 
to  the  general  solution  given  in  Eqs.  2.25  and  2.26.  The  particular  solution  can  be 
obtained  using  recursive  approach  [10].  The  solutions  for  m  =  2  (02  and  02)  are 
given  below. 


=  A2r\^‘^  cos 


+  C2t\^‘^  sin 


5. 

2^' 


4  (35  +  2)  a  5/2  .  5^ 

+  15^  +  2) 


(2.27) 


0.)  (2.28) 

By  substituting  the  expressions  for  0o,  0i,  02,  0o,  0i,  and  02  into  the  Eqs.  2.20- 
2.22  and  representing  the  Eqs. 2. 20  and  2.21  in  terms  of  crack-tip  coordinates,  the 


02  =  B^r^C  sin  cos  ^ 
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displacement  potentials  can  be  represented  as 


(j)  =  cos  sin  cos  (26*;)  +  Cirf  sin  (26*;) 


+  cos  (  i^6*;  )  +  Carj’^^  sin  (  ^6*;  )  + 


5 


5/2 


5 


15  (5  +  2)  (a2_i)^or  sin^y 

(2.29) 


ij  =  5orf  ^  sin  cos  sin  (26*J  +  +>ir^  cos  (26*J 

+  sin  cos 


(2.30) 


where  ri={x‘^  +  tan(6*;)=  r^=(a;2  +  a2|/2)V2^  tan(6*J=  ^ 

2.3.3  Thermo-mechanical  Stress  Fields 

By  substituting  Eq.  2.29  and  Eq.  2.30  into  the  Eq.  2.5,  the  displacement 
held  can  be  obtained.  Further  differentiating  the  displacement  field,  the  in-plane 
strain  field  can  be  obtained.  The  strain  expressions  along  with  Eq.  2.17  can  be 
substituted  into  Eqs.  2.1-  2.3  to  obtain  in-plane  thermo-mechanical  stress  fields 
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a. 


yy 


3  -1/2 

-r,  cos 

4  ' 


(5  -  of  {5  +  2))^  Ao  -  {^r-i  sin  (5  -  {6  +  2))^  Cq 

^  Do  -  4Ai  {S  -  {6  +  2)) 


-  (^r;>/=cos  (|)a, 

-  4Dia,  +  cos  (5  -  af  (5  +  2))^  A2 

(f)  (^-a|(^  +  2)))c,-(^ 


15  1/2  . 
— r,  sm 
4  ' 


e. 


rV"^  cos  (  y  )  )  B2 


9  39 

q;(3(5  +  2)  (  sin  -  +  qir  cos{6)  +  q2r^^‘^  sin  — 


15  _i/2  .  fOs\  \  ^  ^/3(5  +  2  a  1/2  . 

y-' '  =■"( 2  “•  - 2  T-Tsry- '  y  1  "0 


(2.32) 


r. 


/i 


3  -1/2  .  f  A  j_  -1/2  ( ^^  \  \  r 

-aiTi  sm  -  Ao  +  -aiTi  cos  -  Do 


'3 

+  I  4D 


2 

'sill  (  Y  )(1  +  «D  )  Bo 


/^^-i/2 
4  ^ 


cos  (  Y  )  (1  +  )  Do 


+  4Dia/  -  2Di  (1  -  a^)  -  (^air]D  gin  ^  A2 

+  (y“'D~=(|))c2-( 


^r/^sin  fj  )(1  +  q;J)  | 


T'l'  (I)  1  +  “*))  +  2  ( 5 


Qo 

(2.33) 


By  considering  the  singular  term  in  the  above  equation  and  using  the  dehnition 
of  the  dynamic  stress  intensity  factor  K/d  and  Kjjd  for  opening  and  shear  modes 
and  considering  the  crack  face  boundary  conditions  [11]  we  get 

4(1  +  a^)  Kj 


An  — 


^ID 


Aaias  —  (1  +  tts)^ 


p  _ 


(2.34) 


Dn  = 


4q;^ 


K 


I  ID  n  _  1  ^ 

-^0  —  — - v-o 


(2.35) 


°  S^iaiDs  —  {1  +  a'ly)  "  2q;s 

where  p  is  the  shear  modulus  of  the  material  at  the  propagating  crack-tip,  Kjo 
is  mode-I  dynamic  stress  intensity  factor  and  Kud  is  mode-II  dynamic  stress 
intensity  factor. 
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2.4  Results  and  Discussion 

2.4.1  Variation  of  Stress  Components  Near  the  Crack  Tip 

Using  the  developed  stress  fields,  various  components  of  the  stress  tensor  are 
evaluated  and  the  effect  of  temperature  on  these  stress  components  is  discussed. 
The  coefficient  of  thermal  expansion  (a)  of  titanium  (which  is  of  primary  interest 
in  our  experimental  research),  8.9  x  is  used.  In  all  the  plots  only  the  first 

terms  are  considered. 

The  superposition  of  the  temperature  stress  field  with  the  mechanical  stress 
field  results  in  imposing  normal  stress  on  traction  free  crack  faces  as  shown  in 
Fig.  2.2a.  These  stresses  on  the  crack  face  are  generated  due  to  the  thermal 
stresses  produced  by  the  temperature  field.  The  presence  of  these  normal  stresses 
on  the  crack  face  violates  one  of  the  boundary  conditions  and  these  stresses  need 
to  be  removed.  The  removal  of  these  stresses  from  the  crack  face  is  accomplished 
by  superimposing  an  equal  and  opposite  stress  field  on  the  crack  face.  In  particular 
each  point  on  the  crack  face  is  subjected  to  a  line  load  of  specific  magnitude  such 
that  normal  stresses  from  the  crack  face  are  removed.  The  solution  given  by  Jiang 
et  ah  [12]  for  a  single  line  load  is  extended  to  derive  the  stress  field  around  the 
crack  tip  for  multiple  line  loads  applied  on  the  crack  face  line.  Fig.  2.2b  represents 
normal  stress  {(^yy)  near  the  crack  tip  after  the  supper  position  has  been  added.  It 
must  be  mentioned  that  the  crack  face  loads  are  relatively  small  in  this  example 
even  for  large  T  and  have  little  influence  on  the  crack  tip  stress  fields. 

The  thermo-mechanical  stress  fields  these  obtained  by  superposition  of 
the  thermal  and  mechanical  stress  fields  and  suitably  modified  to  ensure  proper 
boundary  conditions  were  used  to  generate  contours  of  in  plane  stresses,  maximum 
principal  stress,  maximum  shear  stress  and  maximum  tangential  stress  around  the 
crack  tip. 


21 


0 


X  10*  0 


X  (m) 


X  10^ 


Figure  2.2.  Contour  of  ayy  in  MPa  for  thermo-mechanical  stress  held  around  the 
crack-tip  (go=250,  c/cg  =  0.5  and  Kjd  =  lMPa-ml/2)(a)  Before  the  opposite 
traction  held  is  applied  (b)  After  the  opposite  traction  held  is  applied. 

Fig.  2.3  shows  the  various  normalized  in-plane  stress  components  (a)  in-plane 
o'xx,  (b)  in-plane  ayy  ,(c)  in-plane  Txy  ,  and  Fig.  2.4  shows  the  various  normalized 
principal  stress  (a)  hoop  stress  agg,  (b)  maximum  principal  stress  ai  and  (c)  max¬ 
imum  shear  stress  Tmax  as  a  function  of  angular  position  9  around  the  crack  tip 
for  a  hxed  value  of  c/cg  and  diherent  temperature  helds.  The  stresses  are  nor¬ 
malized  with  far-held  stress  determined  for  =1  MPa-m^/^,  Kjju  =  0.27^/^ 
-|-  Kjjjj  and  r  =  0.002m.  The  stress  held  is  assumed  to  be  elastic 
for  the  temperature  held  associated  with  qq  =  100  and  go  =  250. 

In  Fig.  2.3  and  Fig.  2.4,  go  =  0  means  no  temperature  held  around  the  crack- 
tip  and  for  this  value  Fig.  2.4  collapse  to  the  isothermal  solution  given  by  Freund  [8] 
if  mode-I  and  mode-II  helds  are  separated  and  plotted  as  shown  in  Figs.  A.l  -  A. 6 
in  Appendix  A. 

The  in-plane  stress  axx  shown  in  Fig.  2.3(a)  has  a  maximum  value  ahead  of  the 
crack  tip,  i.e.  at  6=  0°  and  a  local  maximum  about  ±  135°  for  qo=  0.  These  local 
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e  (degree) 


e  (degree) 


Figure  2.3.  Normalized  mixed  mode  in-plane  stress  fields  as  a  function  of  6  around 
the  crack-tip  for  different  temperature  fields  {c/cg  =  0.5,  r  =  0.002m,  Kj£,  =  IMPa 
Kiid  =0.2  Kjd) 
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Figure  2.4.  Normalized  mixed  mode  principal  stress  fields  as  a  function  of  9  around 
the  crack-tip  for  different  temperature  fields  {c/cg  =  0.5,  r  =  0.002m,  Kjo  =  IMPa 
Kiid  =0.2  Kjd) 
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maximums  increase  in  magnitude  with  an  increase  in  q^.  It  is  also  observed  that 
the  effect  of  temperature  on  a^x  is  predominantly  in  the  region  near  to  the  crack 
face.  The  in-plane  stress  Cyy  shown  in  Fig.  2.3(b)  has  a  maximum  value  around 
-70°  for  (Jo  =  0  and  this  value  decreases  with  an  increase  in  applied  temperature. 
It  also  has  a  local  maximum  at  about  60°  and  this  value  increases  with  an  increase 
in  applied  temperature.  The  in-plane  shear  stress  Txy  shown  in  Fig.  2.3(c)  has  a 
maximum  value  around  -110°  and  a  local  maximum  around  20°.  These  maximum 
values  decreases  as  go  increases.  There  is  no  shift  in  the  position  of  the  maximum 
value  of  in-plane  stresses  as  the  temperature  held  around  the  crack-tip  changes. 

For  qo=  0,  the  hoop  stress  as  shown  in  Fig.  2.4(a)  has  a  single  maximum 
value  along  6  =  -30°.  As  the  temperature  held  increases,  this  stress  develops  a 
maximum  and  local  maximum  values.  For  example  for  qo=  250,  the  hoop  stress 
develops  a  maximum  value  at  an  angle  of  30°  and  a  local  maximum  value  at  angle 
of  -120°  to  the  direction  of  crack  growth.  The  angular  variation  of  the  principal 
stresses  around  the  crack  tip  is  shown  in  Fig.  2.4(b).  For  no  temperature  held,  this 
stress  has  a  maximum  value  at  about  -90°  and  local  maximum  value  at  about  45° 
from  the  direction  of  crack  growth.  As  the  temperature  held  increases,  however, 
the  maximum  value  shifts  towards  the  direction  away  from  the  crack  growth.  In 
Fig.  2.4(c),  the  maximum  shear  stress  as  a  function  of  temperature  held  around  the 
crack  tip  is  shown.  In  the  case  of  no  temperature  held  (i.e.,  qo=  0),  the  maximum 
shear  stress  has  a  maximum  value  at  about  -115°  and  a  local  maximum  value  at 
about  85°  from  the  direction  of  crack  growth.  It  can  be  seen  from  the  hgure  that 
as  the  temperature  increases,  maximum  values  occur  at  four  diherent  locations. 

2.4.2  Direction  of  Crack  Growth 

The  theoretical  prediction  of  crack  extension  angle  is  investigated  by  using 
the  two  well-known  fracture  criterias:  minimum  strain  energy  density  (S-criterion) 


25 


and  maximum  circumferential  stress  {agg  -criterion). 


Minimum  Strain-energy  Density  (MSED)  Criterion 

According  to  this  criterion  [13],  the  crack  initiates  when  the  strain  energy 
density  achieves  a  critical  value  and  propagates  in  the  direction  of  minimum  strain- 
energy  density  value.  The  strain  energy  density  dW/dV  near  the  crack  tip  for  an 
FGM  is  given  as 

^  ^  ((1  -  ^)  +  ^ly)  -  +  <)  (2.36) 

Fracture  takes  place  in  the  direction  of  minimum  S,  and  the  condition  can  be 
obtained  by  using  Eq.  2.37. 

nc  02  c 

-  =  0,  ^>0  at  S  =  S^  (2.37) 

where  Sc  is  the  critical  strain  energy  density 


Maximum  Circumferential-stress  (MCS)  Criterion 

The  maximum  circumferential  stress  criterion  [14]  states  that,  crack  growth 
will  occur  in  the  direction  of  the  maximum  circumferential  stress  and  will  take 
place  when  the  maximum  circumferential  stress  reaches  a  critical  value,  and  it  can 
be  given  as  Eq.  2.38. 


da, 


ee 


89 


=  0, 


d^agg 

89^ 


<  0  at  agg  —  (agg)c 


(2.38) 


where  {agg)c  is  the  critical  circumferential  stress 

Based  on  the  above  two  criteria’s,  the  effects  of  velocity  and  temperature  on  the 
crack  extension  angle  {9)  are  investigated. 

The  crack  extension  angles  as  a  function  of  crack  tip  velocities  as  predicted 
by  the  above  two  criterias  are  shown  in  Fig.  2.5.  For  pure  mode-I  loading 


26 


{Kud/ Kjd  =  0),  the  crack  extends  along  (0)=O  until  the  crack  tip  velocity  reaches 
a  critical  value  at  which  instability  occurs  [15].  When  the  crack  tip  velocity  reaches 
the  critical  value,  the  crack  deviates  and  extends  to  a  different  angle.  For  example 
at  a  crack  tip  velocity  of  c/c^  =  0.7,  the  MSED  criterion  predicts  a  crack  extension 
angle  of  about  —55"  and  the  MCS  criterion  predicts  about  —38°.  As  the  value 
of  Kuo/ Kid  increases  from  0  to  1  and  later  from  1  to  cx)  the  crack  extension 
angle  increases  monotonically.  Broek  [16],  in  his  book  gives  the  crack  extension 
angles  for  mixed  mode  quasi  static  loading  and  these  results  match  well  with  the 
predictions  from  the  current  study. 

The  effect  of  temperature  field  on  the  crack  extension  angle  for  a  crack  tip 
velocity  of  O.bc^  is  shown  in  Fig.  2.6.  The  hgure  shows  the  predicted  of  crack 
extension  angle  by  both  the  minimum  strain  energy  density  and  the  maximum 
circumferential  stress  criterions.  These  two  criterions  predict  same  crack  extension 
angles  at  lower  temperatures.  Both  the  criterion  show  that,  the  crack  extension 
direction  at  room  temperature  is  along  9  =  —30°  and  the  value  decreases  slowly 
with  increase  in  applied  temperature  held.  For  FGM  with  a  >  0,  the  crack  exten¬ 
sion  angle  is  along  9  =  —15°  at  room  temperature  and  again  the  value  decreases 
with  increase  in  applied  temperature  held. 

The  ehects  of  temperature  on  the  crack  extension  angle  as  a  function  of  crack 
tip  velocity  are  plotted  in  Fig.  2. 7. Fig.  2.7(a)  represents  crack  extension  angle  as 
a  function  of  crack  tip  velocity  for  mixed-mode  mechanical  loading  at  room  and 
elevated  temperatures  predicted  by  minimum  strain  energy  density  criterion.  For  a 
crack  tip  velocity  of  c/cg  =  0.3,  the  minimum  strain  energy  density  criterion  predict 
decrease  in  crack  extension  angle  for  increasing  temperature  helds.  Fig.  2.7(b) 
represents  crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed-mode 
mechanical  loading  at  room  and  elevated  temperatures  predicted  by  maximum 
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(b)  Maximum  circumferential  stress  criterion 


Figure  2.5.  Crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed 
mode  thermo-mechanical  loading  (r=0.002m). 


Figure  2.6.  Effect  of  temperature  on  the  crack  extension  angle  for  mixed  mode 
loading  {Kjjd/ Kjj:,  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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circumferential  stress  criterion.  Furthermore,  for  cjcg  >  0.3,  both  criterions  predict 
decrease  in  crack  extension  angle  for  increasing  temperature  held.  The  difference 
in  the  predicted  crack  extension  angles  decreases  as  the  velocity  increases  and 
becomes  insignihcant  at  velocities  c/c^  >  0.55. 

2.5  Summary 

The  stress-helds  near  the  crack  tip  for  mixed-mode  thermo-mechanical  load¬ 
ing  are  developed  using  displacement  potentials  in  conjunction  with  an  asymptotic 
approach.  These  crack-tip  held  equations  can  be  successfully  used  in  extracting 
fracture  parameters  from  experimental  stress  and  deformation  helds  when  ma¬ 
terials  are  subjected  to  thermo-mechanical  loading.  The  following  are  the  key 
observations  of  the  analytical  stress  held  contours  developed  using  the  stress  held 
equations. 

•  The  temperature  held  changes  the  magnitude  of  the  in-plane  stress  compo¬ 
nents  o'yy  and  r^y)  but  the  qualitative  prohle  around  the  crack  tip  is 
not  inhuenced. 

•  The  temperature  held  inhuences  both  the  magnitude  and  the  prohle  of  max¬ 
imum  shear  stress  r^ax,  hoop  stress  agy  and  the  principal  stress  ui  around 
the  crack-tip. 

•  The  crack  growth  angle  decreases  with  increase  in  temperature  for  a  given 
crack-tip  velocity. 

•  The  crack  growth  angle  increases  with  increase  in  crack-tip  velocity  for  a 
given  temperature  conditions  around  the  crack-tip. 

•  The  ehect  of  temperature  on  the  crack  growth  angle  is  insignihcant  at  higher 
crack-tip  velocities  c/cg  >  0.55. 
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(a)  Minimum  strain  energy  density  criterion 


(b)  Maximum  circumferential  stress  criterion 

Figure  2.7.  Crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed 
mode  thermo-mechanical  loading  {Kjjjj/Kjo  =  0.2,  r  =  0.002m). 
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CHAPTER  3 


Effect  of  Dynamic  Thermo-Mechanical  Stress  Field  on  Crack 
Propagation  Direction  in  Functionally  Graded  Materials 

3.1  Abstract 

Asymptotic  analysis  in  conjunction  with  displacement  potentials  has  been 
used  to  develop  thermo-mechanical  stress  helds  for  a  mixed  mode  propagating 
crack  in  a  functionally  graded  material  (FGM).  The  shear  modulus,  mass  density, 
thermal  conductivity  and  coefficient  of  thermal  expansion  of  the  FGM  are  assumed 
to  vary  exponentially  along  the  gradation  direction.  First,  asymptotic  tempera¬ 
ture  helds  are  derived  for  an  exponential  variation  of  thermal  conductivity  and 
later  these  temperature  helds  are  used  in  deriving  stress  helds.  Using  asymptotic 
thermo-mechanical  stress  helds  the  variation  of  maximum  shear  stress,  circum¬ 
ferential  stress  and  strain  energy  density  as  a  function  of  temperature  around  the 
crack-tip  are  generated.  Finally,  utilizing  the  minimum  strain  energy  density  crite¬ 
rion  and  the  maximum  circumferential  stress  criterion,  the  crack  growth  direction 
for  various  crack-tip  speeds,  non-homogeneity  coefficients  and  temperature  helds 
are  determined. 

3.2  Introduction 

Functionally  graded  materials  (FGMs)  are  essentially  non-homogeneous  com¬ 
posites  which  have  characteristics  of  spatially  varying  microstructure  and  mechan¬ 
ical/thermal  properties  to  meet  a  predetermined  functional  performance  [1,  2]. 
Although  their  performance  in  real-life  engineering  applications  are  still  under 
investigation,  FGMs  have  shown  promising  results  when  they  are  subjected  to 
thermo- mechanical  loading  [2] .  Hasselman  and  Youngblood  [3]  were  among  the  hrst 
to  study  thermal  stresses  in  nonhomogeneous  structures  associated  with  thermo- 
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mechanical  loading.  By  introducing  thermal  conductivity  gradient,  they  realized 
signihcant  reductions  in  the  magnitude  of  the  tensile  thermal  stress  in  ceramic 
cylinders.  In  other  studies,  thermal  residual  stresses  are  relaxed  in  metal-ceramic 
layered  materials  by  inserting  a  functionally  graded  interface  layer  between  the 
metal  and  ceramic  [4-6].  In  their  studies,  Kudora  et  ah,  and  Takashashi  et  ah, 
reported  that  when  subjected  to  thermal  shocks,  FGM  coatings  suffer  signihcantly 
less  damage  than  conventional  ceramic  coatings  [7,  8]. 

In  continuation  of  above  studies,  several  studies  on  the  quasi-static  fracture 
of  FGMs  under  thermo-mechanical  loading  have  been  reported.  Assuming  expo¬ 
nential  variation  of  material  properties,  Jin  and  Noda  [9]  investigated  the  steady 
thermal  stress  intensity  factors  in  functionally  gradient  semi-inhnite  space  with  an 
edge  crack  subjected  to  thermal  load.  Later,  Erdogan  and  Wu  [10]  also  determined 
the  steady  thermal  stress  intensity  factor  of  a  FGM  layer  with  a  surface  crack  per¬ 
pendicular  to  the  boundaries.  By  further  assuming  the  exponential  variation  of 
thermal  and  mechanical  properties  of  the  materials,  Jin  and  Batra  [11]  investigated 
stress  intensity  relaxation  problem  at  the  tip  of  an  edge  crack  in  a  FGM  subjected 
to  a  thermal  shock.  Using  both  experimental  and  numerical  techniques,  Kokini 
and  Ghoules  [12]  and  Kokini  and  Gase  [13]  studied  surface  and  interface  crack¬ 
ing  in  FGM  coatings  subjected  to  thermal  shocks.  By  employing  a  hnite  element 
method  (FEM),  Noda  [14]  analyzed  an  edge  crack  problem  in  a  zirconia/titanium 
FGM  plate  subjected  to  cyclic  thermal  loads.  Using  hnite  element  method  and 
boundary  element  method,  Jin  and  Paulino  [15]  studied  transient  thermal  stresses 
in  an  FGM  with  an  edge  crack  and  having  constant  Young’s  modulus  and  Pois¬ 
son’  ratio  but  varying  thermal  properties  along  the  thickness  direction.  Walters  et 
al.  [16]  developed  general  domain  integral  methods  to  obtain  stress  intensity  factors 
for  surface  cracks  in  FGMs  under  mode-I  thermo-mechanical  loading  conditions. 
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The  above  studies  provide  closed  form  solutions  for  stress  intensity  factors  under 
thermo-mechanical  loading,  however  for  extracting  fracture  parameters  from  exper¬ 
imental  studies,  asymptotic  expansion  of  thermo-mechanical  stress  helds  around 
the  crack  tip  are  essential.  In  this  direction,  recently,  Jain  et  ah  [17]  developed 
quasi-static  stress  and  displacement  helds  for  a  crack  in  an  inhnite  FGM  medium 
under  thermo-mechanical  loading. 

In  this  paper,  the  stress  helds  for  a  propagating  crack  at  uniform  speed 
along  the  direction  of  mechanical  and  thermal  property  variation  in  a  FGM  under 
thermo-mechanical  loading  conditions  are  developed.  The  elasto-dynamic  problem 
for  FGM  is  formulated  in  terms  of  displacement  potentials  and  the  solutions  are 
obtained  through  an  asymptotic  analysis.  In  analyzing  this  problem,  we  transform 
the  general  partial  diherential  equation  in  the  dynamic  equilibrium  into  Laplace’s 
equation  whose  solution  involves  harmonic  functions.  Using  the  developed  equa¬ 
tions,  angular  variation  of  the  maximum  shear  stress,  circumferential  stress  and 
minimum  strain  energy  density  are  plotted  as  a  function  of  temperature  around 
the  crack-tip.  Using  both  minimum  strain  energy  density  criterion  and  maximum 
circumferential  stress  criterion,  the  crack  growth  directions  for  various  crack-tip 
speeds,  non-homogeneous  coefficients  and  diherent  temperatures  helds  are  also 
determined. 

3.3  Theoretical  Formulation 

At  a  continuum  level,  the  properties  at  any  given  point  in  an  FGM  can  be 
assumed  to  be  same  in  all  directions;  hence  FGMs  can  be  treated  as  isotropic 
non-homogeneous  solids.  Spatial  variation  of  elastic  properties,  mass  density  and 
thermal  properties  make  analytical  solutions  to  the  elastodynamic  eqnations  ex¬ 
tremely  difficult.  Hence,  an  asymptotic  analysis  similar  to  that  employed  by  Fre- 
nnd  [18]  is  used  to  expand  the  stress  held  around  a  propagating  crack  nnder  thermo- 
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mechanical  loading  conditions.  The  complete  derivation  is  provided  in  Appendix 

C. 

Shear  modulus  (/i),  Lame’s  constant  (A),  density  (p)  ,  thermal  expansion  (a)  and 
heat  conductivity  {k)  of  the  FGM  are  assumed  to  vary  in  an  exponential  manner 
as  given  by  Eq.  3.1,  whereas,  Poisson’s  ratio  (z/)  is  assumed  to  be  a  constant. 

fi  =  X  =  p  =  poe^‘'^\  a  =  aoe^^^\  k  =  koe^^^^  (3.1) 

The  equations  of  motion  for  a  plane  problem  are  given  by  Eq.  3.2 

daxx  drxY  _  d^u  dayy  drxy  _  d^v  ^ 

dX  ^  dY  ~^dH'  dY  ^  dX  ~^dH  ^  ■  ’ 

The  relationship  between  stresses  and  strains  for  a  plane  strain  thermo-mechanical 

problem  can  be  written  as 

(^xx  =  e*'^^^((Ao  -l-  2/io)exx  +  Xeyy  —  (3Ao  -l-  2pQ)aQe^^^^T)  (3.3) 


ayy  —  e^‘’^^((Ao  -l-  2po)eyy  +  Xexx  —  (3Ao  -l-  2po)aoe^^^^T)  (3.4) 


Txy  =  (3.5) 

where  X  and  Y  are  reference  coordinates, and  where  i  =  X,Y  and  j  =  X,Y 
are  in-plane  stress  and  strain  components,  A  and  p  denote  Lame’s  constant  and 
shear  modulus  respectively  and  subscript  ”o”  means  at  X  =  0  as  shown  in  Eig.  C.l. 
a  is  coefficient  of  thermal  expansion  and  T  represents  the  change  in  temperature 
in  the  inhnite  medium  and  (  and  (3  are  nonhomogeneity  constants  that  have  the 
dimension  {length)~^ .  Eor  plane  strain  deformation,  the  displacements  u  and  v 
are  derived  from  dilatational  and  shear  wave  potentials  0  and  These  potentials 
can  be  expressed  as 


(90  (90  (90  (90 


(3.6) 
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Figure  3.1.  Propagating  crack  tip  orientation  with  respect  to  reference  coordinate 
configuration. 


Substituting  for  the  stresses  and  density  from  Eqs.  3.3  -  3.5  and  Eq.  3.1  respec¬ 
tively  into  Eq.  3.2  and  after  simplification,  the  equations  of  motion  become 

/  Q  d  \ 

(2po  +  Ao)  VV  -  C  (  (2/Uo  +  Ao)^  0  +  ^  ) 

—  q;c(3Ao  -I-  2po)  +  (2/3  -f  4>  (3-7) 

f  d  d  \  dT 

Po  -  C  (  (Po^  ^  +  Ao^  0  j  -  ac(3Ao  +  2po)C^  =  Po^  ^ 

(3.8) 


ttc  is  the  coefficient  of  thermal  expansion  in  the  vicinity  of  the  instantaneous  crack 
tip  and  is  assnmed  to  be  constant. 

Eqnations  3.7  and  3.8  can  be  written  as 

(5  +  2)  v"'(--c((5  +  2)A^  +  ^) 

—  q;c(3Ao -I- 2/io)  ^^(/5  +  C)/^^  -I- (2/3 -f  C) A  (^•^) 


v^^-C 


d'lp 


ctc(3Ao  -l-  2po)CA  ^ 


dT 

W 


Po  ^ 
/io  dH 


(3.10) 
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where 


v^’ 


ho 


For  a  propagating  crack  shown  in  C.l,  the  transformed  crack  tip  coordinates  can 
be  written  as  x  =  X  —  ct  where  c  is  constant  crack  tip  speed.  In  the  moving 
coordinate  systems  the  above  Eqs.  3.9  and  3.10  can  be  written  as 

2^  ^  ^  C  djj 

dx‘^  dy"^  dx  5  +  2  dy 

-  (((/5  +  <m-'T  +  (213  +  OA-'H)  +  r)  =  0  (3.11) 


— 1“  - ac{S6  +  2)(A  —  =  0 

ox^  oy^  ox  ay  oy 


(3.12) 


where 


ai  — 


1-  (  - 
Cl 


a.^  =  \  1  -  \  — 


*  V  he  ’  ^  dx"^  dy"^ 


Cl  =  C5a/2(1  —  z/)/(1  —  2z/)  for  plane  strain  and  q  =  Cs a/2/ (1  —  u)  for  plane  stress. 
Q  and  Cs  are  the  elastic  dilatational  wave  speed  and  the  elastic  shear  wave  speed 
of  the  material  at  the  crack  tip. 

It  is  assumed  that  in  the  above  transformation,  the  helds  0  and  ip  do  not 
depend  explicitly  on  time  in  the  moving  coordinate  reference  and  their  time  de¬ 
pendence  is  only  through  the  transformation  x  =  X  —  ct.  To  solve  the  above 
equations  of  motion  for  displacement  potentials  (p  and  ip  using  asymptotic  ap¬ 
proach,  Erst  the  derivation  of  temperature  held  around  the  crack-tip  is  essential. 
Hence  its  derivation  is  discussed  in  the  following  section. 


3.3.1  Temperature  Fields  around  the  Crack  Tip 

In  this  analysis  it  is  assumed  that  the  temperature  held  around  the  crack 
tip  changes  asymptotically.  Also,  the  transient  effects  are  neglected.  The  heat 
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conductivity  is  assumed  to  vary  exponentially  as  given  by  Eq.  3.1.  The  steady 
state  heat  conduction  equation  can  be  written  as 


dX  \^dX ) 


(3.13) 


Where  k  is  the  coefficient  of  thermal  conductivity 

Assuming  that  k  is  constant  in  the  region  considered,  the  above  equation  can  be 
written  as 

V^T  +  /3— =0  (3.14) 

where 

Transforming  the  above  equation  to  the  crack-tip  moving  coordinate  system  {x  = 
Xc  —  t,y  =  Y),  Eq.  3.14  can  be  written  as 

V^T  +  /3— =  0  (3.15) 

where  v'  =  ^ 

As  this  stages  asymptotic  analysis  is  performed  to  solve  the  above  Eq.  3.15.  In 
this  process,  hrst  a  new  set  of  coordinates  is  introduce  as 


^  y  /'Q 1 

hi  =  -  ^2  =  -  (3.16) 

where  e  is  an  arbitrary  parameter  and  is  assumed  to  be  0  <  e  <  1 
Equation  3.15  can  be  written  in  a  new  scaled  coordinates  (?7i,?72)as 

8T 

V^T  +  e/3— =  0  (3.17) 

orji 

wnere  v  -  ^ 

For  the  asymptotic  analysis  T  is  represented  as  a  power  series  expansion  in  e  as 

oo 

T{x,  y)  =  T{riie,  r]2e)  =  e^^T™(77i,  172)  (3.18) 

m=0 
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Substituting  Eq.3.18  into  Eq.3.17  gives  the  following  equation. 


±1)  ax 

VV  dvf  dr^l 


,  (""+3)  ^  dTm 

orji 


(3.19) 


For  Eq.3.19  to  be  valid,  the  partial  differential  equations  corresponding  to  each 
power  of  e  e,  e^/^...)  should  vanish  independently.  This  leads  to  the  set  of 

partial  differential  equations. 

For  m  =  0  and  m  =  1 


For  m  =  2 


d^Tm  _ 

drjl  drjl 


,  d^T2  ,  ^dTo 


+  /9X  =  0 


(3.20) 


(3.21) 


d7]l  drjl  drjl 

Equation  3.20  (i.e.  for  m  =  0  and  m  =  1)  is  ordinary  Laplace’s  equation  in  the 
domain  p  =  rji+irj2  and  the  solution  is  same  as  for  homogenous  material  [19,  20]. By 
assuming  insulated  crack  surface  boundary  condition,  (i.e.  ^  =  0  at  0  =  Ttt  ), 
the  solution  for  Eq.3.20  can  now  be  written  as 


For  m  =  0 


For  m  =  1 


To  =  goP^'^^sin 


Ti  =  cos  (9) 


(3.22) 


(3.23) 


where  p  =  {pf  + 


The  solution  T2,  for  Eq.3.21  corresponding  to  higher  powers  of  e  (m  =  2)  has 
two  parts:  homogeneous  and  particular  solution.  The  particular  solution  can  be 
obtained  using  recursive  approach  [19,  20],  and  the  complete  solution  for  Eq.3.21 


is  given  below. 


T2  =  g2P^/^sin  00^  +  qopp^/^in  (3.24) 
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Transforming  back  to  crack  tip  coordinates  x  and  y,  the  temperature  field  near  the 
crack  tip  is  given  as 


T  =  sin  ^2  +qi^  cos  {6)+q2r'^^‘^  sin  ^^6*^  +  qoP'T^^'^  sin  (3.25) 

where  r={x‘^  +  and  Q=  tan“^ 


3.3.2  Asymptotic  Expansion  of  Crack  Tip  Stress  Fields 

Similar  to  derivation  of  temperature  held,  asymptotic  approach  is  again  used 
in  deriving  solutions  for  displacement  potentials  of  equations  of  motion  Eqs.  3.11 
and  3.12.  Scaled  coordinates  as  discussed  in  Eq.3.16  are  now  applied  to  Eqs.  3.11 
and  3.12  as  given  below 


2^  /  (90  C 

dr]l  drjl  ^  \  drji  S  +  2  dr]2  ) 

acf^0|^(£"r  +  €q2^(  +  OA-‘|l+€-‘(/3  +  C)/3A-‘r)  =0  (3.26) 


2  (9^0  (9^0 

dril  drjl 


><90 

^dr]i 


ae(3(5  +  2)e^CA“^ 


(9T 

dri2 


0 


(3.27) 


At  this  stage  it  is  assumed  that  0,  0,  and  T  can  be  represented  as  a  power 
series  expansion  in  e. 


OO 


0(a;,  y)  =  0(er/i,  er/2)  =  X]  e  2  0m(hi,  V2) 

m=0 

(3.28) 

CXD 

0(a;,  y)  =  0(er/i,  er/2)  =  X]  e  2  0^(r/i,  7/2) 

m=0 

(3.29) 

00 

T{x,y)  =  T(er/i,er/2)  =  ^  e^^T^(r/i,  r/2) 

m=0 

(3.30) 
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As  given  by  Jin  and  Noda  [9]  the  solution  for  the  temperature  field  is  obtained 
by  ensuring  the  heat  flux  or  the  derivative  of  the  temperature  field  near  the  crack 
tip  is  singular.  This  means  that  the  potentials  for  temperature  are  one  power  less 
than  the  displacement  potentials.  Substituting  Eqs.  3.28-3.30  into  Eqs.3.26  and 
3.27  gives  the  following  equations. 


E,  m+. 

(e  2 


m+3  /  2  <9^0™  <9^0 


m=0 


a 


+ 


drjf  drj2 

(3J  +  2)  /  m-|-5 


(J  -h  2) 


—  e  2 


771  +  7 


771  +  5  /  d(f)' 


C 


C 

dr]i  '5  +  2  dr]2 


+ 


Tm  +  e  2  (^2/3  +  C)  A 


771  +  9 

-L  1^0 


i9?7i 


+  62  Tm  j  )  —  0 

(3.31) 


Oe(3i5  +  2)e"f^CA“' =  0  (3.32) 

0712 

For  Eqs. 3. 31  and  3.32  to  be  valid,  the  partial  differential  equations  corresponding 
to  each  power  of  e  ...)  should  vanish  independently.  This  leads  to  the 

following  set  of  partial  differential  equations. 

For  m  =  0  and  m  =  1 


a, 


+ 


dill  dr]. 


=  0 


(3.33) 


For  m  =  2 


!<9^0rr 
’  dril 


+ 


d^jjm 

dril 


(3.34) 


2*9^2  <9^2 

dril  drjl 


-c 


<900  1  000  \ 

dr]i  5  +  2  dri2  J 


+ 


35  +  2 
5  +  2 


CTcTo 


(3.35) 
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_  /  (9^/>o  1  <900 

^  1  r  I  o  c»„ 


^  (9?7f  d7]l  \  (9?7i  5  +  2  drj2 


(3.36) 


Eqs.  3.33  and  3.34  are  similar  to  that  for  homogeneous  material  where  as  the 
partial  differential  equations  Eqs. 3. 35  and  3.36,  associated  with  higher  powers  of 
are  coupled  to  the  differentials  of  the  lower  order  functions  through  the  nonho¬ 
mogeneity  parameters  and  temperature  term.  Eqs.  3.33  and  3.34  (i.e.  for  m  =  0 
and  m  =  1)  can  be  easily  reduced  to  Laplace’s  equation  in  the  respective  complex 
domains  0  =  hi  +  C  =  hi  +  ic(sV2  ,  i  =  and  the  solutions  are  same  as 

that  for  homogenous  material  [18-21]  and  can  be  written  as 


M  (  a  A  (m+3)/2  I  r'  (m+3)/2  .  +  /Q  Q7\ 

4>m[ph0ut)  =  AraPl  COS  I  - - - 6';  1  Cmp]  SlU  I  - - - Oi  1  (3.37) 


,m  +  3) 
2 


i^mips:  Os,  t)  =  siu  COS 


(3.38) 


where 


Pi  =  {pI  +  ,  tan(0i)  =  p,  =  {pi  +  alplY'^  ,  tan(0,)  =  ^ 


,2  I  „2_2\1/2 


and  Am,  Bm,  Cm,  and  Dm  are  real  constants. 

Using  the  dehnition  of  dynamic  stress  intensity  factor  Kjo  and  Kud  for  opening 
mode  and  shear  mode  [21],  the  relation  between  Ao  and  K/b  and  Co  and  Kud  are 


obtained. 


An  — 


4(1  -|-  al)  Kid 

Aaias  -  (1  +  alY  pv^’ 


(3.39) 


C'o  = 


3(4q;zq;s  -  (1  +  pv^’ 


(3.40) 
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where  /i  is  the  shear  modulus  of  the  material  at  the  propagating  crack-tip,  Ki£)  is 
mode-I  dynamic  stress  intensity  factor  and  Kuu  mode-II  dynamic  stress  intensity 
factor.  Now  considering  the  crack  face  boundary  conditions  (J22  =  0  and  (T12  =  0 
we  can  also  obtain  the  following  relationship  between  and  Bo  and  Co  and  Do 

Bo  =  (3.41) 

1  +  af 


1  -|- 

Do  =  -^Co  (3.42) 

2as 

The  solution  for  the  Eqs.  3.35  and  3.36  corresponding  to  higher  powers  of  e  (m  =  2) 
consists  of  two  parts:  solution  for  homogeneous  equation  and  a  particular  solution 
due  to  nonhomogeneity  and  temperature  and  these  can  be  obtained  recursively  [18, 
21].  The  solutions  02  and  02  obtained  are  given  below.  The  solutions  0oi  0Oi  0i  and 
01  automatically  satisfy  the  compatibility  equations  because  these  the  solutions 
are  same  as  those  for  homogeneous  materials.  Since  the  non-homogeneous  specihc 
parts  of  02  and  02  are  obtained  from  0o  and  0o,  they  also  automatically  satisfy 
the  compatibility  equations. 


A2Pi^^  cos  sin 

~  (^klocos  +  (Co  sin 


2^^  (bqCos  +  (T'osin 

l  ^s)  \  \^/  \^// 


5(S  +  2)(a^-ar‘ 

4(3i  +  2)  a, 

TslJTWW^)  2 


(3.43) 


02  =  B2pl^‘^  cos  +  D2pl^‘^  sin 

-  \  (^5osin  +  (Do  cos 

■  (i"')  +  (Co An  (0,))  (3.44) 
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By  assembling  together  the  above  result  for  the  first  few  terms  and  by  transforming 
back  to  the  x  —  y  plane,  the  combined  solution  can  be  written  for  0  and  ip  as  Eq. 
3.45  and  3.46 


0  =  Aorf^  cos  sin  cos  {29i)  +  Cirf  sin  (26'0 


+  A2r0  cos  +  C2ri  sin  -Oi 


(^^ocos  +  (Co  sin 

4(35  +  2)  a,  _5/2„:  5, 


\^ou  -r  ^  U.C  5/2  .  o 

l5  (5  +  2)  ""2^ 


(3.45) 


+  Bir‘1  sin  {29s)  +  Dir^  cos  {29s) 

+  B2r^J^  cos  +  D2rl^‘^  sin 

-  (^5osin  +  (Do cos 

-  ^5  /  2^^  ( Aq  sin  (^9^^  +  (Co  sin  ^  (3.46) 

5  {af-aj)  \  \2  J  J  J 

where 


2„,2\1/2 


n  = 


2„,2\V2 


tan(6'0  = - ,  Vs  =  +  a‘ly‘^) 

X 


tan{9s)  =  — 

X 


The  above  dehnitions  of  the  displacement  potentials  are  now  used  with  Eq.  3.6  to 
get  the  displacements  helds.  These  displacement  helds  are  then  used  to  get  strain 
helds.  These  strain  fields  and  Eq.  3.25  are  substituted  into  Eqs.  3.3  -  3.5  to  obtain 
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in-plane  stress  fields  around  the  crack  tip.  They  are  shown  in  Eq.  3.47  -  3.50. 


G  -r 


exp{Cx)i^o 


—  ac{36  +  sin  -|-  gircos(6')  -|-  q2r^^‘^ 

+  ^go/5r='/"sin(^0) 


sin 


(3.47) 
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avY 


exp{Cx)fio 


-r-i  ^/^cos  l^jJ  (S-a^(S  + 

{ “S  (?')('*-  “?  ('5  + 


Vl6a?  ' 

/3  C  1/2 

— - ^r.  ( 


3  -1/2  . 

-r,  sm 

4  ' 


/3  C  1/2  .  /30z' 

\lQaf 


{5 -of  (5 +  2)))  Co 


(l)  ('^  +  «'('^  +  2)))  C'o 


/  <3  -1/2  I  \  3  C  1/2 

—  cis  COS  —  —  - — rj  cos 


o  -1/2  •  I  \  ^  1/2  •  / 

-r„  '  cisSin  —  —  - — rJ  sm  — — 


Sa^ 

3  C 


8  as 

4Ai  [S  -  af  (5  +  2))  -  45ia, 


^y^cos  (  ^  J  {S-a‘f{S  + 


'8ds 
2 
OOs 


B, 

a 


0 


^2  -  (  —rl^^as  cos 
C2  +  (  ^r]/‘^as  sin 


^25  +  2af  —  cig 
[85  +  2  ac 


6  +  2  af 


+  ac{36  +  2)  {qor^^'^  sin  +  qir  cos(6')  +  q2r^^‘^  '• 
+  sin 
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3.4  Results  and  Discussions 

The  thermo-elasticity  relations  developed  above  are  next  used  to  study  the 
effect  of  temperature  held,  crack  velocity  and  the  non-homogeneity  on  the  variation 
of  maximum  shear  stress,  the  circumferential  stress  and  the  strain  energy  density 
around  the  crack-tip.  The  coefficient  of  thermal  expansion  (oc)  of  titanium  (which 
is  of  primary  interest  in  our  experimental  research)  ,  8.9  X  10  ®/°C  is  used  in  the 
analysis.  The  temperature  coefficient  Qo  is  varied  from  0  to  2000  .  This 

variation  results  in  nominal  temperature  change  of  about  100  °C.  It  is  assumed 
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that  the  resulting  temperature  range  generates  elastic  deformation  around  the 
crack-tip. 

To  plot  the  stresses  described  in  the  next  section,  the  following  fracture  pa¬ 
rameters  are  used:  crack  tip  velocity  c/cs=0.5,  Ki]^=lQQ  MPa,  Kiij:,=Q.2Kjj:,  MPa 
and  keff=\/KiD~+TQiD 

3.4.1  Variation  of  Stress  Components  Near  the  Crack  Tip 

The  variation  of  stress  field  near  the  crack  tip  for  the  given  thermo-mechanical 
conditions  can  be  represented  in  many  different  ways.  In  this  paper,  three  major 
stress  components;  circumferential  tensile  stress,  maximum  principal  stress  and 
maximum  shear  stress  are  plotted  as  a  function  of  angle  around  the  crack-tip.  The 
variation  of  maximum  shear  stress  (Tmax)  for  a  typical  hxed  value  of  Kef f/V 27rr 
and  constant  crack  tip  velocity  (c/c^  =  0.5)  is  plotted  as  shown  in  Fig.  3.2.  is 
plotted  for  different  temperature  helds  around  the  crack-tip.  The  angle  at  which 
maximum  shear  stress  occurs  changes  with  temperature  and  non  homogeneity  pa¬ 
rameter.  Angular  variation  of  the  maximum  principal  stress  and  circumferential 
tensile  stress  with  hxed  value  of  and  constant  crack  tip  velocity  {c/cg  =  0.5)  are 
plotted  in  Fig.  3.3  and  Fig.  3.4.  Similar  trend  as  seen  for  maximum  shear  stress 
is  observed  in  both  the  hgures  and  the  values  and  the  angle  at  which  the  max- 
imums  occur  change  with  temperature  and  non-homogeneity  parameters.  The 
results  obtained  in  the  present  study  collapse  to  the  ones  given  by  Freund  [18]  for 
homogeneous  materials  without  temperature  helds  (go  =  0)  if  mode-I  and  mode-II 
helds  are  separated  and  plotted. 

3.4.2  Crack  Extension  Angle 

A  dynamically  moving  crack  tends  to  deviate  from  its  path  due  to  crack  tip 
instability  conditions.  In  the  present  study,  using  the  derived  thermo-mechanical 
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a)^  =  0 


b)^=0.4  c)<'=-0.4 

Figure  3.2.  Variation  of  normalized  maximum  shear  stress  with  angle  around  the 
crack-tip  for  mixed  mode  thermo-mechanical  loading  in  a  FGM  for  several  values 
of  temperature  coefficients  {Kju:,/ Ki^  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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a^/  (Kg^^/V(2rt')) 


b)^=0.4  c)^=-0.4 


Figure  3.3.  Variation  of  normalized  the  largest  principal  stress  with  angle  around 
the  crack-tip  for  mixed  mode  thermo-mechanical  loading  in  a  FGM  for  several 
values  of  temperature  coefficients  {Kjjd/ Kjd  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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0  (degree) 


a)  ^=0 


b)^=0.4  c)^=-0.4 


Figure  3.4.  Variation  of  normalized  circumferential  stress  with  angle  around  the 
crack-tip  for  mixed  mode  thermo-mechanical  loading  in  a  FGM  for  several  values 
of  temperature  coefficients  {Kjjd/Kid  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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stress  field  equations,  the  effect  of  temperature,  crack-tip  velocity  and  material  non¬ 
homogeneity  on  the  crack  instability  is  presented.  The  theoretical  prediction  of 
crack  extension  angle  is  investigated  by  using  the  two  well-known  fracture  criterias: 
minimum  strain  energy  density  (S-criterion)  and  maximum  circumferential  stress 
{aee  -criterion). 

Minimum  strain- energy  density  (MSED)  eriterion 

According  to  this  criterion  [21],  the  crack  initiates  when  the  strain  energy  density 
achieves  a  critical  value  and  propagates  in  the  direction  of  minimum  strain-energy 
density  value.  The  strain  energy  density  dW/dV  near  the  crack  tip  for  an  FGM  is 
given  as 


dW 

W 


S  = 


d/ie'i* 


((1  -  ^)  +  ^ly 


)  —  2z/(T, 


xx(^yy  +  ^xy) 


(3.50) 


Fracture  takes  place  in  the  direction  of  minimum  S,  and  the  condition  can  be 
obtained  by  using  Eq.(37) 


dS 


=  0, 


&^S 

W 


>  0 


at  S  =  Sr 


(3.51) 


where  Sc  is  the  critical  strain  energy  density 

Variations  of  strain  energy  density  with  angle  6  from  — tt  to  vr  around  the  crack- 
tip  for  mixed  mode  thermo-mechanical  loading  in  an  FGM  for  several  values  of 
the  temperature  coefficient  are  plotted  in  Fig.  3.5.  The  angle  at  which  the  strain 
energy  density  reaches  a  minimum  value  changes  with  temperature  and  the  non¬ 
homogeneity  parameter. 

Maximum  circumferential-stress  (MCS)  criterion 
The  maximum  circumferential  stress  criterion  [22]  states  that,  crack  growth  will 
occur  in  the  direction  of  the  maximum  circumferential  stress  and  will  take  place 
when  the  maximum  circumferential  stress  reaches  a  critical  value,  and  it  can  be 
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0  (degree)  6  (degree) 


a)^=0  b) ^  =  0  ,  enlarged  view 


0  (degree)  0  (degree) 


c)^=0.4  d)^=-0.4 

Figure  3.5.  Variation  of  strain  energy  density  with  angle  around  the  crack-tip  for 
mixed  mode  thermo-mechanical  loading  in  FGM  subjected  to  different  temperature 
helds  (Kjju/Kjd  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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given  as  Eq.3.52 


daee 

de 


d'^aee 

de^ 


<  0 


(J00  —  {U00)c 


(3.52) 


where  {(J00)c  is  the  critical  circumferential  stress 

Variations  of  circumferential  stress  with  angle  9  around  the  crack-tip  for  mixed 
mode  thermo-mechanical  loading  in  an  FGM  for  several  values  of  temperature  co¬ 
efficients  are  plotted  in  Fig.  3.6.  The  angle  at  which  the  circumferential  stress 
reaches  a  maximum  value  changes  with  temperature  and  non-homogeneity  param¬ 
eter. 

Based  on  the  above  two  criteria’s,  the  effects  of  velocity,  non-homogeneity  and 
temperature  on  the  crack  extension  angle  {9)  are  further  investigated. 


Effect  of  Crack  Tip  Velocity 

The  crack  extension  angles  as  a  function  of  crack  tip  velocities  as  predicted 
by  the  above  two  criterias  are  shown  in  Fig.  3.7.  For  pure  mode-I  loading 
{Kiid/Kid  =  0),  the  crack  extends  along  until  the  crack  tip  velocity  reaches  a 
critical  value  at  which  instability  occurs  [23] .  When  the  crack  tip  velocity  reaches 
the  critical  value,  the  crack  deviates  and  extends  to  a  different  angle.  For  example 
at  a  crack  tip  velocity  of  c/cg  =  0.7,  the  MSFD  criterion  predicts  a  crack  extension 
angle  of  about  —55°  and  the  MCS  criterion  predicts  about  —38°.  On  the  other 
hand,  as  the  value  of  Kuo/Kid  increases  from  0  to  1  and  later  from  1  to  cx)  the 
crack  extension  angle  increases  monotonically.  Broek  [24],  in  his  book  gives  the 
crack  extension  angles  for  mixed  mode  quasi  static  loading  and  these  results  match 
well  with  the  predictions  from  the  current  study. 


Effect  of  Non-homogeneity 

The  effect  of  non-homogeneity  parameter  on  crack  extension  angle  for  a  crack 
tip  velocity  of  0.5cs  at  room  temperature  (go  =  0)  is  shown  in  Fig.  3.8.  For  both 
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c)^  =  0.4  d)^  =  -0.4 

Figure  3.6.  Variation  of  normalized  circumferential  stress  with  angle  around  the 
crack-tip  for  mixed  mode  thermo-mechanical  loading  in  FGM  subjected  to  different 
temperature  helds  {Kn^,/ Kjd  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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(a)  Minimum  strain  energy  density  criterion 


(b)  Maximum  circumferential  stress  criterion 


Figure  3.7.  Crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed 
mode  thermo-mechanical  loading  in  homogeneous  material  {(  =0,  r=0.002m). 
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(a)  Minimum  strain  energy  density  criterion 


(b)  Maximum  circumferential  stress  criterion 


Figure  3.8.  Crack  extension  angle  as  a  function  of  non-homogeneity  parameter  for 
mixed  mode  crack  without  heat  (c/c5=0.5,  r=0.002m). 
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homogenous  material  (i.e.  C  =  0  )  a  FGM  with  increasing  stiffness  in  the 
direction  of  crack  growth  (i.e.  C  >  0);  fhe  MCS  criterion  provides  a  maximum 
value  and  the  MSED  criterion  provides  a  minimum  value  along  6  =  0  under  pure 
mode-I  loading  {Km:,/ Kid  =  0).  However,  for  a  FGM  with  decreasing  stiffness  in 
the  direction  of  crack  growth  (i.e.  C  <  0  ),  the  MGS  criterion  provides  a  maximum 
value  at  angle  of  about  —55"  and  the  MSED  criterion  provides  minimum  value  at 
angle  of  about  —28°.  It  can  be  also  observed  that,  for  complete  range  of  Kud/ Kjd 
values,  a  FGM  with  a  <  0  has  larger  crack  extension  angle  compared  to  both 
homogenous  and  a  FGM  with  a  >  0.  This  might  be  attributed  to  presence  of 
compliant  material  ahead  of  the  crack  tip  and  the  propagating  crack  needs  less 
energy  to  be  unstable. 

Effect  of  Temperature 

The  effect  of  temperature  field  on  the  crack  extension  angle  for  a  crack  tip 
velocity  of  0.5cs  is  shown  in  Fig.  3.9.  Both  the  criterion  show  that,  for  homogeneous 
material,  the  crack  extension  direction  at  room  temperature  is  along  9  =  —30°  and 
the  value  decreases  slowly  with  increase  in  applied  temperature  held.  For  FGM 
with  a  >  0,  the  crack  extension  angle  is  along  9  =  —15°  at  room  temperature  and 
again  the  value  decreases  with  increase  in  applied  temperature  held.  In  the  case  of 
FGM  with  a  <  0,  the  crack  extension  angle  is  about  —50°  at  room  temperature  and 
increases  in  magnitude  as  the  temperature  increases.  The  increased  temperature 
held  increases  the  compliance  of  the  already  compliant  material  (in  case  ofa  <  0) 
and  hence  the  propagating  crack  needs  less  energy  to  cause  instability. 

Effect  of  Crack  Tip  Velocity,  Temperature  aud  Nou-homogeueity 

The  ehects  of  temperature  on  the  crack  extension  angle  as  a  function  of  crack 
tip  velocity  in  homogeneous  material,  a  FGM  with  a  >  0  and  a  FGM  with  a  <  0 
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Figure  3.9.  Effect  of  temperature  on  the  crack  extension  angle  for  mixed  mode 
loading  in  FGM  (Kjjd/Kjd  =  0.2,  c/cg  =  0.5,  r  =  0.002m). 
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are  plotted  in  Figs.  3.10,  3.11  and  3.12.  Figure  3.10  represents  crack  extension 
angle  as  a  function  of  crack  tip  velocity  for  mixed-mode  mechanical  loading  in 
homogeneous  material  at  room  and  elevated  temperatures.  For  a  crack  tip  velocity 
of  cjcg  =  0.3,  both  the  criterions  predict  decrease  in  crack  extension  angle  for 
increasing  temperature  helds.  Furthermore,  for  c/c^  >  0.3,  both  criterions  initially 
predict  decrease  in  crack  extension  angle  for  increasing  temperature  held  and  show 
insignihcant  change  at  higher  velocities. 

Figure  3.11  represents  crack  extension  angle  as  a  function  of  crack-tip  velocity 
at  room  and  at  elevated  temperatures  for  mixed  mode  mechanical  loading  in  FGM 
with  a  >  0.  For  a  crack  tip  velocity  of  c/c^  =  0.3,  the  MSED  criterion  predicts  a 
single  value  of  6^  =  0°.  However  for  the  same  crack  tip  velocity,  the  MCS  criterion 
predicts  different  crack  extension  angles  for  different  temperature  helds.  Again 
similar  to  Fig.  3.10,  for  c/c^  >  0.3,  FGM  with  a  >  0  the  crack  extension  angle 
decreases  as  the  temperature  increases.  However,  as  the  crack  velocity  increases 
the  diherence  in  crack  extension  angles  becomes  less  and  less. 

Figure  3.12  represents  crack  extension  angle  as  a  function  of  crack-tip  velocity 
in  the  direction  of  crack  growth  at  room  and  at  elevated  temperatures  for  mixed 
mode  mechanical  loading  in  FGM  with  a  <  0  .  The  variation  of  crack  extension 
angle  in  Fig.  3.12  shows  signihcantly  diherent  trend  compared  to  the  above  two 
cases.  Both  criterions  predict  diherent  crack  extension  angles  at  cjcg  =  0.3  for 
diherent  temperature  helds.  As  the  velocity  increases,  the  crack  extension  angle 
initially  decreases  and  then  increases.  It  is  quite  contrasting  to  above  two  cases 
that  for  any  crack  tip  velocity,  higher  the  temperature  higher  the  crack  extension 
angle.  The  diherent  trend  predicted  for  FGM  with  a  <  0  might  be  again  attributed 
to  compliant  material  ahead  of  the  propagating  crack  tip. 

Previous  work  in  literature  states  that  the  critical  velocity  at  which  the  crack 
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(a)  Minimum  strain  energy  density  criterion 


(b)  Maximum  circumferential  stress  criterion 

Figure  3.10.  Crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed 
mode  thermo-mechanical  loading  for  homogenous  material  {Kjjd/ Kjd  =  0.2,  (  = 
0,  r  =  0.002m). 
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(a)  Minimum  strain  energy  density  criterion 


Figure  3.11.  Crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed 
mode  thermo-mechanical  loading  at  different  temperature  {Kud/Kid  =  0.2,  (  = 
0.4,  r  =  0.002m). 
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(a)  Minimum  strain  energy  density  criterion 


(b)  Maximum  circumferential  stress  criterion 

Figure  3.12.  Crack  extension  angle  as  a  function  of  crack  tip  velocity  for  mixed 
mode  thermo-mechanical  loading  at  different  temperature  {Kud/Kid  =  0.2,  (  = 
—0.4,  r  =  0.002m). 
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instability  starts  in  rnnning  cracks  depends  on  the  extent  of  non-singnlar  stress 
fields  at  the  crack-tip  [18].  It  is  noted  from  the  derived  equations  that  stress-held 
terms  associated  with  non-  homogeneity  and  temperature  held  are  non-singular 
and  this  could  be  the  reason  why  the  crack  extension  angle  is  diherent  for  diherent 
non-homogeneity  parameters  and  temperatures.  However,  the  non-homogeneity 
parameter  has  more  signihcant  role  in  the  crack  extension  angle  compared  to  the 
ehect  of  temperature  held. 

3.5  Summary 

The  stress-helds  near  the  crack  tip  for  mixed  mode  thermo-mechanical  load¬ 
ing  in  graded  material  are  developed  using  displacement  potentials  in  conjugation 
with  an  asymptotic  approach.  Using  insulated  crack  face  boundary  condition  and 
steady  state  heat  conduction  assumption,  hrst  the  temperature  held  near  to  the 
crack  tip  is  developed.  By  incorporating  the  developed  temperature  held  into  the 
mechanical  held,  the  thermo-mechanical  stress  helds  near  to  the  crack  tip  are  devel¬ 
oped.  Using  the  developed  equations,  angular  variation  of  maximum  shear  stress, 
circumferential  stress  and  the  largest  principal  stress  are  plotted  as  a  function  of 
temperature  around  the  crack-tip.  Using  both  minimum  strain  energy  density  cri¬ 
terion  and  maximum  circumferential  stress  criterion,  the  crack  instability  direction 
for  various  crack-tip  speeds  and  non-homogeneous  coefficients  is  also  determined. 
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CHAPTER  4 


Dynamic  Constitutive  Behavior  of  Ti/TiB  FGM  under  Thermo- 

Mechanical  Loading 

4.1  Abstract 

An  experimental  investigation  is  conducted  to  evaluate  the  thermo-mechanical 
constitutive  behavior  of  a  functionally  graded  material  (FGM)  under  dynamic 
loading.  Cylindrical  specimens  are  machined  from  titanium  /  titanium  mono¬ 
boride  (Ti/TiB)  layered  FGM  plate  using  electrical  discharge  machining  (EDM). 
A  Split  Hopkinson  Pressure  Bar  (SHPB)  apparatus  with  infrared  spot  heaters  is 
used  to  investigate  the  effect  of  temperature  on  mechanical  response  of  the  FGM 
material.  A  series  of  experiments  are  conducted  at  different  temperatures  and 
the  stress  strain  relation  for  different  temperatures  is  obtained.  The  material 
showed  high  thermal  softening  at  elevated  temperature  resulting  in  a  reduction  in 
compressive  strength  and  an  increase  in  failure  strain. 

4.2  Introduction 

FGMs  are  materials,  which  have  continuously  or  discreetly  varying  mechanical 
and  thermal  properties  in  a  specific  direction.  Ideally,  by  grading  the  composition 
from  one  surface  to  another,  it  is  possible  to  create  a  material  suitable  for  a  par¬ 
ticular  application.  This  type  of  material  was  first  proposed  in  Japan  in  1987  for 
space  access  vehicles  that  can  operate  reliably  in  combined  thermal  and  mechan¬ 
ical  loadings.  Typically,  FGMs  are  made  of  a  metal  and  a  ceramic  as  opposite 
faces  with  the  intermediate  zones  consisting  of  varying  volume  fractions  of  con¬ 
stituents.  Titanium  /  titanium  mono-boride  (Ti/TiB)  is  one  of  the  FGM  which 
is  currently  in  use.  This  material  consists  of  the  desirable  properties  of  ceramic, 
such  as  hardness,  corrosion  resistance,  and  high  melting  temperature,  without  los- 
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ing  the  required  properties  of  metallic  titanium,  such  as  good  fracture  toughness, 
machinability  and  weldability. 

For  efficient  design  of  structures  using  FGM  materials  a  fundamental  under¬ 
standing  of  the  properties  of  FGM  materials  at  different  loading  conditions  is 
required.  To  date  research  on  FGM  materials  has  focused  mostly  on  quasi-static 
problems  with  very  few  studies  in  the  dynamic  regime.  Also,  most  of  the  mechanics 
studies  relate  to  crack  problems.  Delale  and  Erdogan  [1],  Eischen  [2]  and  Jin  and 
Noda  [3]  solved  crack  problems  for  non-homogeneous  materials  under  quasi-static 
mechanical  loading.  All  these  investigations  concluded  that  the  inverse-square 
root  singularity  at  the  crack  is  not  affected  by  nonhomogeneity.  Eor  propagating 
cracks  in  EGMs,  Parameswaran  and  Shukla  [4]  and  Ghalivendra  et  ah  [5]  devel¬ 
oped  the  structure  of  the  first  stress  invariant  and  the  out  of  plane  displacement.  In 
their  study  they  brought  out  the  effects  of  nonhomogeniety  through  an  asymptotic 
analysis.  Lee  [6]  developed  nonhomogeneity  specific  terms  for  individual  stress 
and  displacement  components  using  displacement  potentials.  Recently,  Shukla 
and  Jain  [7]  and  Ghalivendra  and  Shukla  [8]  developed  transient  held  eqnations 
for  cracks  propagating  at  arbitrary  velocities.  Ghalivendra  [9]  developed  an  asymp¬ 
totic  analysis  of  the  transient  out  of  plane  displacement  helds  for  a  curved  crack 
propagating  at  arbitrary  velocity  in  EGMs.  A  review  paper  by  Shnkla  et  al.  [10] 
presents  a  comprehensive  summary  of  dynamic  fracture  studies  in  FGM’s. 

There  are  few  stndies  related  to  the  stress  and  displacement  helds  dne  to  ap¬ 
plied  loads  in  graded  materials,  Eor  instance  Wang  et  al.  [11]  present  the  solution 
for  displacements  and  stresses  in  an  FGM  snbjected  to  a  vertical  point  load  in 
a  continnously  inhomogeneous  transversely  isotropic  half-space  with  Yonng’s  and 
shear  modnli  varying  exponentially  with  depth.  Horgan  and  Ghan  [12]  investigate 
the  ehects  of  material  inhomogeneity  on  the  response  of  linearly  elastic  isotropic 
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hollow  circular  cylinders  under  uniform  internal  or  external  pressure.  Li  et  al.  [13] 
examined  the  mechanical  behavior  of  layered  plates  made  of  metal-ceramic  com¬ 
posites  with  the  volume  fraction  of  ceramic  reinforcement  varying  through  thick¬ 
ness  direction  under  impulse  loading.  Chi  and  Chung  [14]  studies  the  mechanical 
behavior  of  functionally  graded  material  plate  under  transverse  loading. 

There  are  very  few  studies  on  the  properties  of  FGM  under  thermal  or  thermo¬ 
mechanical  loading.  Jain  et  al.  [15]  developed  the  stress  held  equations  for  quasi 
static  cracks  under  thermo-mechanical  loading  in  FGM’s.  Praveen  and  Reddy  [16] 
carried  out  the  nonlinear  transient  analysis  of  stress  and  dehection  for  a  function¬ 
ally  graded  ceramic-metal  plate  under  thermal  loading  by  using  the  hnite  element 
method  (FEM).  They  found  that  the  response  of  the  plates  with  material  properties 
between  those  of  the  ceramic  and  metal  is  not  intermediate  to  the  response  of  the 
ceramic  and  metal  plates.  Dai  et  al.  [17]  also  presented  thermo-mechanical  anal¬ 
ysis  of  dehections  under  different  loadings  in  FGM  plates  containing  distributed 
piezoelectric  sensors  and  actuators  using  element-free  Galerkin  method. 

In  this  paper  an  experimental  investigation  on  the  dynamic  constitutive  be¬ 
havior  of  Ti/TiB  FGM  under  thermo- mechanical  loading  is  presented.  SHPB 
apparatus  with  infrared  spot  heating  system  is  used  to  investigate  the  constitutive 
properties  of  the  FGM  at  different  temperatures.  The  material  exhibited  thermal 
softening  at  higher  temperatures,  which  results  in  a  decrease  in  flow  stress  and  an 
increase  in  failure  strain.  The  effect  of  machining  conditions  on  the  constitutive 
behavior  is  also  demonstrated  in  these  experiments. 

4.3  Material  and  Specimen  Geometry 

The  layered  (Ti/TiB)  functionally  graded  material  used  in  this  study  is  sup¬ 
plied  by  BAE  Systems  in  the  form  3.175  mm  thick  plates.  Hill  et  al.  [18]  have 
explained  the  fabrication  technique  for  this  material  and  their  procedure  is  pre- 
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sented  below.  Tape  cast  layer  composed  of  varying  mixtures  of  titanium  and 
titanium  diboride  powders  are  placed  on  top  of  commercially  pure  titanium  metal 
plate.  The  assembled  laminate  is  hot  pressed  at  1578  K  at  a  pressure  of  12.8  MPa. 
In  order  to  facilitate  densihcation  at  this  temperature,  a  proprietary  sintering  aid- 
containing  nickel  is  added  to  the  starting  powders.  This  material  creates  a  liquid 
phase  at  1215  K  that  also  catalyzes  the  reaction  of  titanium  and  TiB2  to  form  TiB 
with  virtually  no  residual  TiB2.  The  resulting  FGM  is  composed  of  seven  layers 
ranging  from  pure  Ti  on  one  side  to  85%  TiB  on  the  other  [18].  The  material  used 
in  our  studies  shows  no  clear  or  distinct  interface  between  the  layers  as  shown  in 
Fig.  4.1.  Table  4.1  shows  the  composition  and  the  physical  properties  of  each  layer 
of  the  FGM  as  provided  by  the  vendor  and  [18]. 

■  '0*oTi 

30**TiB 

85*  0  Ti 
15* .  TiB 

100«.  Ti 


X2  7 


Figure  4.1.  Scanning  electron  microscope  image  of  layer  1  (TiB  0%),  layer  2  (TiB 
15%)  and  layer  3  (TiB  30%). 

Gylindrical  specimens  with  diameter  of  7.62  mm  and  thickness  of  3.175  mm  are 
machined  from  the  FGM  plate  discussed  above.  An  electrical  discharge  machining 
(EDM)  technique  is  used  to  machine  the  specimens.  The  thickness  of  the  specimen 
is  selected  from  the  available  FGM  plate  thickness,  where  as  the  diameter  of  the 
specimen  is  selected  based  on  the  SHPB  theory  to  assure  uniform  loading  and 
constant  stain  rate.  The  hrst  batch  of  specimens  is  machined  using  wire  EDM  and 
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Table  4.1.  Composition  and  elastic  properties  of  Ti/TiB  FGM 


Layer  no. 

Vol.  %Ti 

Vol.  % 

TiB 

Thickness 

(mm) 

E  (GPa) 

V 

1 

100 

0 

0.2032 

106 

0.340 

2 

85 

15 

0.3810 

170 

0.278 

3 

70 

30 

0.3810 

227 

0.238 

4 

55 

45 

0.3810 

262 

- 

5 

40 

60 

0.3810 

289 

- 

6 

25 

75 

0.3810 

303 

0.152 

7 

15 

85 

1.0668 

316 

0.140 

cooled  by  flooding  (Fig.  4.2a).  The  second  batch  of  specimens  is  machined  with 
wire  EDM  with  the  material  fully  immersed  in  coolant  (Fig.  4.2b). 


Figure  4.2.  Wire  EDM  machined  specimens  (a)  cooled  by  flooding,  (b)  cooled  by 
fully  immersing  in  coolant. 


4.4  Experimental  Procedure 

SHPB  apparatus  with  infrared  spot  heaters  is  used  to  investigate  the  dynamic 
behavior  of  Ti/TiB  FGM  at  room  and  at  elevated  temperatures.  In  the  present 
study,  the  technique  used  by  Lennon  and  Ramesh,  [19]  and  Shazly  et  ah  [20]  is 
adapted.  The  schematic  description  of  the  SHPB  in  our  lab  is  shown  in  Fig.  4.3 
and  the  set-up  with  infrared  (IR)  spot  heaters  is  shown  in  Fig.  4.4.  The  IR  spot 
heaters  have  a  circular  cut  shields  that  can  concentrate  a  heat  flux  as  high  as  650 
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watts  per  6.45  square  centimeters  at  the  focal  point  with  a  6.35  mm  diameter. 
The  bars  are  made  up  of  12.7  mm  diameter  marraging  steel  having  nominal  yield 
strength  of  2500  MPa.  The  striker  bar  is  101.6  mm  long,  while  the  incident  and 
transmitter  bars  are  1.27  m  long.  The  striker  bar  is  propelled  using  an  air-operated 
gun.  Two  strain  gages  are  placed  on  both  the  incident  and  the  transmitter  bar  at 
equal  distance  form  the  bar-specimen  assembly.  These  gages  are  placed  at  1800  to 
each  other  to  avoid  the  bending  effect  on  the  strain  data.  The  strain  signals  are 
recorded  using  a  Vishay  2301A  signal-conditioning  ampliher  that  connects  with  an 
oscilloscope.  The  specimen  is  sandwiched  between  the  two  bars.  However,  during 
the  elevated  temperature  testing  two-carbide  inserts  are  placed  between  the  two 
bars  and  the  specimen  is  sandwiched  between  the  inserts.  The  inserts  are  used  to 
eliminate  the  temperature  gradient  in  the  bars.  The  inserts  are  impedance  matched 
to  the  bars  and  hence  do  not  disturb  the  incident,  transmitted  and  reflected  wave 
prohles.  The  insert  and  the  bar  are  placed  in  a  sleeve  to  assist  in  alignment. 
The  assembly  of  the  carbide  inserts  with  the  specimen  and  the  position  of  the 
thermocouple  is  shown  in  Fig.  4.5. 

For  elevated  temperature  experiment,  the  specimen  is  heated  up  to  the  de¬ 
sired  temperature  (usually  about  30  °C  -  40  °C  higher  than  the  test  temperature) 
and  the  bars  are  brought  manually  in  contact  with  the  inserts-specimen-inserts 
assembly.  The  temperature  of  the  specimen  is  monitored  by  0.127  mm  chromel- 
alumel  thermocouple,  which  is  attached  into  the  specimen  using  high  temperature 
epoxy.  Digital  camera  is  used  to  record  the  temperature  reading  and  to  monitor 
the  sequence  of  loading  on  the  specimen.  The  exact  temperature  on  the  specimen 
at  which  the  load  applied  is  found  form  the  record.  Boron  nitride  is  used  as  the 
lubrication  between  specimen-inserts  and  inserts-bar  assembly  for  all  elevated  tem- 
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Figure  4.3.  Schematic  of  split  Hopkinson  pressure  bar. 


Figure  4.4.  Split  Hopkinson  pressure  bar  with  heater  assembly. 
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Figure  4.5.  Specimen-carbide  inserts,  thermocouple  and  sleeves  assembly. 

perature  experiments.  In  most  of  the  experiments  it  takes  less  than  two  minutes 
to  heat  the  specimen  to  the  required  temperature.  It  takes  less  than  10  seconds  to 
bring  the  bars  in  contact  with  the  specimen  and  hre  the  gun. 

Using  the  above  conhguration  a  series  of  experiments  was  conducted.  For  all 
experiments  a  clay  pulse  shaper  is  placed  on  the  impact  face  of  the  incident  bar. 
The  pulse  shaper  smoothes  the  relative  sharp  front  of  the  incident  stress  wave, 
thus  allowing  high  strain-rate  experiments  to  be  conducted  at  near  constant  strain 
rates  [20].  This  also  helps  the  premature  failure  of  the  relatively  brittle  Ti/TiB 
FGM,  especially  during  the  early  part  of  stress  wave  loading  [20].  A  typical  strain 
record  from  one  of  the  experiment  is  shown  in  Fig.  4.6. 

4.5  Results  and  Discussion 
4.5.1  Stress-strain  Relation 

A  series  of  eight  experiments  is  conducted  at  four  different  temperatures  on 
specimens  machined  by  wire  EDM  and  cooled  by  flooding.  The  stress  strain  curves 
obtained  from  the  dynamic  experiments  with  strain  rate  of  about  3  ±  10^  per  sec 
are  shown  in  Fig.  4.7.  The  stress-strain  curves  show  a  monotonic  increase  till 
failure  is  reached.  In  all  these  experiments  failure  was  achieved.  The  unloading 
part  of  the  curves  after  failure  is  not  meaningful.  At  room  temperature  a  failure 
stress  of  2150  MPa  and  a  failure  strain  of  1.4  %  is  observed.  As  the  temperature 
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Figure  4.6.  Typical  strain  pulse  profiles  obtained  during  an  experiment 


increases  a  reduction  in  the  failure  stress  and  an  increase  in  the  failure  strain  is 


observed.  For  example  at  600  a  failure  stress  of  910  MPa  and  a  failure  strain 
of  2.6  %  is  observed.  In  all  the  room  temperature  experiments  the  specimens  fail 
completely  leaving  small  fragments  as  shown  in  Fig.  4.8.  It  was  also  observed 
that  the  fragmentation  process  even  at  room  temperature  was  associated  with  a 
flash  of  light.  A  typical  flash  of  light  during  failure  was  recorded  using  a  digital 
camera  and  is  shown  in  Fig.  4.9.  The  occurrence  of  flash  led  us  to  believe  that  high 
residual  stresses  were  present  in  the  material  and  these  could  have  been  created 
during  the  machining  process.  The  material  fabrication  technique  is  believed  to 
create  a  relatively  stress  free  material.  The  coloration  on  the  edge  of  the  specimen 
shown  in  Fig.  4.2a  indicates  a  heat-effected  zone  generated  during  machining  by 
wire  EDM. 


Figure  4.7.  True  stress- True  strain  curve  as  a  function  of  temperature  (specimen 
cooled  by  flooding  during  EDM). 

To  investigate  the  effect  of  machining  technique  on  the  behavior  of  the  ma¬ 
terial  another  batch  of  sample  was  machined  using  EDM  with  the  material  fully 


78 


2  mm 


Figure  4.8.  A  typical  failed  specimen  at  room  temperature  under  dynamic  loading. 


Figure  4.9.  A  flash  of  light  produced  during  room  temperature  dynamic  experiment 
of  specimen  machined  by  EDM  and  cooled  by  flooding. 
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immersed  in  the  coolant.  A  series  of  eight  experiments  was  again  conducted  using 
these  specimens  at  four  different  temperatures.  Figure  4.10  shows  the  dynamic 
compressive  stress  strain  curve  of  the  materials  at  different  temperatures.  At  room 
temperature  a  failure  stress  of  3050  MPa  and  a  failure  strain  Of  1.2  %  is  obtained. 
As  the  temperature  increases,  the  ductility  of  the  material  increases  and  results  in 
an  increase  in  the  failure  strain  and  a  decrease  in  failure  stress.  For  example,  at 
the  testing  temperature  of  730  °C  a  failure  stress  of  2300  MPa  and  a  failure  strain 
of  2.6  %  are  observed  and  at  800  °C  a  failure  stress  of  2100  MPa  and  a  failure 
strain  of  3.4  %  are  obtained.  There  was  no  flash  of  light  observed  during  testing  of 
these  specimens.  Figure  4.11  shows  the  effect  of  temperature  on  the  failure  stress 


Figure  4.10.  True  stress  -  True  stain  curve  as  a  function  of  temperature  (specimen 
fully  immersed  in  coolant  during  EDM). 

of  Ti/TiB  for  both  the  machining  cases  considered.  Both  the  curves  show  similar 
trend  that  the  failure  stress  decreases  as  the  temperature  increase.  In  the  case  of 
specimen  machined  by  EDM  and  cooled  by  flooding  the  failure  stress  at  600  °C 
is  reduced  by  54  %  compared  with  the  room  temperature  failure  stress.  Eor  fully 
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immersed  cooling  EDM  specimen  the  failure  stress  at  800  °C  is  reduced  by  30  % 
compared  with  the  room  temperature  failure  stress. 
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Figure  4.11.  Failure  stress  as  a  function  of  temperature. 

4.5.2  Post  Failure  Analysis 

Photographs  of  the  fracture  specimen  were  taken  for  post  failure  analysis.  In 
all  the  room  temperature  experiments  the  specimens  fractured  completely,  leaving 
small  fragments  as  shown  in  Fig.  4.8.  However  in  the  case  of  high  temperature  ex¬ 
periments  the  specimens  fracture  leaving  the  fragments  loosely  bonded  as  shown  in 
Fig.  4.12.  This  indicates  at  higher  temperature  the  material  is  thermally  softened 
and  shows  a  ductile  failure.  Figure  4.13  shows  a  micrograph  image  of  fragment 
from  a  room  temperature  experiment.  From  the  hgure  it  is  observed  that  the 
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cracks  run  predominantly  in  the  loading  direction,  which  is  perpendicular  to  the 
interface  of  layers,  indicating  the  bonding  of  the  interface  is  strong.  Scanning 


Figure  4.12.  Typical  failed  specimen  during  dynamic  loading  at  elevated  temper¬ 
ature. 

Electron  Microscopic  (SEM)  images  were  taken  of  the  fractured  specimens  tested 
at  four  different  temperatures.  Figure  14  shows  SEM  images  of  four  high  temper¬ 
ature  specimens  taken  at  a  magnihcation  of  1000.  The  post  mortem  examination 
of  the  failed  specimens  showed  that  at  room  temperature  and  at  460  °C  the  spec¬ 
imens  fracture  by  cleavage  (Fig.  4.14a  and  Fig.  4.14b).  The  fracture  surfaces  from 
specimens  tested  at  temperatures  of  730  °C  and  800  °C  are  shown  in  Fig.  4.14c 
and  Fig.  4.14d  respectively,  and  both  of  them  indicate  a  ductile  failure  mode.  A 
number  of  large  smooth  surfaces  observed  in  Fig.  4.14c  and  Fig.  4.14d  are  believed 
to  have  been  caused  by  plastic  shear  [21].  These  smooth  surfaces  are  shown  at 
higher  magnihcation  (5000  times)  in  Fig.  4.15a  and  Fig.  4.15b. 

4.6  Summary 

The  dynamic  constitutive  behavior  of  Ti/TiB  FGM  under  thermo-mechanical 
loading  is  investigated.  SHPB  apparatus  along  with  spot  heaters  is  utilized  to 
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Figure  4.13.  Typical  fracture  surface  of  Ti/TiB  specimen. 


study  the  dynamic  constitutive  behavior  of  the  material  in  the  temperature  range 
from  25  °C  to  800  °C.  The  hndings  are  summarized  and  presented  as  follows: 

•  There  are  no  distinct  interfaces  between  the  layers  of  the  FGM  material 
investigated.  The  bonds  between  the  layers  are  strong  and  the  cracks  run 
predominantly  perpendicular  to  the  layers  during  fracture. 

•  Samples  cut  by  wire  EDM  and  cooled  by  flooding  and  tested  at  25  °C  show  a 
failure  stress  of  2150  MPa  and  a  failure  strain  of  1.4  %.  At  600  °C  the  same 
sample  type  shows  a  failure  stress  of  910  MPa  and  a  failure  strain  of  2.6  %. 
This  represents  a  54  %  reduction  in  failure  stress  and  an  85  %  increases  in 
failure  strain. 

•  Fully  immersed  cooled  EDM  specimens  tested  at  25  °C  show  a  failure  stress 
of  3050  MPa  and  a  failure  strain  of  1.25  %.  At  800  °C  a  failure  stress  of 
2100  MPa  and  a  failure  strain  of  3.4  %  is  observed.  This  represents  a  30  % 
reduction  in  failure  stress  and  a  125  %  increase  in  failure  strain. 
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(c)T=730°C  (d)T=800°C 


Figure  4.14.  SEM  images,  at  a  magnification  of  1000,  of  fracture  surface  of  four 
different  specimens  tested  at  four  different  temperatnres. 
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(a)  T=7  30  °C  (b)  T=800  °C 


Figure  4.15.  SEM  images,  at  a  magnification  of  5000,  of  fracture  surface  of  speci¬ 
mens  tested  as  two  different  temperatures. 

•  The  FGM  material  showed  thermal  softening  at  higher  temperature  with 
a  decrease  in  failure  stress  and  an  increase  in  failure  strain  for  both  the 
machining  cases  considered. 

•  The  specimens  machined  from  fully  immersed  cooling  EDM  show  superior 
compressive  strength  at  room  and  at  elevated  temperature  when  compared 
to  sample  machined  by  EDM  with  flooding.  This  is  due  to  the  presence  of 
residual  stresses  in  the  specimens  that  were  not  fully  immersed  for  cooling 
during  machining. 

•  The  post  mortem  SEM  images  indicate  cleavage  type  failure  in  specimens 
tested  at  460  °C  and  below.  The  failure  mechanism  changes  to  ductile  at 
higher  temperatures. 
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CHAPTER  5 


Quasi-Static  and  Dynamic  Fracture  Initiation  Toughness  of  Ti/TiB 
Layered  Functionally  Graded  Material  under  Thermo-Mechanical 

Loading 

5.1  Abstract 

Quasi-static  and  dynamic  fracture  initiation  toughness  of  Ti/TiB  layered 
Functionally  Graded  Material  (FGM)  is  investigated  using  a  three  point  bend 
specimen.  The  modihed  split  Hopkinson  pressure  bar  (SHPB)  apparatus  in  con¬ 
junction  with  induction  coil  heating  system  is  used  during  elevated  temperature 
dynamic  loading  experiments.  A  simple  and  accurate  technique  has  been  devel¬ 
oped  to  identify  the  time  corresponding  to  the  load  at  which  the  fracture  initiates. 
A  series  of  experiments  are  conducted  at  different  temperatures  ranging  from  room 
temperature  to  800  °G,  and  the  effect  of  temperature  and  loading  rate  on  the  frac¬ 
ture  initiation  toughness  is  investigated.  The  material  fracture  toughness  is  found 
to  be  sensitive  to  temperature  and  the  fracture  initiation  toughness  increases  as  the 
temperature  increases.  Furthermore,  the  fracture  initiation  toughness  is  strain  rate 
sensitive  and  is  higher  for  dynamic  loading  as  compared  to  quasi-static  loading. 

5.2  Introduction 

An  experimental  investigation  is  conducted  to  research  quasi-static  and  dy¬ 
namic  fracture  initiation  toughness  of  Ti/TiB  FGM  under  thermo-mechanical 
loading.  Ti/TiB  is  a  commercially  available  metal/ceramic  FGM  that  has  been 
proposed  for  many  different  applications.  The  basic  idea  behind  the  design  of 
metal/ceramic  FGMs  is  to  reduce  the  problem  associated  with  the  low  toughness 
of  ceramics.  Hence  the  fracture  toughness  is  the  primary  and  limiting  parameter 
to  design  structures  using  metal/ceramic  FGMs.  Due  to  this  fact,  number  of  re¬ 
searchers  have  studied  theoretically  the  fracture  behavior  of  FGMs.  For  example. 


Delale  and  Erdogan  [1],  Eischen  [2]  and  Jin  and  Noda  [3]  solved  crack  problems 
for  non-homogeneous  materials  under  quasi-static  mechanical  loading.  All  these 
investigations  concluded  that  the  inverse-square  root  singularity  at  the  crack  is 
not  affected  by  non-homogeneity.  Similarly,  for  propagating  cracks  in  FGMs,  the 
theoretical  studies  concluded  that  the  crack  tip  helds  are  similar  to  that  of  homoge¬ 
nous  material  provided  that  the  FGMs  are  continuous  and  piece-wise  continuously 
differentiable  [4-6].  Jin  and  Batra  [7]  further  studied  the  effects  of  quasistatic 
loading  conditions,  specimen  size  and  metal  particle  size  on  the  crack  growth  re¬ 
sistance  curve  (R-curve)  and  the  residual  strength  of  a  ceramic-metal  FGM  using 
crack-bridging  concept. 

On  the  other  hand,  there  are  very  limited  experimental  studies  on  the  fracture 
behavior  of  metal/ceramic  FGMs.  Hill  et  al.  [8]  have  studied  the  fracture  behavior 
of  Ti/TiB  FGM  using  single-edged-notched  bend  specimens  and  reported  that  the 
pre-cracking  method  and  residual  stress  have  a  signihcant  effect  on  the  measured 
fracture  toughness.  The  study  is  limited  to  room  temperature  and  quasi-static 
loading.  The  study  by  Kidane  and  Shukla  [9]  on  the  constitutive  behavior  of 
Ti/TiB  FGM  conhrmed  that  the  material  showed  thermal  softening  at  higher 
temperatures  and  this  resulted  in  a  decrease  in  failure  stress  and  an  increase  in 
failure  strain.  But  the  fracture  behavior  of  Ti/TiB  under  dynamic  loading  and  at 
elevated  temperatures  has  not  been  investigated  yet. 

The  purpose  of  this  experimental  study  is  to  £11  that  gap  and  investigate  the 
effect  of  temperature  and  loading  rate  on  the  fracture  initiation  toughness  of  this 
material.  An  experimental  investigation  on  the  fracture  initiation  toughness  of 
Ti/TiB  under  quasi-static  and  dynamic  loading  is  conducted  at  different  tempera¬ 
tures  ranging  from  room  temperature  to  800  °C .  The  fracture  initiation  toughness 
is  found  to  be  temperature  dependant,  increases  with  raise  in  temperature,  and 


Table  5.1.  Composition  and  elastic  properties  of  Ti/TiB  FGM 


Layer  no. 

Vol.  %Ti 

Vol.  % 

TiB 

Thickness 

(mm) 

E  (GPa) 

V 

1 

100 

0 

0.4064 

106 

0.340 

2 

85 

15 

0.7620 

170 

0.278 

3 

70 

30 

0.7620 

227 

0.238 

4 

55 

45 

0.7620 

262 

- 

5 

40 

60 

0.7620 

289 

- 

6 

25 

75 

0.7620 

303 

0.152 

7 

15 

85 

2.1336 

316 

0.140 

also  loading  rate  dependant,  higher  at  dynamic  loading  compared  with  quasi-static 
loading. 

5.3  Material  and  Specimen  Geometry 

The  material  used  in  this  experimental  study  is  a  Ti/TiB  FGM  layered  plate 
supplied  by  BAE  systems.  It  is  a  7.36  mm  thick  pate  with  seven  layers,  no  clear 
or  distinct  interface  between  them,  ranging  from  pure  Ti  on  one  side  to  85  %  TiB 
on  the  other.  The  details  of  the  fabrication  technique  and  material  properties  can 
be  found  in  the  literature  [8,  9].  The  thickness  and  the  properties  of  each  layer  are 
listed  in  Table  5.1. 

Two  different  batches  of  three  point  bend  specimens  are  machined,  from  the 
FGM  plate  discussed  above,  using  an  electrical  discharge  machining  (EDM)  tech¬ 
nique.  The  hrst  batch  of  specimens  is  machined  in  such  a  way  that  the  grada¬ 
tion  is  in  the  thickness  direction,  i.e.  the  crack  is  perpendicular  to  the  gradation 
(Fig.  5.1a).  The  second  batch  of  specimens  is  machined  in  such  a  way  that  the 
crack  is  along  the  gradation  direction  (Fig.  5.1b).  The  thickness  of  the  specimen 
is  selected  from  the  available  FGM  plate  thicknesses,  where  as  the  other  dimen¬ 
sions  are  determined  based  on  ASTM  Test  Methods  for  Measurement  of  Fracture 
Toughness  (E1820-96).  Due  to  the  brittle  nature  of  this  material  no  fatigue  pre¬ 
crack  was  made.  However,  the  crack  was  made  using  EDM  wire  of  0.0508  mm. 
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a) 


Tiside  w=  10.16 

1 

S=  40.64 

b) 


B=3.18 


Figure  5.1.  Three  point  bend  specimen  (a)  FGM  graded  in  the  thickness  direction 
(b)  FGM  graded  in  the  crack  direction  (all  dimensions  are  in  mm) 
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resulting  in  a  very  small  notch  root  radius. 


5.4  Experimental  Procedure 

5.4.1  Quasi-static  Fracture  Initiation  Toughness 

The  quasi-static  fracture  initiation  toughness  is  investigated  using  a  three- 
point  bending  experiment.  The  three  point  bend  specimen  is  placed  between  the 
loading  bar  and  supporting  block  as  shown  in  Fig.  5.2.  The  experiment  is  con¬ 
ducted  at  different  temperatures  ranging  from  25  °C  to  800  °C.  A  tungsten  carbide 
loading  tip  is  used  to  eliminate  the  deformation  of  the  loading  bar  at  higher  tem¬ 
perature.  The  specimen  is  heated  to  the  desired  temperature  using  an  induction 
heating  system.  The  induction  heater  has  a  precise  controller  that  enables  us  to 
control  the  temperature  to  the  desired  value.  The  experiment  is  conducted  under 
a  hxed  loading  rate  of  Imm/min  at  room  temperature.  To  avoid  the  creep  effect 
during  elevated  temperature,  a  loading  rate  of  5  mm/min  is  used. 


Figure  5.2.  Quasi-static  three  point  bed  experimental  setup. 
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5.4.2  Dynamic  Fracture  Initiation  Toughness 

A  modified  split  Hopkinson  pressure  bar  (SHPB)  apparatus  with  induction 
heating  system  is  used  to  investigate  the  dynamic  fracture  initiation  toughness 
of  Ti/TiB  FGM.  The  modified  SHPB  is  a  well-developed  method  to  study  the 
dynamic  fracture  initiation  toughness  of  engineering  materials.  The  principle  of 
modified  SHPB  is  presented  in  brief  and  the  detailed  theory  can  be  found  in  the 
literature  [10,  11].  The  schematic  of  the  modified  SHPB  apparatus  is  shown  in 
Fig.  5.3.  As  displayed  in  the  figure,  the  apparatus  mainly  consists  of  an  incident 
bar,  striker  bar  and  pressure  gun.  To  reduce  the  impendence  mismatch  between 
the  specimen  and  the  incident  bar,  a  T6061  aluminum  material  is  used  for  the 
incident  bar.  To  generate  a  well-defined  loading  pulse,  the  same  material  that  is 
used  for  the  incident  bar  is  used  for  the  striker  bar.  Two  semiconductor  strain 
gages  are  attached  in  the  middle  of  the  incident  bar  diametrically  opposite  to  one 
another.  The  semiconductor  strain  gages  are  about  75  times  more  sensitive  than 
the  foil  type  strain  gages  and  are  preferred  to  capture  very  small  strain  signals. 


Figure  5.3.  Schematic  representation  of  modified  SHPB  apparatus  with  induction 
heater. 

During  loading,  the  specimen  is  sandwiched  between  the  incident  bar  and  the 
rigid  frame.  For  an  elevated  temperature  experiment,  the  bar  is  first  kept  apart 
and  the  specimen  is  heated  to  the  desired  temperature  (usually  about  20  °C  to  50 
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°C  higher  than  the  test  temperature)  and  later  the  bar  is  brought  manually  into 
contact  with  the  specimen.  The  temperature  of  the  specimen  is  monitored  by  0.127 
mm  chromel-alumel  thermocouple,  which  is  spot  welded  onto  the  specimen.  Once 
the  specimen  is  in  contact  with  the  incident  bar,  the  striker  bar  is  propelled  towards 
the  incident  bar  using  an  air-operated  gun.  The  impact  generates  a  compressive 
stress  wave  in  the  bar  which  propagates  toward  the  bar/specimen  interface.  When 
the  wave  reaches  the  specimen,  some  of  the  wave  is  reflected  back  and  part  of  the 
wave  is  transmitted  into  the  specimen.  The  incident  and  reflected  strain  signals 
are  recorded  using  a  Vishay  2301A  signal-conditioning  ampliher  that  is  connected 
with  an  oscilloscope.  The  load  history  at  the  specimen/bar  interface  is  obtained 
from  the  recorded  strain  data  using  a  one  dimensional  elastic  wave  theory  [10] 
given  by  Eq.  5.1. 

F{t)  =  {e,{t)  +  er{t))EA  (5.1) 


where  F  is  the  force,  and  are  the  incident  and  reflected  strain  pulses,  E  is  the 
Young’s  modulus,  and  A  is  the  cross  sectional  area  of  the  bar. 

In  order  to  avoid  the  transient  effects  all  dynamic  experiments  are  conducted 
at  a  striker  speed  of  1.0  m/sec.  When  the  time  of  fracture  is  sufficiently  long,  the 
dynamic  stress  intensity  factor  can  be  calculated  from  the  input  load  as 


JFLf(—\ 
B\/w’  VM'V 


(5.2) 


where  Kj  is  the  stress  intensity  factor,  B  the  specimen  thickness,  W  the  specimen 
width,  a  the  initial  crack  length  and  f\a/W)  is  a  geometric  factor. 

The  dynamic  fracture  initiation  toughness  {Kid)  corresponds  to  the  stress 
intensity  factor  at  the  time  of  crack  initiation,  i.e.  Kjd  =  Kj  {t initiation)  ■ 

One  of  the  challenges  in  the  dynamic  fracture  initiation  toughness  experiment 
is  to  hnd  a  method  to  accurately  determine  the  crack  initiation  time.  Researchers 
[12-14]  used  a  strain  gage  to  detect  the  fracture-initiation  time  at  high  loading 
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rate.  Usually  the  gage  is  placed  on  the  specimen  near  the  crack  tip,  and  when  the 
fracture  initiates,  the  sudden  decrease  of  strain  as  a  function  of  time  corresponds  to 
the  crack  initiation  point.  However,  the  time  that  the  stress  wave  travels  from  the 
crack  tip  to  the  strain  gage  has  to  be  deducted  to  obtain  the  correct  crack  initiation 
time.  This  is  difficult  in  the  case  of  materials  with  unknown  stress  wave  speed.  In 
the  present  study,  a  simple  and  accurate  technique  is  developed  to  identify  the  time 
corresponding  to  the  load  required  to  initiate  the  crack.  A  conductive  silver  paint 
is  placed  ahead  of  the  crack  tip  on  a  non-conductive  substrate,  in  the  direction  of 
expected  crack  path.  The  paint  is  connected  to  a  strain  gage  which  is  attached  to 
the  loading  bar  at  the  same  position  where  the  two  semiconductor  strain  gages  are 
located.  When  the  fracture  initiates,  the  sudden  jump  in  the  strain  signal  indicates 
the  crack  initiation  time,  and  this  is  used  later  to  determine  the  load  associated 
with  the  fracture  initiation  time. 

A  schematic  of  strain  signals  from  the  semiconductor  strain  gages  and  the 
foil  type  strain  gage  connected  with  the  silver  paint  are  shown  in  Fig.  5.4.  In  the 
hgure,  ti  is  the  duration  of  the  incident  wave,  t2  is  the  duration  of  the  reflected 
wave,  ts  is  the  time  required  by  the  wave  to  travel  from  the  location  of  the  strain 
gages  to  the  bar-specimen  interface  and  to  return  back  to  the  location  of  the  strain 
gages  and  t  is  the  total  time  duration  starting  from  the  wave  passing  the  location 
of  the  strain  gages  until  the  fracture  occurs  (silver  paint  broken).  Hence,  t^/2  is 
the  time  required  by  the  wave  to  travel  from  the  location  of  the  strain  gages  to  the 
bar-specimen  interface.  Therefore,  t  -  ^3/2  is  the  time  required  to  initiate  the 
fracture.  Also  the  location  of  this  time  {ti)  in  the  incident  and  reflected  signal  is 
shown  by  the  dotted. 

Using  the  above  conhguration  a  series  of  experiments  are  conducted  at  dif¬ 
ferent  temperatures.  For  all  experiments  a  pulse  shaper  is  placed  on  the  impact 
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Figure  5.4.  A  schematic  of  strain  signals  from  the  semiconductor  strain  gages  and 
the  foil  type  strain  gage  connected  with  silver  paint. 

face  of  the  incident  bar.  The  pulse  shaper  smoothes  the  relative  sharp  front  of  the 
incident  stress  wave,  thus  allowing  high  strain-rate  experiments  to  be  conducted  at 
near  constant  strain  rates  [9] .  This  also  helps  the  premature  failure  of  the  relatively 
brittle  Ti/TiB  FGM,  especially  during  the  early  part  of  stress  wave  loading  [9]. 

5.5  Results  and  Discussions 

5.5.1  Quasi-static  Fracture  Initiation  Tonghness  of  FGM  Graded  in 
the  Thickness  Direction 

Fig.  5.5  shows,  a  typical  load  history  curve  for  a  quasi-static  three  point  bend 
experiment  of  FGM  graded  in  the  thickness  direction.  From  the  plot,  it  can  be 
seen  that  there  is  a  clear  and  well  dehned  peak  load  at  which  the  crack  initiates. 
The  load  history  shown  above  represents  the  integrated  structure  resistance,  not 
the  load  in  each  separate  layer.  The  fracture  resistance  across  the  thickness  in 
each  layer  may  vary,  as  observed  by  Wadgaonkar  and  Parameswaran  [15].  This 
phenomenon  is  also  observed  from  the  failed  specimen  subjected  to  quasi-static 
loading  shown  in  Fig.  5.6.  As  shown  in  the  hgure,  the  fracture  topography  of  the 
surface  is  different  from  one  end  to  the  other.  The  fracture  surface  near  the  TiB 
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rich  layer  is  smooth  and  the  surface  near  the  Ti  rich  layer  is  rough.  This  indicates 
that,  the  energy  required  to  fracture  the  surface  is  different  for  the  two  extreme 
ends,  and  the  fracture  may  first  start  from  the  TiB  rich  end  and  later  propagate 
towards  Ti  rich  side. 


Figure  5.5.  Typical  load-time  plot  for  FGM  graded  in  the  thickness  direction  under 
quasi-static  loading  (25  °C). 

However,  the  intention  of  the  present  study  is  to  determine  the  fracture  ini¬ 
tiation  toughness  of  the  integrated  structure,  and  not  the  fracture  resistance  of 
each  layer.  Hence,  the  single  load  history  obtained  from  the  integrated  structure 
is  used  to  calculate  the  fracture  initiation  toughness.  By  following  a  similar  proce¬ 
dure,  a  series  of  experiments  are  conducted,  and  the  fracture  initiation  toughness 
at  different  temperatures  is  obtained. 

Fig.  5.7  shows  the  quasi-static  fracture  initiation  toughness  of  FGM  graded  in 
the  thickness  direction  as  a  function  of  temperature.  A  total  of  three  experiments 
are  conducted  at  each  temperature  and  the  bars  indicate  the  range  of  values  ob- 
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Figure  5.6.  Typical  failed  specimen  of  FGM  graded  in  thickness  direction  under 
quasi-static  loading  (25  °C). 

tained.  The  room  temperature  quasi-static  fracture  initiation  toughness  is  about 
4.75  ±  0.07  MPa  and  the  value  increase  to  10.43  ±  0.52  MPa  at  800 
°C.  The  quasi-static  fracture  initiation  toughness  increases  linearly  as  the  test 
temperature  raises.  It  can  be  argued  that  at  elevated  temperature  the  material  is 
thermally  softened  and  this  results  in  an  increase  in  fracture  toughness. 

5.5.2  Quasi-static  Fracture  Initiation  Tonghness  of  FGM  Graded  in 
the  Crack  Direction 

In  this  experiment  the  specimen  is  first  pre-cracked  to  2.54  mm  from  TiB  rich 
side,  which  means  up  to  the  middle  of  75  %  TiB  and  25  %  Ti  layer,  and  the  crack 
is  allowed  to  run  towards  Ti  rich  end  as  shown  in  Fig.  5.8.  The  crack  is  chosen 
to  start  from  the  TiB  side,  as  the  fracture  resistance  of  TiB  is  much  lower  than 
that  of  Ti.  Fig.  5.9  shows  a  typical  load  history  curve  for  a  quasi-static  fracture 
experiment  at  25  °C.  The  figure  displays  a  definite  and  clear  peak  load  followed 
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Temperature  (°C) 

Figure  5.7.  Effect  of  temperature  on  the  fracture  initiation  toughness  of  FGM 
graded  in  the  thickness  direction  under  quasi-static  loading. 

by  yet  another  peak  load  that  occurred  after  the  decline  of  the  Erst.  The  first  peak 
load  corresponds  to  the  load  required  to  initiate  the  crack,  in  75  %  TiB  and  25 
%  Ti  layer,  and  the  second  peak  load  corresponds  to  the  load  required  to  crack 
the  next  layer,  i.e.  60  %  TiB  and  40  %  Ti,  which  is  0.35  mm  away  from  the 
initial  position  of  crack  tip.  Fig.  5.10  shows  a  typical  fracture  surface,  investigated 
at  room  temperature.  From  the  hgure  it  is  clear  where  the  second  peak  load 
occurred.  Furthermore,  it  is  observed  the  fracture  surface  near  to  the  TiB  rich 
layer  is  smooth  and  is  rough  towards  the  Ti  rich  side.  Even  though  the  crack  has 
already  started  at  the  first  peak  load,  still  more  energy  is  required  for  the  crack  to 
propagate  to  the  next  layer. 

Fig.  5.11  shows  a  typical  load  time  curve  for  a  quasi-static  fracture  experiment 
at  elevated  temperatures.  Unlike  the  room  temperature  experiment,  there  is  no 
repeated  peak  load,  the  load  decreased  slowly  after  it  reached  the  critical  value. 
From  the  load  history  it  can  be  clearly  seen  that  the  temperature  reduces  the 


99 


Figure  5.8.  A  schematic  representation  of  crack  on  FGM  specimen  graded  in  the 
crack  direction. 
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Figure  5.9.  Typical  load-time  plot  for  FGM  graded  in  the  crack  direction  under 
quasi-static  and  room  temperature  loading  (25  °G). 
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Figure  5.10.  Typical  failed  specimen  of  FGM  graded  in  the  crack  direction  under 
quasi-static  loading  (25  "C). 

brittleness  of  the  material.  For  all  the  cases  considered,  the  hrst  peak  load  is  used 
to  calculate  the  fracture  initiation  toughness. 

Fig.  5.12  shows  the  fracture  initiation  toughness  as  a  function  of  temperature 
for  FGM  graded  in  the  crack  direction.  A  total  of  three  experiments  are  conducted 
at  each  temperature  and  the  bars  indicate  the  range  of  values  obtained.  The  quasi¬ 
static  fracture  initiation  toughness  at  room  temperature  (25  °C)  is  about  3.75  ± 
0.05  MPa  and  increases  to  about  9.55  ±  0.44  MPa  at  800  °C.  The 
trend  on  the  quasi-static  fracture  initiation  toughness  of  FGM  graded  in  the  crack 
direction  is  same  as  that  of  FGM  graded  in  the  thickness  direction.  With  both,  the 
fracture  initiation  toughness  increases  consistently  with  an  increase  in  temperature. 


5.5.3  Dynamic  Fracture  Initiation  Toughness  of  FGM  Graded  in  the 
Thickness  Direction 

A  typical  strain  record  from  one  of  the  dynamic  experiments  of  FGM  graded 
in  the  thickness  direction  is  shown  in  Fig.  5.13.  The  sudden  jump  in  strain  signal 
from  the  strain  gage,  which  was  attached  to  the  conductive  silver  paint,  indicates 
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Figure  5.11.  Typical  load  -  time  plot  for  FGM  graded  in  the  crack  direction  under 
quasi-static  elevated  temperature  loading. 


Figure  5.12.  Effect  of  temperature  on  the  fracture  initiation  toughness  of  FGM 
graded  in  the  crack  direction  under  quasi-static  loading. 
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the  crack  initiation  time.  As  shown  in  the  hgure,  the  crack  initiates  140  /is  after  the 
stress  wave  enters  the  specimen  and  this  allows  the  wave  to  reverberate  sufficiently 
to  create  an  equilibrium  loading  condition  before  the  crack  initiates.  Fig.  5.14 


Figure  5.13.  Typical  incident  and  reflected  strain  plot  of  FGM  graded  in  the 
thickness  direction  subject  to  room  temperature  dynamic  loading. 

shows  the  dynamic  fracture  initiation  toughness  of  FGM  graded  in  the  thickness 
direction  as  a  function  of  temperature.  A  total  of  three  experiments  are  conducted 
at  each  temperature  and  the  bars  indicate  the  range  of  values  obtained.  At  room 
temperature  the  dynamic  fracture  initiation  toughness  is  about  5.2  ±  0.20  MPa 
about  11.96  ±  0.77  MPa  at  800  °C .  It  is  also  observed  that  the 
fracture  initiation  toughness  increases  linearly  with  increase  in  temperature.  This 
can  be  attributed  to  the  fact  that,  at  higher  temperature  the  material  is  softened 
and  results  in  increased  fracture  toughness.  This  phenomenon  agrees  with  the 
one  observed  by  the  authors  in  their  previous  study  during  the  investigation  of  the 
dynamic  constitutive  behavior  of  this  material  [9].  Fig.  5.15  shows  the  SEM  images 
of  typical  fracture  surface  of  FGM  graded  in  the  thickness  direction  subjected  to 
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dynamic  loading  at  different  temperatnres.  The  images  are  taken  near  the  crack 
tip  at  the  Ti  rich  end.  It  can  be  seen  that  the  fractnre  snrface  becomes  more 
dnctile  as  the  temperatnre  increases. 


Temperature  (°C) 

Fignre  5.14.  Effect  of  temperatnre  on  the  fractnre  initiation  tonghness  of  FGM 
graded  in  the  thickness  direction  nnder  dynamic  loading. 

5.5.4  Dynamic  Fracture  Initiation  Toughness  of  FGM  Graded  in  the 
Crack  Direction 

A  typical  strain  signal  obtained  from  the  dynamic  fractnre  experiment  for 
FGM  graded  in  the  crack  direction  is  shown  in  Fig.  5.16.  As  discnssed  in  the 
previous  section,  the  conductive  silver  paint  is  used  to  identify  the  load  which 
corresponds  to  the  fracture  initiation  time.  The  fracture  initiates  after  60  fis 
and  this  allowed  the  wave  to  reverberate  about  20  times  in  the  specimen  before 
crack  initiates.  This  is  sufficient  to  create  an  equilibrium  loading  condition  before 
the  crack  initiates.  High  speed  digital  imaging  is  also  used  to  study  the  contact 
history  of  the  bar  and  specimen  and  the  fracture  process.  Fig.  5.17  shows  typical 
high  speed  images  taken  during  a  dynamic  fracture  experiment  of  FGM  graded  in 
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c)  T=  600  °C 


d)  T=  800  °C 


Figure  5.15.  SEM  images  of  typical  fracture  surface  of  FGM  graded  in  the  thickness 
direction  subjected  to  dynamic  loading  at  different  temperatures. 
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the  crack  direction  at  a  temperature  of  800  °C.  It  can  be  clearly  seen  that  the  bar 
is  always  in  contact  with  the  specimen  during  the  failure  process.  Hence,  the  use 
of  load  corresponding  to  the  crack  initiation  time  is  justihed  in  the  calculation  of 
calculate  the  dynamic  fracture  initiation  toughness. 


Figure  5.16.  Typical  strain  signal  obtained  from  the  dynamic  fracture  experiments 
of  FGM  graded  in  the  crack  direction  subjected  to  dynamic  loading. 

Fig.  5.18  shows  the  dynamic  fracture  initiation  toughness  as  a  function  of 
temperature  for  FGM  along  the  crack  direction.  A  total  of  three  experiments  are 
conducted  at  each  temperature  and  the  bars  indicate  the  range  of  values  obtained. 
A  similar  trend  as  observed  in  the  previous  batch  of  specimen  is  seen  here,  the 
dynamic  fracture  initiation  toughness  increases  with  increase  in  temperature.  At 
room  temperature  the  dynamic  fracture  initiation  toughness  is  about  5.0  ±  0.07 
MPa  and  increases  to  10.55  ±  0.85  MPa  at  800  °C.  SEM  images  of  the 
micro  structure  at  different  temperatures  at  the  crack  tip  (75  %  TiB,  25  %  Ti)  are 
shown  in  Fig.  5.19.  Some  evidence  of  stretching  associated  with  material  softening 
appears  at  high  temperatures. 
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Figure  5.17.  Typical  high  speed  digital  images  of  FGM  graded  in  the  crack  direc¬ 
tion  subjected  to  dynamic  loading  at  800  °C. 


Figure  5.18.  Effect  of  temperature  on  the  fracture  initiation  toughness  of  FGM 
graded  in  the  crack  direction  under  dynamic  loading. 
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Figure  5.19.  SEM  images  of  typical  fracture  surface  of  FGM  graded  in  the  crack 
direction  subjected  to  dynamic  loading  at  different  temperatures. 
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5.5.5  Effect  of  Loading  Rate  on  the  Practnre  Initiation  Tonghness  of 
FGM 

The  quasi-static  and  dynamic  fracture  initiation  toughness  of  FGM,  through 
the  thickness  direction  and  along  the  crack  direction,  as  a  function  of  tempera¬ 
ture  are  shown  in  Fig.  5.20  and  Fig.  5.21.  In  the  hgure,  it  is  clear  the  fracture 
initiation  toughness  is  load  rate  dependant.  The  fracture  initiation  toughness  at 
all  temperatures  is  higher  in  dynamic  loading  (2.5  MN/sec)  as  compared  with  the 
corresponding  quasi-static  loading  (100  N/Sec). 


Temperature  (°C) 

Figure  5.20.  Effect  of  temperature  on  the  fracture  initiation  toughness  of  FGM 
graded  in  the  thickness  direction  under  quasi-static  and  dynamic  loading. 


5.6  Summary 

The  quasi-static  and  dynamic  fracture  initiation  toughness  of  Ti/TiB  FGM  is 
investigated  under  thermo-mechanical  loading.  A  simple  and  accurate  technique 
has  been  developed  to  identify  the  time  corresponding  to  the  load  at  which  the 
crack  initiates.  The  Endings  are  summarized  below: 

•  The  fracture  initiation  toughness  is  temperature  dependent,  and  it  increases 
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Temperature  (°C) 

Figure  5.21.  Effect  of  temperature  on  the  fracture  initiation  toughness  of  FGM 
graded  in  the  crack  direction  under  quasi-static  and  dynamic  loading. 

as  the  temperature  increases.  For  FGM  graded  in  the  thickness  direction, 
as  the  temperature  increase  from  room  temperature  to  800  "G,  the  quasi¬ 
static  fracture  initiation  toughness  increased  by  160  %  and  the  dynamic 
fracture  initiation  toughness  increased  by  130  %.  For  FGM  graded  in  the 
crack  direction,  as  the  temperature  increased  from  room  temperature  to  800 
°G,  the  quasi-static  fracture  initiation  toughness  increased  by  154  %  and  the 
dynamic  fracture  initiation  toughness  increased  by  110  %. 

•  The  fracture  initiation  toughness  is  rate  dependent;  it  is  higher  at  dynamic 
loading  as  compared  with  quasi-static  loading.  For  FGM  graded  in  the  thick¬ 
ness  direction,  the  room  temperature  fracture  initiation  toughness  increased 
by  30  %  and  the  fracture  initiation  toughness  at  800  °C  increased  by  15  %  as 
the  loading  rate  increased  from  100  N/sec  to  2.5  MN/sec.  For  FGM  graded 
in  the  crack  direction,  the  room  temperature  fracture  initiation  toughness 
increased  by  33  %  and  the  fracture  initiation  toughness  at  800  °C  increased 
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by  10  %  for  the  same  increasing  in  the  loading  rate. 
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CHAPTER  6 


Dynamic  Crack  Propagation  in  Transparent  Functionally  Graded 

Material 

6.1  Abstract 

A  model  transparent  graded  materials  is  nsed  to  investigate  the  steady  state 
and  transient  crack  propagation  in  a  fnnctionally  graded  material.  A  transparent 
graded  material  was  fabricated  using  two  different  materials  systems  of  polymer 
resin,  which  were  cast  using  a  special  mold;  to  give  continuously  graded  properties 
along  the  crack  propagation  direction.  High-speed  digital  photography  combined 
with  photoelasticity  technique  is  used  to  record  the  full-held  stress  data  around  the 
propagating  crack.  By  analyzing  the  photoelastic  fringe  patterns  the  propagation 
velocity  of  the  crack  tip  are  obtained. 

6.2  Introduction 

In  this  study,  dynamic  fracture  in  transparent  FGM  has  been  investigated  for 
different  crack  conhguration  and  specimen  geometries  to  completely  characterize 
the  behavior  of  a  moving  crack  in  FGMs.  The  specimen  geometry  used  in  this 
study  was  Single  Edge  Notch  Tension  (SENT).  Dynamic  photoelastic  studies  along 
with  high-speed  photography  were  applied  to  obtain  full-held  stress  data  for  a 
propagating  crack  in  above-mentioned  FGMs.  The  real  time  images  were  analyzed 
to  calculate  the  crack  length,  crack  velocity  and  dynamic  stress  intensity  factor. 
FGMs  are  materials,  which  have  continuously  varying  mechanical  properties  in  a 
specihc  direction.  Even  though  FGM  is  a  preferable  material  for  diherent  specihc 
applications,  few  studies  are  available  on  the  dynamic  fracture  behavior  of  these 
materials.  This  calls  for  further  investigation  of  fracture  behavior  of  FGM  under 
diherent  loading  conditions.  Photoelasticity  technique  with  high  speed  digital 
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imaging  system  is  one  of  the  most  widely  used  technique  to  study  the  fracture 
behavior  of  materials  and  was  used  in  this  study. 

Analytical  work  on  functionally  graded  materials  goes  back  as  early  as  the 
late  1960s  when  soil  was  modeled  as  a  nonhomogeneous  material  by  Gibson  [1], 
More  recently,  for  propagating  cracks  in  FGMs,  Parameswaran  and  Shukla  [2]  and 
Ghalivendra  et  ah  [3],  developed  the  structure  of  the  first  stress  invariant  and  the 
out  of  plane  displacement.  In  their  study  they  brought  out  the  effects  of  non- 
homogeniety  through  an  asymptotic  analysis.  Lee  [4]  developed  nonhomogeneity 
specihc  terms  for  individual  stress  and  displacement  components  using  displace¬ 
ment  potentials.  A  review  paper  by  Shukla  et  ah  [5]  presents  a  comprehensive 
summary  of  dynamic  fracture  studies  in  FGM. 

In  contrast  there  are  only  few  experimental  studies  on  the  dynamic  fracture  of 
FGM.  Parameswaran  and  Shukla  [6]  used  photoelasticity  to  investigate  dynamic 
fracture  in  FGMs  with  discrete  property  variation.  Li  et  ah  [7]  have  investigated 
the  quasi  static  fracture  properties  of  FGM  crated  by  selective  ultraviolet  (UV) 
irradiation  of  a  polymer.  Rousseau  and  Tippur  [8]  studied  the  crack-tip  defor¬ 
mation  in  FGMs  subjected  to  low  velocity  impact.  Recently,  Jain  and  Shukla  [9] 
investigated  the  transient  crack  growth  behavior  in  FGM  fabricated  as  a  partic¬ 
ulate  composite  with  continuously  varying  particle  volume  fraction  along  a  single 
dimension. 

Apart  from  a  few  studies  reported  above  there  is  no  sufficient  experimental 
work,  explaining  the  stress  helds  near  the  crack  tip.  In  this  experimental  study, 
dynamic  photoelasticity  and  high-speed  digital  imaging  were  used  to  study  the 
dynamic  fracture  behavior  of  FGM.  The  crack  tip  stress  helds  were  recorded  at 
different  time  intervals.  The  crack  velocity,  acceleration  and  the  dynamic  stress 
intensity  factor  have  been  calculated. 
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6.3  Material  and  Specimen  Preparation 

A  transparent  graded  material  was  fabricated  and  the  mechanical  properties 
was  determined,  at  Sangju  National  University,  South  Korea,  using  two  different 
materials  systems  of  polymer  resin.  This  was  accomplished  by  pouring  the  two 
resins  in  a  specially  prepared  mold.  Resin  1  was  poured  first  in  an  inclined  mold 
and  was  let  to  cure  till  gelation  occurred.  Resin  2  was  then  poured  at  top  of 
Resin  1  and  let  to  cure.  These  two  polymers  had  different  Young’s  modulus, 
different  fringe  constants  but  the  same  density.  The  resulting  composite  created 
an  inclined  interface  as  shown  in  Fig.  6.1  and  produced  graded  properties.  The 
Young’s  modulus  and  Poisson’s  ratio  of  these  two  resins  were  obtained  using  a 
cross  strain  gage  in  tensile  load.  The  stress  fringe  constant  was  determined  using 
a  compressive  diamond  test.  The  materials  and  photoelasticity  properties  of  these 
resins  are  shown  in  Table.  6.1. 
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Figure  6.1.  In-house  fabricated  typical  SENT  specimen:  a)  FGM  with  increasing 
Young’s  modulus  in  the  direction  of  crack  propagation,  b)  FGM  decreasing  in 
young’s  modulus  in  the  direction  of  crack  propagation  (Dimensions  W=127mm, 
H/W=2,  h/W=1.45) 
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Table  6.1.  Material  properties  of  polymer  resins 


Polymer  type 

Young’s 

Poisson’s 

Stress 

Density 

modulus 

(GPa) 

ratio 

fringe 

constant 

(kN/m) 

(P) 

Resin  1  (hard) 

2.45 

0.38 

17.30 

1200 

Resin  2  (soft) 

0.46 

0.50 

5.10 

1200 

6.4  Experimental  Procedure 

6.4.1  Dynamic  Fracture  at  Constant  Velocity 

A  schematic  of  the  photoelastic  conhguration  used  for  studying  dynamic  crack 
propagation  is  shown  in  Fig.  6.2.  High-speed  digital  imaging  was  employed  along 
with  dynamic  photoelasticity  to  obtain  real  time,  full-field  quantihcation  of  the 
dynamic  failure  process.  Two  circular  polarizers  were  placed  on  either  side  of  the 
specimen  to  form  a  light-held  circular  polariscope.  Power  Light  2500DR  xenon 
hash  lamps  (Photogenic  Inc.)  were  used  as  light  sources  to  illuminate  the  speci¬ 
men  and  a  monochromatic  hlter  was  placed  just  before  the  camera  to  ensure  that 
the  imaged  isochromatic  fringe  patterns  correspond  to  a  single  wavelength  of  light. 
The  isochromatic  images  were  captured  at  preset  time  interval  using  an  Imacon 
200  ultra-high-speed  digital  camera.  This  CCD  based  camera  provides  16  indepen¬ 
dently  programmable  digital  images  of  dynamic  events  up  to  a  maximum  framing 
rate  of  200  XIO®  frames/s.  The  camera  has  a  facility  to  choose  diherent  range  of 
interframe  times.  Interframe  times  were  chosen  so  as  to  ensure  that  the  dynamic 
fracture  event  was  captured  within  the  16  frames. 

Single  Edge  Notched  Tension  (SENT)  specimens  were  machined  from  FGM 
material  prepared  in  the  procedure  explained  above.  A  series  of  experiments  was 
conducted  with  two  different  SENT  specimen  geometries.  The  conhgurations  con¬ 
sidered  in  this  study  were  (a)  crack  propagating  towards  increasing  fracture  tough¬ 
ness  direction,  Fig.  6.1a  and  (b)  crack  propagating  towards  decreasing  fracture 
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toughness  directions,  Fig.  6.1b.  During  the  experiment  the  specimens  were  stat¬ 
ically  loaded  to  a  predetermined  value  of  initiation  stress  intensity  factor  using 
INSTRON  5585  apparatus.  The  crack  was  initiated  by  tapping  a  sharp  razor 
blade  on  the  specimen  notch.  Upon  initiation  of  crack,  a  conductive  silver  paint 
lined  near  the  crack  tip  in  the  direction  of  crack  growth  and  connected  with  an 
electronic  circuit  triggered  a  high-speed  camera  IMACON  200.  The  camera  was 
set  to  take  photographs  at  a  rate  of  50,000  frames  per  second  and  an  exposure 
time  of  200  ns.  The  camera  recorded  16  frames  of  the  crack  propagation  event  (0- 
400ps). 

6.4.2  Transient  Crack  Propagation 

Dynamic  fracture  in  the  FGM  material  explained  above  for  transient  crack  tip 
loading  was  also  investigated.  The  modified  SENT  specimen  shown  in  Fig.  6.3  was 
used.  The  same  experimental  procedure  as  explained  above  was  followed.  However, 
in  this  case  the  camera  was  set  at  a  rate  of  140,000  frames  per  sec.  During  loading, 
first  the  crack  propagates  towards  the  circular  hole  and  arrests  until  the  load  reach 
its  critical  value  to  propagates  into  the  second  crack.  Once  again  the  load  reach 
the  critical  value  it  suddenly  accelerates  resulting  in  a  transient  crack  propagation. 


6.5  Results  and  Discussion 
6.5.1  Crack  tip  velocity 

Typical  isochromatic  fringes  obtained  from  the  dynamic  fracture  experiments 
with  constantly  moving  crack  tip  loading  are  shown  in  Fig.  6.4  and  6.5.  From 
isochromatic  fringes  shown  above,  it  was  observed  that,  there  is  a  closure  force  act¬ 
ing  behind  the  crack  tip.  This  can  be  clearly  seen  in  the  case  of  crack  propagating 
towards  decreasing  fracture  toughness  (increasing  Young’s  modulus)  direction  as 
shown  in  Fig.  6.5.  This  is  expected  because  at  the  left  end  of  the  specimen,  the 
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Figure  6.3.  Modified  SENT  specimen  for  transient  crack  propagation  in  FGM 
(Dimensions  W=127mm,  H/W=2,  h/W=1.45,  ai=12mm,  02=  5  mm,  D=5  mm  1 
=20mm). 
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beginning  of  the  crack,  the  softer  part  of  the  material  is  thick  and  has  a  tendency 
to  stay  together  as  the  crack  propagates  through  the  harder  layer.  The  crack  tip 
distance  with  the  corresponding  time  is  determined  from  the  images.  Using  these 
images  the  crack  length  as  a  function  of  time  is  plotted  in  Fig.  6.6.  The  data 
is  later  used  to  determine  the  crack  tip  velocity  as  a  function  of  time  as  shown 
in  Fig.  From  these  figures  it  can  be  seen  that  for  both  the  cases,  for  increasing 
and  decreasing  stiffness  in  the  crack  propagation  direction,  the  crack  velocity  is 
constant  and  about  350  m/s. 


Frame  1 ,  a=0. 048m  Frame  2,  a=0.054m  Frame  3,  a=0.062  m 


Frame  4,  a=0.068m  Frame  5 ,  a=0.075m  Frame  6,  a=0.082  m 


Figure  6.4.  The  dynamic  isochromatic  fringes  at  the  crack  tip  for  FGM  with 
decreasing  Young’s  modulus  along  the  crack  propagating  direction.  Time  interval 
between  each  frame  is  20  s 

Typical  isochromatic  images  obtained  from  the  dynamic  fracture  of  transient 
crack  propagation  loading  experiments  are  shown  in  Fig.  6.7.  Using  these  images, 
the  crack  length  and  velocity  as  a  function  of  time  is  calculated  and  plotted  as 
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Frame  1,  a=0.047m  Frame  2,  a=0.054m  Frame  3,  a=0.061  m 


Frame  4,  a=0. 068m  Frame  5,  a=0. 075m  Frame  6,  a=0.083  m 


Figure  6.5.  The  dynamic  isochromatic  fringes  at  the  crack  tip  for  FGM  with 
increasing  Young’s  modulus  along  the  crack  propagating  direction.  Time  interval 
between  each  frame  is  20  s 
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Figure  6.6.  Crack  length  and  velocity  as  a  function  of  time  for  linearly  varying 
FGM  under  dynamic  loading 
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shown  in  Fig.  6.8.  It  is  observed  that  the  crack  tip  velocity  increases  at  high 
acceleration  until  it  reaches  a  maximum  of  550  m/s  and  470  m/s  for  increasing 
and  decreasing  stiffness  respectively.  Then  the  velocity  oscillates  with  an  average 
value  of  300  m/s  for  the  rest  of  crack  propagation  event. 


Frame  1,  a=0. 0  28m  Frame  2,  a=0. 028m  Frame  3,  a=0.028  m 


Frame  4,  a=0. 033m  Frame  5,  a=0. 034m  Frame  6,  a=0.036  m 


Figure  6.7.  The  dynamic  isochromatic  fringes  for  transient  crack  propagation 
on  FGM  with  increasing  Young’s  modulus  along  the  crack  propagating  direction. 
Time  interval  between  each  frame  is  7  s 

6.5.2  Dynamic  Stress  Intensity  Factor 

The  isochromatic  fringes  obtained  during  the  dynamic  experiment  were  used 
to  extract  the  dynamic  stress  intensity  factors.  It  can  be  seen  that  the  isochro¬ 
matics  are  symmetric  about  the  crack  plane,  which  indicates  the  opening  mode 
(mode  I)  dominant  loading  at  the  crack-tip.  In  the  present  analysis,  the  dynamic 
stress  helds  for  FGM  with  linearly  varying  properties  developed  by  Lee  [4]  is  used. 
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Figure  6.8.  Crack  length  and  velocity  as  a  function  of  time  for  linearly  varying 
FGM  under  transient  dynamic  loading 

Isochromatics  from  each  of  the  frames  were  analyzed  using  the  aforementioned 
photoelastic  procedure  to  obtain  the  stress  intensity  factor  history.  A  detailed 
analysis  of  the  isochromatic  fringes  associated  with  the  dynamic  fracture  behavior 
of  the  crack,  which  contains  crack  initiation  and  propagation  can  be  found  [3,  4,  6]. 

However  the  available  stress  held  does  not  support  the  gradation  across  the 
thickness  direction.  In  the  present  experiment,  it  is  observed  that  the  gradation 
through  the  thickness  direction  has  effect  on  the  fracture  behavior  especially  for  the 
case  of  decreasing  Young’s  modulus  in  the  crack  propagation  direction.  Hence,  fur¬ 
ther  investigation  is  required  to  accurately  determine  the  dynamic  stress  intensity 
factor  . 


6.6  Summary 

Simple  procedure  for  preparing  a  transparent  model  FGMs  with  continuous 
gradation  of  properties  was  developed  using  two  different  polymer  resins.  The 
physical,  elastic  and  fracture  properties  of  the  prepared  FGMs  were  characterized. 
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Dynamic  photoelasticity  along  with  high  speed  photography  was  used  to  investi¬ 
gate  the  behavior  of  moving  crack  in  FGMs.  The  propagation  crack  tip  velocity 

and  displacement  are  obtained. 
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CHAPTER  7 


Conclusions  and  Recommendations 

7.1  Conclusions 

During  this  project  a  fundamental  study  was  conducted  to  evaluate  dynamic 
fracture  and  failure  in  functionally  graded  materials  (FGMs)  at  room  and  at  high 
temperatures.  Room  temperature  studies  included  experimental  and  analytical 
evaluation  of  dynamic  fracture  in  model  materials.  Analytical  studies  are  con¬ 
ducted  to  develop  thermo-mechanical  stress,  strain  and  displacement  helds  around 
stationary  cracks  and  around  cracks  propagating  under  steady  state  conditions. 
An  experimental  study  is  also  investigated  on  the  dynamic  constitutive  and  frac¬ 
ture  initiation  toughness  of  the  commercially  available  functionally  graded  material 
(Ti/TiB).  The  hndings  from  the  presented  study  are  summarized  bellow. 

•  The  stress-helds  near  the  crack  tip  for  mixed  mode  thermo-mechanical  load¬ 
ing  in  homogeneous  and  graded  material  are  developed  using  displacement 
potentials  in  conjugation  with  an  asymptotic  approach.  Using  the  devel¬ 
oped  equations,  angular  variation  of  maximum  shear  stress,  circumferential 
stress  and  the  largest  principal  stress  are  plotted  as  a  function  of  tempera¬ 
ture  around  the  crack-tip.  Using  both  minimum  strain  energy  density  cri¬ 
terion  and  maximum  circumferential  stress  criterion,  the  crack  instability 
direction  for  various  crack-tip  speeds  and  non-homogeneous  coefficients  is 
also  determined.  It  observed  that  the  crack  extension  angle  for  mixed  mode 
thermo-mechanical  loading  depends  on  the  crack  tip  velocity,  temperature 
and  non-homogeny  parameter. 

•  The  dynamic  constitutive  behavior  of  Ti/TiB  FGM  under  thermo- mechanical 
loading  is  investigated.  The  FGM  showed  thermal  softening  at  higher  tem- 
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perature  with  a  decrease  in  failure  stress  and  an  increase  in  failure  strain. 
Finally  it  is  observed  that,  the  failure  mechanism  changes  to  ductile  at  higher 
temperatures. 

•  The  quasi-static  and  dynamic  fracture  initiation  toughness  of  Ti/TiB  FGM 
is  investigated  under  thermo-mechanical  loading.  A  simple  and  accurate 
technique  has  been  developed  to  identify  the  time  corresponding  to  the  load 
at  which  the  crack  initiates.  Further  it  is  observed  that  the  fracture  initiation 
toughness  is  temperature  and  rate  dependent. 

•  A  model  transparent  graded  material  is  used  to  investigate  the  steady  state 
and  transient  crack  propagation  in  a  functionally  graded  material.  It  ob¬ 
served  that  that  for  both  the  cases,  for  increasing  and  decreasing  stiffness  in 
the  crack  propagation  direction,  the  crack  velocity  is  constant.  On  the  other 
hand;  for  transient  crack  propagation  experiment,  it  is  observed  that  the 
crack  tip  velocity  increases  at  high  acceleration  until  it  reaches  a  maximum 
value  and  then  remain  constant  for  the  rest  of  crack  propagation  event. 

7.2  Recommendations 

The  current  research  represents  a  fundamental  and  first  step  in  understanding 
the  dynamic  fracture  and  failure  of  functionally  graded  materials  under  thermo¬ 
mechanical  loading.  In  the  future  work,  investigation  should  be  conducted  to 
exploit  the  properties  of  graded  multi-functional  materials  under  different  loading 
conditions  and  to  predict  their  performance  and  structural  integrity  under  extreme 
mechanical  and  thermal  environments.  The  high  temperature  facilities  and  the 
dynamic  stress  helds  that  has  been  developed  as  a  part  of  the  present  work  can 
be  utilized.  The  specihc  deliverables  of  the  proposed  project  are  summarized  as 
follows: 
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•  Analytical  development  of  transient  stress,  strain  and  displacement  fields 
for  propagating  cracks  in  FGM  subjected  to  combined  and  rapidly  changing 
thermo-mechanical  loading  under  a  variety  of  boundary  conditions. 

•  Dynamic  response  of  FGM  when  subjected  to  controlled  shock  loading  at 
various  temperatures  utilizing  3D  digital  image  correlation  techniques  and 
high  speed  imaging  under  different  boundary  conditions. 

•  Development  of  constitutive  models  for  various  spatially  tailored  materials  at 
high  temperatures  and  high  strain  rates  of  loading  with  and  without  conhning 
stresses. 

•  Low  and  high  speed  impact  response  of  FGM  at  various  temperatures  and 
boundary  conditions. 

•  Finally,  development  of  physics-based  failure  models  for  the  combined  en¬ 
vironment  response  of  spatially  tailored  materials  and  structures  that  show 
good  correlation  between  simple  laboratory  experiments  and  structural  ap¬ 
plications. 
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APPENDIX  A 


Validation  Plots 


Various  stress  fields  obtained  from  the  analytical  study  presented  in  chapter 
2,  if  mode-I  and  mode-II  helds  are  separated  and  plotted,  are  shown  below.  These 
plots  perfectly  matches  with  the  isothermal  solution  given  by  Freund. 


Figure  A.l.  Variation  of  the  circumferential  tensile  stress  with  angle  9  around  the 
crack  edge  for  several  value  of  normalized  crack  speed.  Crack  growth  is  in  mode 
I.  (r  =  0.002m,  Kid  =  IMPa 
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Figure  A. 2.  Variation  of  maximum  principal  stress  with  angle  6  around  the  crack 
edge  for  several  valne  of  normalized  crack  speed.  Crack  growth  is  in  mode  I.  (r  = 
0.002m,  Kir,  =  IMPa 
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Figure  A. 3.  Variation  of  the  maximum  shear  stress  with  angle  9  around  the  crack 
edge  for  several  value  of  normalized  crack  speed.  Crack  growth  is  in  mode  II.  (r  = 
0.002m,  Kjjd  =  IMPa 
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Figure  A. 4.  Variation  of  the  circumferential  tensile  stress  with  angle  6  around  the 
crack  edge  for  several  value  of  normalized  crack  speed.  Crack  growth  is  in  mode 
II.  (r  =  0.002m,  Kud  =  IMPa 
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Figure  A. 5.  Variation  of  maximum  principal  stress  with  angle  6  around  the  crack 
edge  for  several  value  of  normalized  crack  speed.  Crack  growth  is  in  mode  II.  (r  = 
0.002m,  Kjin  =  IMPa 


132 


APPENDIX  B 


Matlab  Code 


0/  0/  0/  0/ 0/  0, 


0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


'/o'/o  Program  to  calculate  and  plot  Thermo-Mechanical  Stress  and 
yX  Displacement  Fields  for  Propagating  Crack  Tip  in 
yXo  Functionally  Graded  Materials 

'/o'/o  Program  used  to  plot  figures  given  in  Chapter  2  and  3 


y  0  /  0  /  0  /  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


clear  all; 
close  all; 
clc ; 

'/ogeneral  constants 
mue=le09 ; 
rho=1200 ; 
nue=0.3; 

k=2*nue/ (l-2*nue) ; 

ac=8 . 5e-6; 

q0=250 ; 

ql=0; 

q2=0; 


'/oThe  constants  in  the  stress  functions 
KlD=le6;yoinput(’ Enter  the  value  of  KIC); 

K2D=0.2*K1D; 

Kef=sqrt(KlD~2+K2D''2)  ; 

'/oNonhomozeneous  constants 
beta=0 ; 

zeta=0;yoinput( ’Enter  the  value  of  nonhomogeneity  factor:’); 

TD=0.003;7o  input(’the  domain  for  fringe  generation=  ’); 

N0P=51;  ZinputC’ Number  of  points  =’); 

ds=2*TD/(N0P-l) ; 

xc=[-TD:ds:TD]  ; 

yc=[-TD:ds:TD]  ; 

count_i=0; 

idata=l ; 
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mue_c=exp (TD*zeta) *mue ; 
rho_c=exp (TD*zeta) *rho ; 

cs=sqrt (mue_c/rho_c) ; 
cl=cs*sqrt (2+k) ; 

sr=0.5;  "/oratio  of  crack  speed  to  shear  wave  speed  at  crack  tip 
c=sr*cs ; 


al=sqrt (1- (c/cl) "2) ; 
as=sqrt (1- (c/cs) "2) ; 


count_i=0; 

A0=(  (4*  (l+(as"2) ) )/( (4*al*as)-(  (l+(as''2) )  ~2) ) )  *  (KID/ (3*mue*sqrt  (2*pi) ) )  ; 
B0=( (-2*al) / (1+ (as~2) ) ) *A0 ; 

C0=(8*as)/ (3*(4*as*al-(l+as"2) "2) ) * (K2D/ (mue*sqrt (2*pi) ) ) ; 
D0=((l+as~2)/(2*as))*C0; 

H=k*(l-al~2)+2; 

G=k*(l-as~2)+2; 

L=k*(l+al~2)+2; 

M=k*(l+as~2)+2; 

N=k-(al''2)*(k+2) ; 

0=k-(as''2)*(k+2) ; 

P=k+(al''2)*(k+2) ; 

Q=k+(as''2)*(k+2) ; 

Pl=(zeta*k*(al~2) )/  (al~2-as''2)  ; 

P2=(zeta*as/ ( (k+2) * (al~2-as''2) ) ) ; 

Ql=( (3*k+2) / (k+2) ) * (ac/ ( (al~2) -1) ) ; 

Ll=(zeta/(al~2) ) ; 

L2=zeta/as; 

'/oExpresseion  for  Stress 

for  j=l:l:N0P 
r=0.002; 

t(j)=-pi+(j-l)*2*pi/(N0P-l) ; 
y.  t(j)  =  (j-l)*2*pi/((N0P-l))  ; 
x(j)=r*cos(t(j)) ; 
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y(j)=r*sin(t(j)) ; 


yoql=-(l+as-2)/(l-as~2)*r~(-l/2)*qO; 

rl(j)=sqrt((x(j)~2)+((al~2)*(y(j)~2))); 

rs(j)=sqrt((x(j)~2)+((as~2)*(y(j)~2))); 


tl(j)=atan2((al*y(j)) ,x(j)) ; 
ts(j)=atan2((as*y(j)) ,x(j)) ; 

°/ot(j)=atan2(y(j) ,  x(i)); 

x_bar ( j )  =  (x(j ) )  ;°/o*cos(zeta)-(y(j)*sin(zeta)) ; 

sll_l(j)  =  (3*H/4)*A0*(rl(j)''(-l/2))*cos(tl(j)/2)  +  (-3*Ll*H/16) 
*A0* (rl ( j ) ~ (1/2) ) *cos (3*tl ( j ) /2) + (-3*Ll/4) 
*A0*L*(rl(j)-(l/2))*cos(tl(j)/2) ; 
sll_2(j)=-(3*H/4)*C0*(rl(j)~(-l/2))*sin(tl(j)/2)+(3*Ll*H/16)* 
C0*(rl(j)-(l/2))*sin(3*tl(j)/2)+(-3*Ll/4) 

*CO*L* (rl ( j ) - ( 1/2) ) *sin(tl ( j ) /2) ; 
sll_3(j)=(3*as/2)*B0*(rs(j)~(-l/2))*cos(ts(j)/2)-(3*L2/8) 
*B0*(rs(j)~(l/2))*cos(3*ts(j)/2) ; 
sll_4(j)=(3*as/2)*D0*(rs(j)~(-l/2))*sin(ts(j)/2)-(3*L2/8) 
*D0*(rs(j)~(l/2))*sin(3*ts(j)/2) ; 
sll_5(j)=-3*Pl*A0*(rl(j)~(l/2))*cos(tl(j)/2)-(3*P2/2)*G 
*B0*(rs(j)~(l/2))*cos(ts(j)/2) ; 
sll_6(j)=3*Pl*C0*(rl(j)~(l/2))*sin(tl(j)/2)-(3*P2/2)*G 
*D0*(rs(j)~(l/2))*sin(ts(j)/2) ; 
sll_7(j)=2*Ql*qO*(r~(l/2))*sin(t(j)/2) ; 
sll_8( j )=-(3*k+2) *ac* (qO^r" (1/2) *sin(t (j )/2)+ql* 
r*cos(t  (j))+q2*r''  (3/2)*sin(3*t(j)/2)) ; 
sll (j )=(mue*exp(zeta* (TD+x_bar (j ))))*(sll_l(j)+sll_2(j) 

+sll_3(j)+sll_4(j)+sll_5(j)+sll_6(j)+sll_7(j)+sll_8(j)) ; 

s22_l(j)=(3*N/4)*A0*(rl(j)~(-l/2))*cos(tl(j)/2)-(3*Ll*N/16) 
*A0*(rl(j)~(l/2))*cos(3*tl(j)/2)-(3*Ll*P/4) 
*A0*(rl(j)~(l/2))*cos(tl(j)/2); 
s22_2(j)=-(3*N/4)*C0*(rl(j)~(-l/2))*sin(tl(j)/2)+(3*Ll*N/16) 
*C0*(rl(j)~(l/2))*sin(3*tl(j)/2)-(3*Ll*P/4) 
*C0*(rl(j)~(l/2))*sin(tl(j)/2) ; 
s22_3(j)=-(3*as/2)*B0*(rs(j)~(-l/2))*cos(ts(j)/2)+(3*L2/8) 
*B0*(rs(j)~(l/2))*cos(3*ts(j)/2) ; 
s22_4(j)=-(3*as/2)*D0*(rs(j)''(-l/2))*sin(ts(j)/2)  +  (3*L2/8) 
*D0*(rs(j)~(l/2))*sin(3*ts(j)/2) ; 
s22_5(j)=3*Pl*A0*(rl(j)~(l/2))*cos(tl(j)/2)-(3*P2*0/2) 
*B0*(rs(j)''(l/2))*cos(ts(j)/2) ; 
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s22_6(j)=-3*Pl*C0*(rl(j)~(l/2))*sin(tl(j)/2)-(3*P2*0/2) 

*D0* (rs ( j ) ~ (1/2) ) *sin(ts ( j ) /2) ; 
s22_7 ( j ) =-2*Ql*qO* (r- (1/2) ) *sin(t ( j ) /2) ; 

s22_8 ( j ) =- (3*k+2) *ac* (qO^r" (1/2) *sin (t ( j ) /2) ) +ql*r*cos (t ( j ) ) 
+q2*r''(3/2)*sin(3*t(j)/2)) ; 

if  j==NOP  y.i  j==l 

s22_9 ( j )=-ac* (3*k+2) * ( (l+as''2) / (l-as~2) ) *qO*sqrt (r) ; 
elseif  j==l 

s22_9  ( j )  =ac*  (3*k+2)  *  ( ( l+as~2)  /  ( l-as''2) )  *qO*sqrt  (r) ; 

else 

s22_9(j)=0; 

end 

s22(j)=(mue*exp(zeta*(TD+x_bar(j))))*(s22_l(j)+s22_2(j)+s22_3(j) 
+s22_4(j)+s22_5(j)+s22_6(j)+s22_7(j)+s22_8(j)+s22_9(j)) ; 

sl2_l(j)=(3/2)*A0*al*(rl(j)~(-l/2))*sin(tl(j)/2) 

-(3*Ll*al/8)*A0*(rl(j)~(l/2))*sin(3*tl(j)/2) ; 
sl2_2(j)=(3/2)*C0*al*(rl(j)~(-l/2))*cos(tl(j)/2) 

-(3*Ll*al/8) *C0* (rl(j)"(l/2)) *cos(3*tl (j )/2) ; 
sl2_3(j)=(3/4)*(l+(as-2))*B0*(rs(j)~(-l/2))*sin(ts(j)/2)-(3*L2/(as*16)) 
*BO* (rs (j)~(l/2)) *sin(3*ts ( j )/2)+(3*L2/ (as*4) )*(l-(as~2)) 
*B0*(rs(j)~(l/2))*sin(ts(j)/2) ; 

sl2_4(j)=-(3/4)*(l+(as~2))*D0*(rs(j)~(-l/2))*cos(ts(j)/2)-(3*L2/(as*16)) 
*D0* (rs (j)~(l/2)) *cos (3*ts ( j )/2)-(3*L2/ (as*4) )*(l-(as~2)) 
*D0*(rs(j)~(l/2))*cos(ts(j)/2) ; 
sl2_5(j)=2*Ql*qO*(r~(l/2))*cos(t(j)/2) ; 
sl2(j)=(mue*exp(zeta*(TD+x_bar(j))))*(sl2_l(j)+sl2_2(j) 
+sl2_3(j)+sl2_4(j)+sl2_5(j)) ; 

sed( j )=(l/(4* (mue*exp(zeta*x_bar ( j ) ) ) ) )* ( (1-nue) * ((  sll(j))~2 
+(s22(j))~2)-2*nue*  sll(j)*  s22(j)+2*(  sl2(j))~2); 
sigv(j)=sll(j)+s22(j)  ; 

sigr(j)=sll(j)*(cos(t(j)))~2+s22(j)*(sin(t(j)))~2 
+2*sl2(j)*(cos(t(j))*sin(t(j))) ; 
sigt(j)=sll(j)*(sin(t(j)))"2+s22(j)*(cos(t(j)))''2 
-2*sl2(j)*(cos(t(j))*sin(t(j))) ; 

sigrt(j)=-sll(j)*(cos(t(j))*sin(t(j)))+s22(j)*(sin(t(j))*cos(t(j))) 
+sl2(j)*((cos(t(j)))~2-(sin(t(j)))~2) ; 
siginl(j)  =  ((sll(j)+s22(j))/2)+sqrt((((sll(j)-s22(j))/2)~2)+sl2(j)"2)  ; 
tmax( j )=sqrt ( ( (s22( j)-sll (j ) ) /2) ~2+sl2(j ) “2) ; 
sedr ( j )=(!/ ( (1+nue) *4* (mue*exp(x_bar (j )))))*( (sigr ( j ) ) ~2 
+(sigt (j ) ) ~2-2*nue*sigt ( j ) *sigr ( j ) -nue~2* (sigr ( j ) 
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+sigt(j))+2*(l+nue)*(sigrt(j))~2) ; 


'/o  To  calculate  the  stress  filed  due  to  opposit  traction  at  the  crack  face 

for  i=l:l:N0P 

%  x(i)=xc(i) ; 

c(i)=i*0. 00006;  "/oO. 00006;  "/opoints  where  the  load  is  taken 

y.  p(i)  =  (5.5el8*c(i)~6-5.6el6*c(i)''5+2.3el4*c(i)''4-5ell*c(i)~3+ 
y7e8*c(i)~2+  6.8e4*c(i))/22;yc=100,  5  order  poly 
p(i)=0;  yc=0 

7.  p(i)  =  (1.2el3*c(i)~4-lell*c(i)~3+4.25e8*c(i)''2+1.3e5*c(i))/10; 
yforce  distribution  c=100  ,  3  ord  poly 
if  t(j)  <=  0 

y  p(i) =- (5. 5el8*c(i)~6-5.6el6*c(i)~5+2.3el4*c(i) "4-5011*0 (i) -3 
+7e8*c(i) "2+6 . 8e4*c(i) )/12; yc=100 ,  5  order  poly 
p(i)=-(3el3*c(i)"4-2.53ell*c(i)"3+1.05e9*c(i)"2+3.3e5*c(i))/12; 
yforce  distribution  c=250 

else 

y  p(i)=(5.5el8*c(i)"6-5.6el6*c(i)"5+2.3el4*c(i)"4-5ell*c(i)"3 
+7e8*c(i) "2+6 . 8e4*c(i) )/12;  yc=100,  5  order  poly 
p(i)=(3el3*c(i)"4-2.53ell*c(i)"3+1.05e9*c(i)"2+3.3e5*c(i))/12; 
yforce  distribution  c=250 

end; 

yp(i)=-(le25*c(i)"8-1.3143e23*c(i)"7+7.162e20*c(i)"6-2el8*c(i)"5 

-3 . 5el5*c (i) "4-3 . 667el2*c (i) "3+2 . 6e9*c (i) "2+8 . 6e4*c (i) ) /740 ; 
yc=250,  7th  poly 

yp(i)=2*(3el3*c(i)"4-2.53ell*c(i)"3+1.05e9*c(i)"2+3.3e5*c(i))/10; 
yforce  distribution  c=500 

'/•y(j)=yc(j) ; 
yc(j)=j *0.00006; 

r2(j ,i)=sqrt((x(j)+c(i))"2+y(j)"2) ; 
t2(j ,i)=atan2(y(j) , (x(j)+c(i))) ; 
pc=(p(i)*sqrt(c(i))) ; 
y  pc=p*sqrt(b) ; 

y  sigy(j)=(pc/(pi*rl(j) .*sqrt(r)))*(cos(tl(j)+t(j)/2) 
+(l/2)*sin(t(j))*sin(tl(j)+3*t(j)/2) 
+(r/rl(j))*sin(t(j))*sin(2*tl(j)+t(j)/2)) ; 
sigxCj ,i)=(pc/(pi*r2(j ,i) . *sqrt (r) ) ) * (cos (t2( j , i)+t ( j ) /2) 
-(l/2)*sin(t(j))*sin(t2(j,i)+3*t(j)/2) 
-(r/r2(j,i))*sin(t(j))*sin(2*t2(j ,i)+t(j)/2)) ; 
sigyCj ,i)=(pc/(pi*r2(j ,i) . *sqrt (r) ) ) * (cos (t2( j , i)+t ( j ) /2) 

+  (l/2)*sin(t(j))*sin(t2(j4)+3*t(j)/2) 

+(r/r2(j ,i))*sin(t(j))*sin(2*t2(j ,i)+t(j)/2)) ; 
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sigxyCj ,i)=(pc/(pi*r2(j ,i) .*sqrt(r)))*((l/2)*sin(t2(j ,i))*cos(t2(j ,i) 
+3*t(j)/2)+(r/r2(j ,i))*sin(t(j))*cos(2*t2(j ,i)+t(j)/2)) ; 

end 

end 

for  j=l:l:N0P 
sumx=0 ; 
sumy=0 ; 
sumxy=0 ; 
for  i=l:l:N0P 

sumx=sigx(j , i)+sumx; 
sumy=sigy(j , i)+sumy; 
sumxy=sigxy( j ,i)+sumxy; 

end 

sigll(j)=sumx; 

sig22(j)=sumy; 

sigl2(j)=sumxy; 

sedb ( j ) = ( 1/ (4* (mue*exp (x_bar (j)))))*( ( 1-nue) * ( (  sigl 1 ( j ) ) ~2 

+(sig22(j))~2)-2*nue*  sigll(j)*  sig22( j )+2* (sigl2(j ) ) "2) ; 
sigvb(j)=sigll(j)+sig22(j) ; 

sigrb(j)=sigll(j)*(cos(t(j)))''2+sig22(j)*(sin(t(j)))''2 
+2*sigl2(j)*(cos(t(j))*sin(t(j))) ; 
sigtb ( j )  =sigll  ( j )  *  (sin(t  ( j ) ) )  ''2+sig22  ( j )  *  (cos  (t  ( j ) ) )  "2 
-2*sigl2(j)*(cos(t(j))*sin(t(j))) ; 
sigrtb(j)=-sigll(j)*(cos(t(j))*sin(t (j))) 

+sig22(j)*(sin(t(j))*cos(t(j)))+sigl2(j)* 
((cos(t(j)))~2-(sin(t(j)))~2) ; 
siginlb(j)  =  ((sigll(j)+sig22(j))/2)+sqrt((((sigll(j) 
-sig22(j))/2)~2)+sigl2(j)~2) ; 
tmaxb(j)=sqrt(((sig22(j)-sigll(j))/2)''2+sigl2(j)~2)  ; 
sedrb( j )=(!/ ( (l+nue)*4* (mue*exp(x_bar ( j) ) ) ) ) * ( (sigrb(j ) ) ~2 
+  (sigtb(j))''2-2*nue*sigtb(j)*sigrb(j) 
-nue''2*(sigrb(j)+sigtb(j))+2*(l+nue)*(sigrtb(j))"2)  ; 

end 

for  j=l:l:N0P 

sigmal(j)=(((sll(j)+sigll(j))+(s22(j)+sig22(j)))/2)+sqrt(((((sll(j) 
+sigll(j))-(s22(j)+sig22(j)))/2)~2)+(sl2(j)+sigl2(j))-2) ; 
taumax(j)=sqrt((((s22(j)+sig22(j))-(sll(j)+sigll(j)))/2)~2 
+(sl2(j)+sigl2(j))~2); 

sigmat(j)=(sll(j)+sigll(j))*(sin(t(j)))~2+(s22(j)+sig22(j))* 

(cos(t(j)))''2-2*(sl2(j)+sigl2(j))*(cos(t(j))*sin(t(j))) ; 
sigmart(j)=-(sll(j)+sigll(j))*(cos(t (j))*sin(t(j)))+(s22(j) 

+sig22(j))*(sin(t(j))*cos(t(j)))+(sl2(j)+sigl2(j))* 
((cos(t(j)))''2-(sin(t(j)))''2) ; 
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sigmar(j)  =  (sll(j)+sigll(j))*(cos(t(j)))''2+(s22(j)+sig22(j))* 

(sin(t(j)))~2+2*(sl2(j)+sigl2(j))*(cos(t(j))*sin(t(j))) ; 

end 


plot(t*(180/pi) , (s22/(Kef/sqrt(2*pi*r))) , 'moO ;  %  combined 
xlabeK ’\theta  (degrees)’); 
ylabeK’u  (N.  mm  /mm"3)’); 
hold  on; 

plot (t* (180/pi) ,( (taumax)/(Kef /sqrt (2*pi*r) )), ’b’ ) ;  %  raction 

xlabeK ’\theta  (degree)’); 

ylabeK ’\sigma_y_y/  (K_e_f_f  /  \surd(2\pir) ) ’ ) ; 
hold  on; 

plot (t* (180/pi) , (sig22) /(Kef /sqrt (2*pi*r) ), ’r’ )  %  before  raction 
hold  on; 

xlabeK ’\theta  (degrees)’); 

ylabeK  ’  \tau_m_a_x  /  (K_l  /  \surd(2\pir) )  ’ ) ; 

title(’  \sigma_y’); 

legend(’q_o  =0’,  ’q_o  =  1000’,’q_o=  2000’); 
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APPENDIX  C 


Theoretical  Formulation 

Thermo-mechanical  Displacement,  Strain  and  Stress  Field  for  Mixed 
Mode  Loading  in  FGM 

The  equations  of  motion  for  a  plane  problem  are  given  by  Eq.  C.l 

daxx  drxY  _  d^u  dayy  drxy  _  d'^v 

dX  ^  dY  ~^dH'  dY  ^  dX  ~^dH  ^  ■  ’ 

The  relationship  between  stresses  and  strains  for  a  plane  strain  thermo-mechanical 


problem  can  be  written  as 

<Jxx  =  (A  -|-  2/i)exx  +  Aeyy  —  q;(3A  -|-  2/i)T  (C.2) 

<Jyy  =  (A  -|-  2/i)eyy  -|-  X^xx  ~  o(3A  -|-  2/i)T  (C'.3) 

Txy  =  p.exy  (C.4) 

where  X  and  Y  are  reference  coordinates, cTjj  and  where  i  =  X,Y  and  j  =  X,Y 


are  in-plane  stress  and  strain  components,  A  and  /i  denote  Lame’s  constant  and 
shear  modulus  respectively,  a  is  coefficient  of  thermal  expansion  and  T  represents 
the  change  in  temperature  in  the  infinite  medium.  For  non  homogenious  materials 
Shear  modulus  (/x).  Lame’s  constant  (A),  density  (p)  ,  thermal  expansion  (a)  and 
heat  conductivity  {k)  of  the  FGM  are  assumed  to  vary  in  an  exponential  manner 
as  given  by  Eq.  C.5,  whereas,  Poisson’s  ratio  (z/)  is  assumed  to  be  a  constant. 

/i  =  A  =  \Qe^^^\p  =  pQe^^^\  a  =  aQe^^^\  k  =  (C.5) 

subscript  ”o”  means  at  X  =  0  as  shown  in  Fig.  C.l.  and  C  and  (3  are  nonhomogene¬ 
ity  constants  that  have  the  dimension  {length)~^ .  For  plane  strain  deformation. 
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Figure  C.l.  Propagating  crack  tip  orientation  with  respect  to  reference  coordinate 
configuration. 

the  displacements  u  and  v  are  derived  from  dilatational  and  shear  wave  potentials 
0  and  0.  These  potentials  can  be  expressed  as 

(90  (90  (90  (90 

The  strain  can  be  derived  from  the  displacement  relation  giveb  by  Eq.  C.6 

du  (9^0  (9^0 

^  ^  dXdY 
dv  (9^0  (9^0 

~  dXdY 

du  dv  ^  d'^<p  (9^0  (9^0 

=  ^  +  ^  =  ^QxdY  ~  ^  ^  ^  ^  ^ 

Substituting  the  relation  for  /i,  A,  p,  o;  and  k  from  Eq.  C.5  the  Hookes’s  law  for 
non-homogenous  material  can  be  written  as 


o'xx  —  +  2fio)exx  +  Aoeyy  —  (3Ao  +  2fj,o)aoe^^^^T)  (C-8) 


o'YY  —  e^^^^((Ao  +  2po)eyy  +  Aoe^x  —  (3Ao  +  2po)ttoe*'^^^T)  (C-9) 

TxY  =  (C.IO) 
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By  substituting  the  strain  relation  given  by  Eq.  C.7  into  Eqs.  C.8  -  C.IO,  the 
Hookes’s  law  for  non-homogenious  material  in  terms  of  displacement  potentials 
can  be  written  as 


—  (3Ao  +  2/io)Q;oe^^'^^T) 


(C.ll) 


ctyy  e  ((  0+  ho)^y2  dXdY^^'^dX^^  dXdY 
—  (3Ao  +  2po)Q^oe^^^^T) 


txy  = 


/  d‘^(j)  d‘^ip  d‘^ip\ 

yww 


(C.12) 


(C.13) 


Substituting  the  stresses  from  Eqs.  C.ll  -  C.13  into  Eq.  C.l  and  after  simplihca- 
tion,  the  equations  of  motion  become 


d 


d 


(2/io  +  V  0  +  V  0  ~  ttot  (3Ao  +  2/io)  (  PT  + 


dX 


'BY 


dT\ 


C  0  +  +  ‘^QxdY^^  ~  aoe^^(3Ao  +  2fio)T] 


J 


(9^  /  (90  (90 


Q  Q  dT 

(2po  +  Ao)^  0  -  0  -  aoe^^(3Ao  +  2po)^ 


(C.14) 


^  (C.15) 


0o  0^0  0^0  0^0  \  (9^  /  (90  dip 

The  above  equation  can  be  further  simplihed  by  rearranging  Eqs.  C.14  and  C.15 
in  the  following  way 

^  {Eq.  C.14)  +  ^  {Eq.  C.15)  after  simplification  gives 

/  g  d  \ 

(2po  +  Ao)  V^0  -  C  (  (2po  +  Ao)^  0  +  0  ) 


—  Q;oe^'^(3Ao  +  2po)  (  (  (0  +  OPT  +  (20  +  C)77T7  +  \/'^T  )  —  PoT^  V 


dT 


dX 


02 


’dH 


(C.16) 
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Similarly  ^  {Eq.C.lA)  +  ^  {Eq.C.lh)  after  simplification  gives 


/io  VV  -  C  (  (/Wo^  ^ 


dT  (9^ 

Q;oe^'^(3Ao  +  2/io)C^^  =  t/’ 


(C.17) 


The  above  equation  can  be  further  simplihed  to 

(5  +  2)  0  -  C  +  2)^0  + 

—  ac(3Ao  +  2/io)  ^^(0  +  C)0^  ^T+(20  +  C)A 

(C.18) 


vV-C 


(90 

0X 


,  00 


ttc(3Ao  +  2/io)CA  ^ 


(9T 

0F 


Po 

/io  dH 


(C.19) 


where 

A-‘  =  0, 

V  Po 

Qic  =  Qioe^^  is  the  coefficient  of  thermal  expansion  in  the  vicinity  of  the  instanta¬ 
neous  crack  tip  and  is  assumed  to  be  constant. 


Transforming  to  the  Crack  Tip  Coordinates 

For  a  propagating  crack  shown  in  Fig.  C.l,  the  transformed  crack  tip  coordi¬ 
nates  can  be  written  as  x  =  X  —  ct  where  c  is  constant  crack  tip  speed.  In  the 
moving  coordinate  systems  the  above  Eqs.  C.18  and  C.19  can  be  written  as 


a. 


dx"^  dy"^  dx^  5  +  2  dy 


—  o;, 


.|0|(((;3  +  C)(3A-r  +  (2/3  +  OA-‘§)  + 


=  0 


(9a;2  dy"^  dx  dy 


0^(35  +  2)CA-^ 


dT 

dy 


0 


(C.20) 

(C.21) 
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where 


Q  =  Csa/2(1  —  v)/(1  —  2z/)  for  plane  strain  and  q  =  (1  —  v)  for  plane  stress. 

Q  and  Cg  are  the  elastic  dilatational  wave  speed  and  the  elastic  shear  wave  speed 
of  the  material  at  the  crack  tip. 

Asymptotic  Expansion  of  Crack  Tip  Stress  Fields 

Asymptotic  approach  is  used  in  deriving  solutions  for  displacement  potentials 
of  equations  of  motion  Eqs.  C.20  and  C.21.  A  new  set  of  coordinates  is  introduce 
as 

r/i  =  -,  =  -  (C.22) 


d  _  1  1  d  _l  I 

dx  e  drji  ’  dy  e  dri2  ’ 

where  e  is  an  arbitrary  parameter  and  is  assumed  to  be  0  <  e  <  1 
Equation  C.20  and  C.21  can  be  written  in  a  new  scaled  coordinates  (pi,  '/72)as 


dx"^ 


1  1 


and 


dy"^ 


1  1 

dr]\ 


(C.23) 


2^  /  (90  C 

drjf  drjl  ^  \  <9^1  5  +  2  dr]2 ) 

a=+0+(€"r  +  £q2/3  +  OA-‘|l  +  £‘‘(^(  +  O/3A-‘T)  =0  (C.24) 


2(9^0  (9^0 

drji  ^  dril 


.dip 
^  dry 


0^(35 +  2)e^CA“^ 


(9T 

drj2 


0 


(C.25) 


At  this  stage  it  is  assumed  that  0,  0,  and  T  can  be  represented  as  a  power  series 
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expansion  in  e. 


E(m  +  d ) 

e^^(pm{vi,V2) 

m=0 

oo 

ij{x,y)  =^{erii,eri2)  =  ^  1/2) 

m=0 

00 

E(m  +  1) 

e^^Tm{vi,V2) 


m=0 


(C.26) 


Substituting  Eq.  C.26  into  Eqs.  C.24  and  C.25  gives  the  following  equations 


^+3  f  2^  0m  0  0m  \  m+6  0u00m  C  00' 


E  m+j  I  2'^  'r  ill  I  ^  Y""!  1  I  '  "  /  /-^-r  iii  I 

m=0 


Clr 


^  drjl  '  dril 
(35  +  2)  /  m+5  m-|-7 


0?7i  5  +  2  dr]2 


(5  +  2) 


(20  +  C)  (0  +  C)  0A-X  j  )  =  0 

(C.27) 


m+3  f  20Vm  ,  0Vm\  ,  iri+5  0  00m  ,  ^.00 

Z^(e  ^  +  )  +e  ^  (C— +  C<5 


m=0 


077f  0775 

Q',(3i  +  2)£'^CA-'^)  =  0 
or]2 


dT]i 


dr]2 


(C.28) 


For  Eqs. C.27  and  C.28  to  be  valid,  the  partial  differential  equations  corresponding 
to  each  power  of  e  (e^A^g2^g5/2^  should  vanish  independently.  The  solution  for 
the  hrst  few  terms  of  the  series  are  considered. 

For  m  =  0,  the  above  equation  Eqs.  C.27  and  C.28  gives 

2  f  ,  2  ,  C  S'lpa' 

a:x  +  icx  +«MC— + 


drjf  dril 


drjl  5  +  2  dri2 


( etTo  +  ei  (20  +  C)  +  ei  (0  +  C)  )  =  0  (C.29) 


(5  +  2)  ^  --  '  --  dr], 

I  0  2  <9^0  ,  0Vo^  ,  +  0^000  ,  ^,000 


62  a 


+ 


+  +  -«e(35  +  2)6^CA 


^  dr]f  dr]2  )  \^  dr],  '  dr]2  J  '  “7"' 


_i0To 


=  0 


(C.30) 
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For  m  =  1,  the  above  equation  Eqs.  C.27  and  C.28  gives 


2  /  2  <9^1  ,  <9^1  ^  ,  3  [ ,  C  <9^1 


a. 


drjf  dri2 
{36  +  2) 


drji  5  +  2  dr]2 


(6  +  2) 


e^Ti  +  (2/3  +  C)  A 


_i^ 

dvi 


+  e^(/3  +  0/3A-^Tr  =0  (C.31) 


2  /  2 
€  a. 


ir^  +  -  ttc(35  +  2)e  CA  =  0  (C.32) 

oTj^  uTj2  j  \  <y^2  /  *y^2 


For  m  =  2,  the  above  equation  Eqs.  C.27  and  C.28  gives 


£2  a 


,(9^02  ,  (9^02 


+ 


a, 


dr/l  dr/2 
(35  +  2)  f  7 


+  eMC 


,<902 


+ 


C  002 


(9?7i  5  +  2  (9772 


(5  +  2) 


e^T2  +  ei  (20  +  C)  A'^l^  +  (0  +  C) 

dr/i 


=  0  (C.33) 


e2  a 


,(9^2  ,  0^2 


+ 


dr]{  dr/. 


0771 


-002 

dr/2 


-a035  +  2)eiCA-'^  =  O 


0772 


(C.34) 


By  collecting  the  partial  differential  equations  corresponding  to  each  power  of  e 
(e^A^g2^g5/2^  leads  to  the  following  set  of  partial  differential  equations. 

3 

6  2  terms(77r  =  0) 


6^  terms  (m  =  1) 


„2  0Vo  ,  0Vo  _  n 
dr/l  dr/l 

(C.35) 

20^0  0Vo_. 

dr/l  dr/l 

(C.36) 

2  0Vi  ,  0Vi  . 

dr/l  dr/l 

(C.37) 

2  0^1  0201  _  p 

dr/l  dr/l 

(C.38) 
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€2  terms  {m  =  2) 


2(9^2  <9^2  _  A *900  1  (9^o\  3(5  +  2 

(9?]^  (97^2  V  (5  +  2  (97^2  /  (5  +  2  ° 


(C.39) 


2(9^2  (9^2  ^  ^  ,  1  (900  \ 

(97]^  dri2  \  drji  <5  +  2  dri2  J 


(C.40) 


Eqs.  C.35  -  C.38  are  similar  to  that  for  homogeneous  material  where  as  the  par¬ 
tial  differential  equations  Eqs.  C.39  and  C.40,  associated  with  higher  powers  of 
are  coupled  to  the  differentials  of  the  lower  order  functions  through  the  nonho¬ 
mogeneity  parameters  and  temperature  term.  Eqs.  C.35  -  C.38  (i.e.  for  777  =  0 
and  777  =  1)  can  be  easily  reduced  to  Laplace’s  equation  in  the  respective  complex 
domains  0  =  hi  +  ^011^2,  C  =  hi  +  *Q^sh2  ;  i  =  and  the  solutions  are  same  as 
that  for  homogenous  material  and  can  be  written  as 

A  f  Q  +\  A  (m+3)/2  +  i  (m+3)/2  .  /"  (^  +  3)  ^  \  \ 

0m(PZ,  01,  t)  =  AmPl  COS  I  - ^ - 9i  1  CmPl  SlU  I  - ^ - 0i  1  (C.41) 


i^m{Ps,  Os,  t)  =  Sin  +  D,„p(-+3)/2 

(C.42) 

where 


Pi  =  {pI  +  a^plY^^  ,  tcin{0i)  =  ps  =  {pI  +  oIpIY'^  ,  tan(0,)  = 

hi  hi 

and  Am,  Bm,  Cm,  and  Dm  are  real  constants. 

The  solution  for  the  Eqs.  C.39  and  C.40  corresponding  to  higher  powers  of  e  (m 
=  2)  consists  of  two  parts:  solution  for  homogeneous  equation  and  a  particular 
solution  due  to  nonhomogeneity  and  temperature  and  these  can  be  obtained  re¬ 
cursively.  Eq.  C.39  can  be  divided  into  four  different  problems  as  shown  below 


aiP2 


2  I  „2_2\1/2 


asV2 
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2  0^2  <9^2  _  . 

dril  dril 

(C.43) 

2  <9^2  <9^2  _  ^<900 

drjf  dr]l  drji 

(C.44) 

2.9^2  ,  <9^2  C  <900 

r>;  -k  = 

(9?]^  drj2  (5  +  2  (97^2 

(C.45) 

„2'9V2  ,  <9^2  _  3<^  +  2^  ^ 

(9?7f  dr]l  (5  +  2  ° 

(C.46) 

By  solving  the  above  partial  differential  equations  separately  and  later  by  adding 
all  together,  the  solution  for  02  can  obtained  as 


02  =  A2Ppcos  +  C2Ppsin  (^di 


(C.47) 


Similarly  Eq.  C.40  can  be  divided  into  three  different  problems  as  shown  below 


a, 


,  (9^02  <9^0 


+ 


drji  dr]. 


=  0 


(C.48) 


2*9^2  <9^2  ^  ^<900 

drjf  dr]l  drji 

^2  <9^2  ,  <9^2  _  C  <900 
drjf  drj2  <5  +  2  dri2 


(C.49) 

(C.50) 


149 


By  solving  the  above  partial  differential  equations  separately  and  later  by  adding 
all  together,  the  solution  for  1^2  can  be  obtained  as 


^2  =  cos  +  D2pl^^  sin 

“  +/^ocos 

The  solutions  0o,  t/’o,  4>i  and  ^jJl  automatically  satisfy  the  compatibility  equations 
because  these  the  solutions  are  same  as  those  for  homogeneous  materials.  Since 
the  non- homogeneous  specific  parts  of  02  and  02  are  obtained  from  0o  and  ipo, 
they  also  automatically  satisfy  the  compatibility  equations. 

By  assembling  together  the  above  result  for  the  first  few  terms  and  by  transforming 
back  to  the  x  —  y  plane,  the  combined  solution  can  be  written  for  0  and  ip  as  Eq. 
C.52  and  C.53 


=  Aorf^cos  (  1  -h  Cor^^sin  ( 


3/2 


yl2r/  cos  [^2^1  j  +  C2r/  sin 
-  I^AqCos  +  C'osin 


c 


a. 


5  (5  +  2)  («2  -  a2) 
4  (35  +  2)  ttc 


)0  cos 


15  (5  +  2) 


+  Aivf  cos  {29i)  +  Cirf  sin  (26';) 


(C.52) 
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j  +  Birl  sin  {26 s)  +  Dir]  cos  (26'^) 

+  cos  +  D2rl/‘^  sin 

-  \  (^Bosm  +  Do  cos 

“  f  2^^  2\P^i^‘^  ( ^0 +  ^0 1  (C.53) 

5  (tt«  -  «s)  V  V2  /  V2  // 

where 


2„,2\l/2 


n  =  +  afy"^) 


2„,2\V2 


tan(6'i)  = - ,  Vs  =  +  a^y"^) 

X 


tan(0,)  =  ^ 

X 


The  above  displacement  potentials  are  now  used  with  Eq.C.6  to  get  the  dis¬ 
placements  helds.  The  in-plane  displacement  can  be  given  as 
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/3i/2  312  f  ^  a  \  j_0  f  ^  A 

V2  '  V2  7  8^2  '  \^2  7  V2  7/ / 

/3i/2.  1C  3/2/  o-  q\  \o  ■  aW\ 

V2  V2  7  Saf  '  \^2  7  V2  ))) 

+  f  cos  ( ^  — ry2  (^3  cos  r -  2  cos  r ^ 


f— -a^ry^cos  —  -— fssin  +2 sin  Dq 

V  2  "  V2  7  8  a,  "  V  V2  7  V2  7  /  / 

+  2Ain  cos  (6^;)  +  2Ciri  sin  {9i)  +  2i?ir,a,  cos  (26*,)  +  2Zi)ir,a,  sin  (26*,) 
+  \^2rf^  cos  f  +  ^52a,ry2  cos  f  ^ 


+  72ry  sin  -6*«  -  -D2asrl^‘^sm  -6*, 


AoV  cos  +  ^orf  sin  -6*z 


a2-a2  V  '  ^^^V2  7  ' 

72^7 

yTj  y„,f  ^os  -  c 

2°^"  2  COS  r  1^, 

(5  +  2af  -a2  V  V2 

3  (3(^  +  2) _ Oc _  3/2  .  3 

2  (5  +  2)  (a2-l)^°  ®  2 


sin  (  -6*, 


Aorf  cos  -6*z  -  Cory  sin  -6*z 


|^5ory2  cos  cos  00, 


(C.54) 
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V  = 


1  -  (  3sin  (  1  +  2sin  (  ^9, 


1  C  3/2 


+  ^-azr;/^cos  ^-^zj  -  ^3  cos  \^-9,  j  -2  cos  ^-^z 

+  (3 sin  -  2 sin 

+  (^-^’^y^cos  1^-3  cos  -  2  cos 

-  2Airiai  sin  {9i)  -  25ir^  sin  {9s)  +  2Ciairi  cos  {9i)  -  2Dir^  cos  (6^^ 

-  ^Asazr-y^sin  -  ^^arf^cos 

+  ^C'attzry^sin  -  ^Da^y^  cos 


8  az 
1  C  3/2 


Cn 


5n 


+ 


C^ai 
of  -  al 

c 


Aory^sin  (  )  +  Cq^^^cos  (  ^6*, 


3/2 


tt; 


5  +  2  q;2  _  0-2 

3(35  +  2)  ttc  3/2  -s 
2  7^K3iy^or/  cos-^ 


5ory2  sin  cos  ^ 


D, 


(C.55) 


Temperature  Fields  around  the  Crack  Tip 

In  this  analysis  it  is  assumed  that  the  temperature  held  around  the  crack 
tip  changes  asymptotically.  Also,  the  transient  effects  are  neglected.  The  heat 
conductivity  is  assumed  to  vary  exponentially  as  given  by  Eq.  C.5.  The  steady 
state  heat  conduction  equation  can  be  written  as 


^  (k^\  +  ^ 

dX\  dx)  BY 


(C.56) 


Where  k  is  the  coefficient  of  thermal  conductivity 

Assuming  that  k  is  constant  in  the  region  considered,  the  above  equation  can  be 
written  as 

BT 

V^T  +  /3— =0  (C.57) 

where 

Transforming  the  above  equation  to  the  crack-tip  moving  coordinate  system  {x  = 
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Xc  —  t,  y  =  Y),  Eq.  C.57  can  be  written  as 


v'T  +  /3— =  0 


(C.58) 


wtipre  V7^  —  4- 

wneie  v  —  dx^  ^  Oy^ 


As  this  stages  asymptotic  analysis  is  performed  to  solve  the  above  Eq.  C.58. 

Eq.  C.58  in  a  new  scaled  coordinates  {r]i,  772)  given  by  Eq.  C.22can  be  written 


v'T  +  e/3— =  0 


(C.59) 


where 

V  dri(  dri^ 


For  the  asymptotic  analysis  T  is  represented  as  a  power  series  expansion  in  e  as 


T{x,y)  =  T(r/ie,r/2e)  =  ^  r/2) 


(C.60) 


Substituting  Eq.C.60  into  Eq.C.59  gives  the  following  equation. 


±1)  .  d‘^T„ 


("1+3) 


+  62/3 


dril  drjl  )  Orji 


(C.61) 


For  Eq.C.61  to  be  valid,  the  partial  differential  equations  corresponding  to  each 
power  of  e  e,  e^/^...)  should  vanish  independently.  This  leads  to  the  set  of 

partial  differential  equations. 


For  m  =  0,  the  above  equation  Eq.  C.61  gives 


1  9,,Y^  'dr,. 


For  m  =  1,  the  above  equation  Eq.  C.61  gives 

^  \  drjl  drjl  )  ^  drji 


(C.62) 


(C.63) 


For  m  =  2,  the  above  equation  Eqs.  C.61  gives 


d^T2  d^T2\  .r2ndT2 
Y  +  Yfx)+  ^'13^  =  0 
I  dyl  )  dyi 


(C.64) 
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By  collecting  the  partial  differential  equations  corresponding  to  each  power  of  e 
...)  leads  to  the  following  set  of  partial  differential  equations. 
ternis(m  =  0) 


e  terms  (m  =  1) 


drjf  dril 


d^Ti  d^Ti 

dril  dr]l 


(C.65) 


(C.66) 


5 

62 


terms  {m 


2) 


fd^T2  d^T2\ 
V  dvj  ^  dT,l  ) 


+  /3|^  =  0 

orji 


(C.67) 


Eqs.  C.65  and  C.66  are  similar  to  that  for  homogeneous  material  where  as  the 
partial  differential  equation  Eq.  C.67  associated  with  higher  powers  of  m  is  cou¬ 
pled  to  the  differentials  of  the  lower  order  functions  through  the  nonhomogeneity 
parameters  and  temperature  term.  Eqs.  C.65  and  C.66  (i.e.  for  m  =  0  and  m  =  1) 
can  be  easily  reduced  to  Laplace’s  equation  in  the  respective  complex  domains 
p  =  Tji  ir]2  ,  and  i  =  and  the  solutions  are  same  as  that  for  homogenous 

material.  By  using  crack  tip  insulated  boundary  condition  i.e  ^  =  0,  at  6*  =  tt, 
the  solution  for  above  laplace  equation  can  written  as 
For  m  =  0 

To  =  qop^^'^  sin  (^9^  (C.68) 

For  m  =  1 

Ti  =  gipl/2 

where  p  =  {pf  + 

The  solution  T2,  for  Eq.C.67  corresponding  to  higher  powers  of  e  (m  =  2)  has  two 
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parts:  homogeneous  and  particular  solution  and  can  be  written  independently  as 
follows 


d^T2  d^T2  , 

dT]l  dri2 


(C.70) 


=  -P 


on 


(C.71) 


dr]l  drjl  J  drji 

Eq.  C.70  is  similar  to  that  for  homogeneous  material  and  the  particular  solution 
can  be  obtained  using  recursive  approach,  and  the  complete  solution  for  Eq.  C.70 
and  Eq.  C.71  is  given  below. 

Tj  =  sin  sin  (C.72) 

By  adding  the  above  terms,  the  temperature  held  for  the  hrst  2  terms  can  be  given 
as 

T  =  sin  +  gipcos  (P)  +  q2p^^‘^  sin  06^^  +  sin  Qp 

(C.73) 

where  p={pf  +  and  6=  tan“^  j 

Transforming  back  to  crack  tip  coordinates  {x  and  y),  the  temperature  held  near 
the  crack  tip  is  given  as 

T  =  qor^^'^  sin  +9'i^  cos  {9)+q2r^^‘^  sin  QoP'r^^'^  sin  (C.74) 

where  r={x‘^  +  and  6=  tan“^ 
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In-plane  Strain  Fields 

By  substituting  the  eqs.  C.54  and  C.55  in  to  Eq.  C.7,  and  after  simplification 
the  in  plane  strain  can  written  as 


+  I 


sin  — - ^r, '  6  sin  — - sin  — 

2/  Saf'V  v2/  2  V2 


'3  -1/2 

+  I  '  cos 

'3  -1/2  . 

+  I  !  sin  f  y 


3  C  1/2  ^3d 

—  —rj  cos 

lb  do 


3  C  1/2  • 

—  —rJ  sin 
lb  as 


2 


Bq  +  2yli 
-Do  +  ‘2Biaa 


+  cos  A2  +  y  cos  B2 

^  15  1/2  .  1/2  • 

+  — r/  sin  (  -  1^2  -  — r/  a,sin  (  -  IZI2 


_  3  ^1/2 


3  C 


q; 


2W 


lo  cos  I  -  1  -  Co  sin  I  - 


a. 


2  5  +  2  af  _  ^2 

3(5  +  2  (Uc 


y„cos(|)-H)„sm  (ly 


(5  +  2  a?  -  1 


sin 


Dn 


(C.75) 
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^yy  ~ 


3  _i 


8’  ‘  V 


--Vi  (  I  j  -  (t)  ^ 


3  -1/2  .  f  1  > 

-r;  sin  (  -  )  -  -Cr; 


'30 


-^a.r;V2cos  j^fJJ  Bo  +  2Aaf 

-5„,.;‘/^in  (ij  +  (^))  Bo  +  2B.a, 

- -afrl'^^  cos  f ^  ^2 - cos  f ^  -B2 

15  1/2  .  r0i\^  ^15  1/2  .  f0s]^ 

r,/  sin  I  -  1  ^2  +  — r/  a,  (  ^  )  ^2 


4«r' 


26 +  2af 
3(5  +  2  ac 


1/2, 


—  go  sin 


(C.76) 
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In-plane  Stress  Fields 

The  strain  fields  given  by  Eqs.  C.75  -  C.77  and  the  temperature  field  given 
by  Eq.  C.74  are  substituted  into  Eqs.  C.2  -  C.4  to  obtain  the  in-plan  stress  fields 
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around  the  crack  tip  presented  given  by  Eqs.  C.78  -  C.80 


G  -r 


exp{Cx)i^o 


—  ac{36  +  sin  +  qircos{9)  +  q2r^^‘^ 

+  ^go/5r='/"sin(^0) 


sin 


(C.78) 
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avY 


exp{Cx)fio 


-r-i  ^/^cos  l^jJ  (S-a^(S  + 

{ “S  (?')('*-  “?  ('5  + 


Vl6a?  ' 

/3  C  1/2 

— - ^r.  ( 


3  -1/2  . 

-r,  sm 

4  ' 


/3  C  1/2  .  /30z' 

\lQaf 


{5 -of  (5 +  2)))  Co 


(l)  ('^  +  «'('^  +  2)))  C'o 


/  <3  -1/2  I  \  3  C  1/2 

—  cis  COS  —  —  - — rj  cos 


o  -1/2  •  I  \  ^  1/2  •  / 

-r„  '  cisSin  —  —  - — rJ  sm  — — 


Sa^ 

3  C 


8  as 

4Ai  [S  -  af  (5  +  2))  -  45ia, 


^y^cos  (  ^  J  {S-a‘f{S  + 


'8ds 
2 
OOs 


B, 

a 


0 


^2  -  (  —rl^^as  cos 
C2  +  (  ^r]/‘^as  sin 


^25  +  2af  —  cig 
[85  +  2  ac 


6  +  2  af 


+  ac{36  +  2)  {qor^^'^  sin  +  qir  cos(6')  +  q2r^^‘^  '• 
+  sin 
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APPENDIX  D 


Split  Hopkinson  Bar  Theory 


Longitudinal  wave  in  a  homogeneous  bar  with  density  p  and  Young’s  modulus 
El,  can  be  given  as 

=  (D-1) 


where  c  is  longitudinal  wave  velocity  in  the  bar  given  by  c  =  m  is  the 


dE  dx"^ 

_  lEi 

VC  VCIUCIL^  Hi  LliC  Udi  g,lVCli  C  - 

displacement  and  t  is  time. 

The  general  solution  for  the  above  ID  partial  differential  equation  can  be  written 
as 

u  =  f{x  —  ct)  +  g{x  +  ct)  (D-2) 

where  /  and  g  are  arbitrary  functions  determined  by  initial  conditions. 

From  Eq.  D-2,  the  strain  (e),  stress  (a)  and  particle  velocity  (h)  in  the  bar 
can  be  derived  as 

dx  dx  dx 


a  =  Ee  =  E{f'  +  g' 


(D-4) 


dt  dt 


-cf  +  eg'  =  c{-f  +  g') 


(D-5) 


Split  Hopkinson  bar  apparatus  consists  of  two  elastic  bars  with  small  specimen 
attached  in  between  the  bars  as  shown  in  Fig.  D.l.  In  the  figure;  the  incident, 
reflected  and  transmitted  waves  are  indicated  with  the  corresponding  directions. 
Using  Eq.  D-2  the  displacements  of  the  bar  (m1  and  u2)  at  the  interfaces  can  be 
given  as, 

ui  =  fi{x  —  ct)  +  gi{x  +  ct)  =  Ui  +  Ur  (D-6) 
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Incident  Bar  Transmitter  Bar 


Figure  D.l.  Schematic  representation  of  SHPB 


U2  =  f2ix  -  ct)  =  Ut  (D-7) 

where  u*,  Ur  and  Ut  are  incident  reflected  and  transmitted  waves 
Following  similar  procedure  and  using  the  above  equations,  the  particle  velocity, 
strain  and  stress  at  the  interface  of  the  bar  are  derived  as  follows 
Particle  velocity 

Ml  =  c{-f[  +  g[)  =  c{-ei  +  er)  (D-8) 

U2  =  c(-/2)  =  -cet  (D-9) 

Strain 

(^1  =  fi+  g'l  =  Ui  + =  -ei  +  er  (D-10) 

62  =  /'  =  m'  =  et  (D-11) 

Stress 

cTi  =  Eei  =  E{ei  +  e^)  (D-12) 

a2  =  Ee2  =  Eet  (D-13) 
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So  long  as  the  stress  in  the  bars  remains  under  elastic  limits,  specimen  stress, 
strain,  and  strain  rate  can  be  calculated  from  the  recorded  strain  histories. 

Strain  rate  in  the  specimen  can  be  given  as 


Strain  in  the  specimen: 


Stress  in  the  specimen 


U2  -  Ul 


—  j  (ci  et  e^) 


6..  = 


igdi 


=  -  /  (e*  -  et  -  er)dt 

./n 


Fa 

dl. 


(D-14) 


(D-15) 


(D-16) 


where  Fa  average  force  of  two  interfaces,  and  A  and  Ag  are  bar  and  specimen  cross 
sectional  areas  respectively. 

The  average  force  Fa  can  be  given  as 


Fa 


Fi  +  F2 

2 

(7  lA  +  (72A 
2 


(D-17) 


Substituting  Eq.  D-17  in  to  Eq.  D-16  results  in 

(t^A  -j-  (T2A 

EA 

=  ^^(D  +  et  +  er)  (D-18) 

The  stress  and  strain  of  the  specimen  must  be  equilibrium  during  the  loading 
process,  i.e 


Ci  —  Ct  -f  er 


(D-19) 
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By  using  the  above  equilibrium  relation,  the  equation  for  strain  rate,  strain  and 
stress  in  the  specimen  can  be  reduce  to 


2c 

(D-20) 

2c  /■*  , 

=  -  j-  /  Crdt 

C  Jo 

(D-21) 

EA 

(D-22) 

=  .4, 

Conditions  to  Satisfy  the  Equilibrinm  of  the  Specimen 

•  The  elastic  waves  in  the  bars  must  be  one-dimension  longitudinal  wave  i.e 
the  wave  length  3>  radii  of  the  bar.  This  can  be  maintained  generally  by 
making  the  length  of  the  striker  bar  5  times  larger  than  the  radii  of  the  bar. 


•  The  specimen  must  deform  uniformly.  This  can  be  achieved  by 

a.  Making  the  wave  pulse  long  enough 

b.  Choosing  the  specimen  dimension  using  Eq.  D-24 


(D-23) 


where  Ig,  dg  and  Ug  are  length,  diameter  and  poison’s  ration  of  the  specimen, 

c.  Lubricating  the  interfaces  of  the  bars  and  the  specimen 


Dimension  of  the  Tungsten  Carbide  Inserts 

During  the  elevated  temperature  testing  two-carbide  inserts  are  placed  be¬ 
tween  the  two  bars  and  the  specimen  is  sandwiched  between  the  inserts  as  shown 
in  Fig.  D.2.  The  inserts  are  used  to  eliminate  the  temperature  gradient  in  the 
bars.  The  inserts  are  impedance  matched  to  the  bars  and  hence  do  not  distnrb  the 
incident,  transmitted  and  reflected  wave  prohles. 

The  dimension  of  the  tungsten  carbide  insert  is  chosen  by  using  the  impedance 
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Heater 


Figure  D.2.  Schematic  description  of  SHPB  with  heaters 


Table  D.l.  Material  properties  of  marraging  steel  and  tungsten  carbide 


Material  type 

Young’s  modulus  (GPa) 

Density  (p  ) 

Marraging  steel 

200 

8080 

Tungsten  carbide 

566 

143800 

relation  given  by  Eq.  D-25 

PsAsCs  =  PcAcCc  (D-24) 

where  ps,  Ag  and  c*  are  density,  cross  sectional  area  and  wave  speed  of  the  bar 
respectively  and  pc,  Ac  and  Cc  are  density,  cross  sectional  area  and  wave  speed  of 
the  inserts. 

The  mechanical  and  physical  properties  of  the  marraging  steel  and  the  tung¬ 
sten  carbide  is  given  in  the  following  Table.  D.l 

All  variables  in  Eq.  D-24,  except  the  the  cross-sectional  are  of  the  tungsten  car¬ 
bide  inserts  is  known.  The  cross-sectional  area  of  the  tungsten  carbide  insert  is 
determined  from  Eq.  D-24  as 

j  _ PsAgCg 

Pc^c 

=  56.6mm^  (D-25) 

i.e 

Dc  =  8.38mm  (D-26) 
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Therefore  a  tungsten  bar  with  diameter  of  -Dc=  8.38  mm  is  chosen.  To  make  sure 
that  the  impedance  is  properly  matched,  a  calibration  experiment  is  conducted. 
The  calibration  is  made  by  attaching  the  two  tungsten  carbide  inserts  and  sand¬ 
wiched  them  between  the  bars  as  shown  in  Fig.  D.3.  During  loading  the  projectile 
hits  the  incident  bar,  and  the  stress  wave  propagates  from  incident  bar  to  the 
transmitter  bar  through  the  tungsten  carbide  inserts.  The  strain  on  the  incident 
and  transmitter  bars  is  measured  using  the  stain  gages  (sg)  attached  on  the  bars. 

Fig.  D.4  shows  typical  incident  and  transmitted  signals  from  the  calibration  ex- 


Sg 

Sg 

Projectile  Incident  Bar  Inserts  Transmitter  Bar 

Figure  D.3.  Schematic  description  of  calibration  experimental  setup 

periment.  As  shown  in  the  hgure  the  incident  and  transmitted  signals  are  identical 
and  there  is  no  reflected  signal  observed.  This  can  be  seen  clearly  in  the  Fig.  D.5, 
where  the  incident  signal  is  plotted  on  the  top  of  the  reflected  signal. 

Specimen  Dimension 

The  length  of  the  specimen  is  determined  from  the  available  plate  thickness 
Length  of  the  specimen  is  3.175  mm 

To  determine  the  diameter  of  the  specimen  the  criteria  relation  given  be  Eq.  D-23 
is  considered  and  a  10%  axial  strain  is  assumed 
The  axial  strain  in  the  specimen  can  be  given  as 

AL 

ea  =  — =  0.1  (D-27) 
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Figure  D.5.  Incident  and  transmitted  signals  plotted  at  the  same  position 
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The  poison’s  ration  for  the  material  considered  is  v  is  0.3 
Therefore  the  lateral  strain  on  the  specimen  cab  be  written  as 


where  Df  =  Diameter  of  the  specimen  after  deformation 
Do  =  Diameter  of  the  specimen  before  deformation 
Rearranging  Eq.  D-28  will  results  in 


Di 

Do 


1.03  ^ 


(D-28) 


(D-29) 


At  this  point  it  should  be  notice  that,  the  diameter  of  the  deformed  specimen 
should  be  less  than  the  diameter  of  the  inserts,  i.e  Dj  <  8.38  mm 
Using  the  relation  given  by  Eq.  D-29  and  considering  the  maximum  limit  of  the 
deformation,  a  specimen  with  initial  diameter  Do  =  7.62  mm  is  obtained. 
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